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ABSTRACT 



This thesis contains programs written in APL and 
documentation for performing certain nonpar ametr ic 
tests and computing nonparametr ic confidence intervals. 
These methods of inference are particularly useful in 
dealing with Department of Defense related problems as 
illustrated in the several military examples worked in 
Appendix C. The following nonparametr ic tests are 
considered: Sign Test, Wilcoxon Signed-rank Test, Mann- 
Whitney Test, Kruskal-Wailis Test, Kendall’s B, 
Spearman’s R, and Nonparametr ic Linear Regression. The 
tests are based on the exact distributions of the 
respective test statistics unless a large sample 
approximation is determined to provide at least a three 
decimal place accuracy. The software consists of two 
APL workspaces; one, which is designed for 
microcomputers (IBM PC’s or compatibles) and is menu- 
driven, and the other, without menus, is designed for 
the mainframe computer (IBM 3033) at the Naval 
Postgraduate School. 
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I . INTRODUCTION 



Although nonparametr ic procedures are powerful 
tools to the analyst, they are currently underused and 
often avoided by potential users. Perhaps one reason 
for this is the difficulty in generating the exact 
distributions of the test statistics, even for moderate 
sample sizes. Consequently, tables of these 
distributions are only available for very small sample 
sizes and normal theory based approximations must then 
be used. 

The purpose of this thesis is to make a variety of 
nonparametric procedures quick, easy and accurate to 
apply using menu driven computer programs in APL . 1 
These programs use enumeration, recursion, or 
combinatorial formulas to generate the exact null 
distribution of the various nonparametric test 
statistics. This allows hypothesis testing and 
confidence interval estimation to be based on exact 
distributions without the use of tables. For larger 
sample sizes. the normal, F, and T distributions are 



1 APL was chosen because it is an interactive 
language that is especially powerful at performing 
calculations dealing with rank order statistics and 
vector arithmetic. Menus are not included in the 
workspace designed for the mainframe. 
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used to approximate the distributions of the test 
statistics with three decimal place accuracy. 

Section II addresses workspace design issues, to 
include, workspace requirements and assumptions 
regarding its use. Section III discusses the methods 
used to assess the accuracy of different asymptotic 
approximations, and the sample sizes required for an 
approximation to yield three decimal place accuracy. 
Section IV gives background information and discusses 
programming methodology for nonpar ametr ic tests based 
on single and paired sample data. In Section V, 
nonpar ametr ic tests for two or more independent samples 
are considered. Section VI discusses nonpar ametr ic 
tests for association; and, Section VII deals with 
nonpar ametr ic simple linear regression. Section VIII 
recommends other nonpar ametr ic tests that may be added 
to the workspace and areas for further work. 

To show application of nonparametric statistical 
methods to Department of Defense problems, several 
military examples are worked in Appendix C. 
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II. WORKSPACE DESIGN ISSUES 



This section presents a brief overview of the 
design considerations used in developing the APL 
workspace for both the mainframe and microcomputer. 

A. EQUIPMENT AND SOFTWARE REQUIREMENTS 

The microcomputer must be an IBM PC or AT 
compatible, equipped with 512 kilobytes of RAM and the 
APL*?LUS/?C system software, release 3.0 or later, and 
IBM’s DOS, version 2.00 or later . 1 The 8087 math 
coprocessor chip is not required to run this software, 
but will increase the computational speed. 

B. KNOWLEDGE LEVEL OF THE USER 

The user is expected to have had some exposure to 
APL and a working knowledge of nonpar ametr ic 
statistics. Familiarity with microcomputers or the 
Naval Postgraduate School mainframe computer is 
assumed . 



1 The APL system software requires 144 kilobytes of 
RAM while the NONPAR workspace requires an additional 
1 90 kilobytes . 
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SELECTION OF TESTS 



C . 

The nonparametr ic tests chosen for this workspace 
are some of the more widely known, and are considered 
basic material for any nonparametr ic statistics course. 
More information about the tests can be found in any of 
the textbooks that are referred' to in this document. 



D. MENU DISPLAYS 

The microcomputer’s workspace is designed around 
the use of menus. This was accomplished using the 
software package PC TOOLS from STSC . These menus are 
designed to guide the user through the selection of the 
tests without an excessive amount of prompting. The 
main menu displays the choices available in the 
workspace, while the test menus give the background 
information and options available for each test. Help 
menus to provide additional information about the 
tests are also available. 



E. ORGANIZATION OF WORKSPACE DOCUMENTATION 

Separate documentation is included for the 
microcomputer’s and mainframe computer’s workspaces 
(see Appendices A and B, respectively). These 
appendices explain the organisation and operation of 
the workspaces. Appendix C, which provides example 
problems for each nonparametr ic test, is applicable to 
both workspaces. 
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Ill . GENERAL SAMPLE SIZE CONSIDERATIONS AND ASYMPTOTIC 
APPROXIMATIONS 



In this thesis, the term alpha value is used in 
a general sense, and refers to the probability of 
rejecting a true null hypothesis. The term ?-value 
refers to the probability that a test statistic will 
exceed (or not exceed in the lower-tailed test) the 
computed value, when the hypothesis being tested is 
:rue . 

For selected values, the exact cumulative 
distribution functions (C.D.F.) of the test statistics 
are compared with those obtained from normal based 
asymptotic approximations. The results of the compari- 
sons are used as a basis for assessing the accuracy of 
the approximations. In those cases where more than one 
asymptotic approximation has been suggested in the 
literature, the accuracy of each approximation is 
compared over a range of desired C.D.F. values and 
sample sizes. From the results, the most consistently 
accurate approximation, and the sample size for which 
chat approximacion provides ax least three decimal 
place accuracy is determined. 

Once the accuracy comparisons were completed for a 
specific nonparametric test, microcomputer capabilities 
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IV. TESTS FOR LOCATION BASED ON SINGLE AND PAIRED- 
SAMPLE DATA 



The tests assume that the data consists of a 
single set of independent observations or paired 
observat ions ( Xj_ , Y i ) , i=1,2,...,N, fr om a continuous 
distribution. For the single and paired-sample cases, 
the null hypotheses are concerned with the median of 
the X and the median of the differences Xj_- Y^, 
respectively. The tests considered are the Ordinary 
Sign Test and the Wilcoxon Signed-Rank Test. 

A. ORDINARY SIGN TEST 

The Sign test can be used to test various 
hypothesis about the population median (or the median 
of the population of differences). Confidence 
intervals for these parameters can also be constructed. 
As a final option, nonparametr ic confidence intervals 
for the quantiles of a continuous distribution are 
offered . 

1 . Computation of the Test Statistic 

For single-sample data, the test statistic X 
is computed as the number of observations Xi greater 
than the hypothesized median Mg. For the paired-sample 
case, K is the number of differences X^ - Y ^ that exceed 
Mg). All observations X^ (or X^- Yi) that are equal to 



15 



size 



decreased 



Mg are ignored and the sample 
accordingly. As long as the number of such ties is 
small relative to the size of the sample, the test 
results are not greatly affected. Gibbons [Ref. 2:pp. 
1 08] . 

2 . The Null and Asymptotic Distribution of K 

The null distribution of X is binomial 
with p = .5. In Table 1 , the exact values of the 
C.D.F are compared with the corresponding approximate 
values using a normal approximation with and without 
continuity correction. 

TABLE 1 . C.D.F. COMPARISONS FOR THE SIGN TEST 



PROBE X < 


kJ ; FOR 


SAMPLE SIZE EQUAL 


TO 


24. 






TEST STAT. VALUE 


5 


i 6 
1 


1 7 

1 


3 


"T' 

1 

i 


9 

i 


i 10 i 

L . .. . 1 


11 


EXACT C.D.F. 


.00331 


I 

I .01133 

l 


r 

1 .03136 

i 


.07579 


1 

1 

i 


.15373 1 

' 


I 1 

1 .27063 1 


.41941 


ERROR; NORMAL 


.00117 


1 

1 .00417 

i 


1 

1 .01134 
1 


.02456 


r- 

i 

i 


1 

.04339 1 

i 


t I 

! .06352 1 

i 


.07736 


ERROR; NORM. U/CC 


*.00063 


1 

1 ".00104 
| 


1 

1 *.00114 
1 


*.00073 


—+~ 

! 

i 


1 

.00001 1 


i r 

.00043 1 
1. 


. 00023 


PROBE X < 


k3; FOR : 


SAMPLE SIZE E0.UAL 


TO 


25. 






TEST STAT. VALUE 


I 6 


1 7 


1 3 I 


9 


I 

1 

I 


I 

10 1 
l 


i 

11 1 

... . J_ 


12 


EXACT C.D.F. 


1 

1 .00732 

i 


n 

1 . 02164 


t 1 

1 .05333 I 

i t 


.11476 


~r~ 

i 

1 


.21213 ! 

t 


I 

.34502 1 

- L. 


. 50000 


ERROR: NORMAL 


! . 0026 o 


! .00774 


l i 

1 .01735 1 


. 03400 


1 


.05352 1 


. j707? 


. 07926 


ERROR; NORM. V/CC 1 

i 


! '.00033 


1 *.00111 
i 


1 *.00092 1 

_i l 


*.00031 


I 

__L_ 


.00032 ! 

L 


.00044 1 

i. 


. 00000 
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As can be seen, for sample sizes greater than 
25, a normal approximation with continuity correction 
is accurate to at least three decimal places. 

5 . Hypothesis Testing 

P-values are computed for three basic 
hypotheses comparing the median of the population or 
the median of the population of differences M with a 
hypothesized median Mg. P-values are taken from the 

cumulative distribution of the binomial for the 

following tests of hypothesis. 

a. One-sided Tests 

(1) H0 : M = M 0 Versus HI: M < M 0 . The 

P-value equals Pr [X _<_ k] , where k is the computed 
value of the test statistic. 

(2) H0 : M = M 0 Versus HI: M > M 0 . The 

P-value equals Pr [K _> k] . 

b. Two-sided Test. 

(1) H0 : M = M 0 Versus HI: M * M 0 . The 

P-value equals twice the smaller value of a( 1 ) or 

a(2), but does not exceed the value one. 

For sample sizes greater than 25, a normal 
approximation with continuity correction is used. 

4 . Confidence Interval Estimation 

Confidence intervals for the population 
median are based on the ordered observations in the 
sample. For paired-sample data, confidence bounds are 
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obtained from the ordered differences of the pairs of 
data. A 100(1- a )% confidence interval is determined 
in the following manner. Let k be the number such that 
Pr [K k] <_ ( a /2). Then, the (k+1 )th and (N-k)th 
order statistics constitute the end points of the 
confidence interval. Gibbons [Ref. 1: pp.104]. 

For computing confidence intervals when 

sample size N is greater than 25, a normal 

approximation with continuity correction is used. 

Also included under this test is an option to 
generate nonpar ametr ic confidence intervals for any 
specified quantile given a sample size N from a 
continuous distribution. The end points of the 
intervals are sample order statistics. 

B. WILCOXON SIGNED-RANK TEST 

The signed-rank test requires the added assumption 
that the underlying distribution is symmetric. This 
test uses the ranks of the differences X^- Mg (or 
Yi~ Mg) together with the signs of these differences to 
determine the test statistic. Confidence intervals for 
the median can also be constructed. 

1 . Computation of the Test Statistic 

For single-sample data, the test statistic W 
is computed as follows. 
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Lex Zi 



| 1 if Xj.- M 0 > 0 

( 0 if Xj.- M 0 <_ 0 



N 

and let rj = rank(|Xi~ M 0 | ) . Then, W = ]TZj.rj.. 

i = 1 

For paired sample data, W is calculated in 
the same manner , except the differences to be ranked 
are the paired-differences minus the hypothesized 
median. Zero differences are ignored and the sample 
size is decreased accordingly. When ties occur between 
ranks, the average value of the ranks involved are 
assigned to the tied positions. It has been shown that 
a moderate number of ties and zero differences has 
little effect on the test results. 1 

2 . The Null and Asymptotic Distribution of W 

The exact null distribution of W is given by: 
Pr [W = w] = un(w)/ 2 n , w = 0 , 1 , 2 , . . . , N( N+1 ) /2 , where 
un( w ) is the number of ways to assign plus and minus 
signs to the first N integers such that the sum of the 
positive integers equals w. It can be shown (see 
Gibbons [Ref. 1:pp. 112]) that un(w) , for successive 

values of N, can be computed using the recursive 
relationship : 

un ( w ) = un- 1 ( w-N ) t un-i(w) 



1 For more information on the effects that zeros and 
tied ranks have on the Wilcoxon Signed-Rank Test, see 
Pratt and Gibbons [Ref. 3]. 
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Exact C.D.F. values were compared with those 
obtained using the following asymtotic approximations: 
student’s T with ( N— 1 ) degrees of freedom (T), 
student’s T with continuity correction (TC), normal 
(Z), normal with a continuity correction (ZC), the 
average of T and Z as suggested by Iman [Ref. 4], and 
the average of TC and ZC. 

As can be seen in Table 2 below, the average 
of TC and ZC gives the most consistently accurate 
results with three decimal place accuracy when the 
sample size exceeds 9. 



TABLE 2. C.D.F. COMPARISONS FOR THE WILCOXON 
SIGNED-RANK TEST 

PROBE y i w] ; FOR SAMPLE SIZE EQUAL 10 9. 



TEST STAT. VALUE 


r 

3 


5 


1 

6 1 

l 


i 1 

3 


: 1 

9 


1 1 

12 1 

................. .-..i 


i 

1 14 


EXACT 


C.D.F. 


.0097? 


.01953 


i 

.02734 1 
1 


.04383 


.06445 


i 

.12500 1 
\ 


.17969 


ERROR; 


NORMAL 


". 00067 


.00046 


1 

.00204 1 
1 


.00591 


.00958 


1 

.01324 1 

- . j 


.02272 


ERROR; 


norm, y/cc 


‘.00243 


‘.0024? 


1 

•.00167 1 

1 


.00023 


.00269 


1 

.00593 1 
1 


.00306 


ERROR; 


T DIST 


.00518 


.00608 


.00673 1 

l 


.00701 


.0081? 


I 

.00326 1 

i 


.00313 


ERROR; 


T U/CC 


.00355 


. 00276 


I 

.00233 1 

i 


.00012 


‘.00010 


1 

‘.00437 1 

l 


'.00725 


ERROR; 


AOE T/Z 


.00225 

1 


.0032? 


1 

.00433 1 
» 


. 00646 

> 


• O’jSSS 


.01325 1 


.01545 


ERROR; 


AVE TC /ZC 


l .00056 1 


.00014 1 


.00033 1 


.00017 1 


.00129 


1 

.00128 1 


.00041 
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TA3L3 2. (Continued) 



PROBC U L ui] ; FOR SAMPLE SIZE EQUAL TO 10. 



TEST 3 


TAT. VALUE 


5 

» 


7“ ” ““1 

7 

1 


1 

8 1 
i 


10 


T 

1 12 

L 


1 15 

I 


17 


EXACT 


C.D. ?. 


.0097? 


.01355 


.02441 1 

i 


.04199 


1 .06543 


1 .11621 
I 


.16113 


ERROR; 


NORMAL 


*.00115 


; 

.00023 

L 


.00099 1 

i 


! .00478 


1 .0033? 

i 


1 .01490 

l 


.01333 


ERROR; 


NORM. U/CC 


1 '.00270 


f 1 1 

1 '.00219- 


I 

'.00193 1 


.00043 


1 .00229 


1 .00558 


.00710 


ERROR; 


T D 1 3 T 


t .00399 


1 .00512 


.00525 1 


.00665 


1 .00661 


1 .00651 

1 


.00631 


ERROR; 


T U/CC 


1 .00249 

i 


t .00243 

t . . _ . . 


.00132 1 
! 


.00151 


1 ".00052 

! .. . . 


1 *.00376 
1 


'.00571 


I2PQ2; 


AVE T/Z 


t 

1 .00142 


1 .0026? 


1 

.00312 1 


.00571 


P 

1 .00749 


1 .01075 


.01234 


W 41 *% W f 


wi :c/:c 


! '.00010 
» 


i .00012 


*.00008 1 

i 


.00097 


1 .00039 

t 


1 .00091 

J L 


.00070 



3 . Hypothesis Testing 

P-values are computed for three basic 
hypotheses comparing the median of the population or 
the median of the population of differences M with a 
hypothesized median M 0 as shown below. 

a. One-sided Tests 

(1) H0: M = M 0 Versus HI: M < M 0 . The 

P-value equals Pr [W <_ w] , where w is the computed 
value of the test statistic W. 

(2) H0 : M = M 0 Versus HI: M > M 0 . The 

P-value equals Pr [W >_ w] . 

b. Two-sided Test 

(1) H0 : M = M 0 Versus HI: M f M 0 . The 

P-value equals twice the smaller value of a( 1 ) or 

a(2), but not exceeding the value one. 
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For sample sizes greater than 9, an average 
of the normal and student’s T approximations, each with 
continuity correction, is used. Computations of the P- 
value for each of the alternative hypotheses are: 

a. HI : M < M 0 

Let Pzc = Pr[Z_£(w+.5- / ^w] and 



let ?TC 



r T (N -n < |w ' ” wl ' - 5 




N 0^2 [ | w — /i v 1 ” • 5 ] ^ 


• 

U 1 

1 


N - 1 



where Z is standard normal, T(^_i ) has a student’s T 
distribution with (N-1) degrees of freedom, u w = 
N ( N+ 1 ) / 4 and cr w 2 = ( N( N+1 ) ( 2N+ 1 ) / 24- ) ) . Then, the 

P-value for the test is (Pzc + O ~ p TC))/2 if w is 
less than and ( Pzc + p TC)/ p > otherwise. The above 
formulas are obtained from those given by Iman [Ref. 4] 
after inclusion of a continuity correction. 

b . HI : M > M 0 

The P-value equals ((1 - Pzc) + p TC )/ 2 if 
w is less than and ((1 - Pzc) + ( 1 ” p TC))/ 2 » 
otherwise. The computation of Pzc anci - P TC is similar 

to the above except the sign of the continuity 
correction is changed. 

c . HI : M t Mq 

The P-value equals twice the smaller 
value of a or b above, but not exceeding the value one. 
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Confidence Interval Estimation 



4 . 

For single-sample data, the confidence 
interval for the population median is based on the 
ordered averages of all pairs of observations (Xj_ + 
Xi ) / 2 such that i < 1. A 100(1- a )% confi dence 

u — 

interval is determined in the following manner. Let w 
be the number such that Pr [W £ w] <_ ( a/2). Then, the 
(w+1 )th and (m-w)th order statistics, where m = 
N(N+1)/2 or the total number of paired-averages, 
constitute the end points of the confidence interval. 
A confidence interval for paired-sample data is 
computed in the same manner , except the end points are 
taken from the paired-averages of the differences Xj_- 
Yi- Gibbons [Ref. 1:pp. 114-118]. 

For computing confidence intervals when 
sample sizes are greater than 9, a normal approximation 
with continuity correction is used. 
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V. TESTS BASED ON TWO OR MORE SAMPLES 



The tests assume that the data consists of 
independent random samples from two or more continuous 
distributions. The general null hypothesis is that the 
samples are drawn from identical populations. The 
Mann-Whitney and Kruskal-Wallis tests are considered. 

A. MANN-WHITNEY TEST 

The Mann-Whitney test is based on the distribution 
of the test statistic U, which can be used to compare 
the equality of the population medians or variances 
for two samples. 1 The Mann-whitney test with a 
modified ranking scheme can be used to test for 
equality of variances if the population means or 
medians are assumed to be equal (Conover [Ref. 5:pp. 
229-230]). If the medians differ by a known amount, the 
data can be adjusted before applying the test. A 
confidence interval for the difference in the medians 
of the two populations can also be estimated. 

1 . Computation of the Test Statistic 

For the comparison of population medians, the 
test statistic U is computed from the combined ordered 

1 The test statistic U and the method used to 
compute it are taken from Gibbons [Ref. 1:pp. 140-141]. 
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arrangement of observations Xj_ and Yj , i = 1,2,...,N; 

j = Let ri = rank(Xi) in the combined 

N 

ordered sample and Rx = Z] r i* Then, 

i = 1 

U = R X - M(M+1 )/2 . 



For testing the equality of variances, the 
computation of U is similar except for the method of 
assigning ranks to the ordered sample. This method 
ranks the smallest value 1 , largest value 2, second 
largest value 3, second smallest value 4, and so on , 
by two’s, until the middle of the combined ordered 
sample is reached. For either test, tied ranks for the 
combined sample are assigned the average value of the 
ranks involved. A moderate number of ties has little 
effect on the test results . 

2 . The Null and Asymptotic Distribution of U 

The exact null distribution of U is 
determined using a recursion algorithm due to Harding 
[Ref. 6]. 



Exact C.D.F. values were compared with 
approximate values obtained from the following 
asymtotic distributions: student’s T with (n-2) degrees 
of freedom where n = N + M, the total number of 
observations in both samples (T), student’s T with 
continuity correction ( TC ) , normal (Z), normal with 
continuity correction ( ZC ) , the average of T and Z 
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( Iman [Ref. 7]), and the average of TC and ZC . The 
results for various sample sizes are given in Table 3. 



TABLE 3. C.D.F. COMPARISONS FOR THE MANN-WHITNEY 
TEST 



P20BCU < u] ; 708 SAMPLE 31 113 N I0.UAL TO 9 AND M EQUAL TO 3. 



TEST '3TAT. VALUE 


i 

1 14 


1 17 


—I r- 

1 13 1 

i i 


21 


1 23 


1 27 


1 

1 29 


EXACT C. 0.7. 


1 .00933 


1 .01999 


1 .02515 1 


.04696 


1 .06796- 


1 .12904 


1 .17005 


£8202; NORMAL 


1 *.00025 


1 .00100 

j. _ 


1 .00163 1 

I - [ 


.00441 


1 .00632 


1 .01243 
1 


1 .01511 
; - 


£2202; NORM. U/CC 


1 ".00146 


i *.00114 


1 *.00033 i 


.00026 


i .00130 


i .00354 


r 

1 .00436 


£2202; T DI37 


1 .00237 


1 .00337 


1 .00372 1 


.00464 


1 .00525 

r 


1 .00676 


1 .00779 


Z2B02; T U/CC 


1 .00119 


1 .00108 


1 .00095 1 

| L 


.00007 1 

.. . . - ! 


! *.00073 


1 *.00259 

i 


1 *.00323 

r 


£2202; AVI T/Z 


1 .00105 


1 .00219 


1 .00270 1 


.00453 1 


.00603 


1 

1 .00959 


1 .01145 


£2202; A VS TC/ZC 


1 *.00014 


1 *.00003 


1 .00004 1 

f ! 


.00017 1 

I 


.00026 


1 .00043 
1 


1 .00054 
[ 


P20BCU < 


u]; FOR 


SAMPLE SIZES N EQUAL 


TO 7 AND M EQUAL TO 12. 




T23T STAT. VALUE 


14 


i 

1 17 

! . 


1 19 1 

i i 


21 1 


24 


1 27 1 

i i 


1 30 


EXACT C.D.F. 


.00353 


1 

1 .01792 
> 


I 1 

1 .02732 1 

t t 


.04156 1 


.07111 1 

1 


! 1 

1 .11342 1 

I | 


.17012 


E2P.02; NORMAL 


*. 00045 : 


1 .00062 

i . . 


1 .00137 1 

t ... i 


.00359 1 


.00701 1 

i 


1 .01097 1 

[ . .. ] 


.01437 


ERROR; NORM. U/CC 


*.00152 


l ~ - 

1 -.00123 
! 


i 1 

r. 00079 r 

l ! 


.00003 1 


i 

.00154 1 

j 


! .00322 1 
; ( 


.00453 


ER208; T DIST 


.00199 


1 

1 .00292 

r 


1 1 

i .00356 1 

r i 


.00422 1 


1 

.00527 1 
. | 


! .00643 1 

i 


.00796 


ERROR; T U/CC 


. 00094 


i : — 

I .00092 

i 


1 1 

1 .00063 1 

! 1 


.00026 1 


1 

*.00066 1 

l 


! I 

*.00173 1 

1 


*.00252 


ERROR; A VS T/Z 


.00077 


i 

1 .00177 


i t — 

1 .00272 1 


.00391 1 


i 

.00614 1 


l 

.00370 1 


.01142 


Z22 02: -i l /i 7C/3C 


*. 00029 


: *.00013 


*.00005 1 


.00011 1 


.00044 i 

L 


.00072 1 

L 


. 00093 



26 



TABLE 3. (Continued) 



?£03C’J 1 uJ; FOR 'SAMPLE SHIS’ M EQUAL 70 5 AND M EQUAL 70 17. 



7IS7 37A7. VALUE 


1 13 


l 

1 15 

1 


1 1 

! 13 1 

» 


20 


' 1 1 

1 23 1 


27 


1 

1 30 


EXAC7 C.D.F. 


I .00963 


1 .01933 

j 


I .02316 I 

! 1 


.04245 


1 '.07013 1 

• . . i 


. 12440 


1 .17935 


22202; NORMAL 


1 '.00077 

? 


1 .00039 

! 


1 .00170 1 
? 1 


.00349 


r - 11 — ■ l 

I .00633 I 

' __ t 


.01210 


1 .01564 


22202; NORM. U/CC 


J 

( '. 00190 


1 ".00130 


1 ".00037 1 


.00007 


' 1 " ' ' i 

1 .00133 1 


. 00444 


I .00573 


ERROR; 7 0137 


1 .00130 


1 

1 .00222 


I .00299 1 
1 


.00390 


1 .00544 1 


.00774 


1 .00955 

i 


ERROR; 7 U/CC 


1 .00017 

f 


1 .00023 
1 


l l 

1 .00023 1 


.00021 


1 1 
i .00007 i 
f— - f 


'.00025 1 

1 


1 *.00051 


ERROR; AVI 7/1 


I 

1 .00027 


1 .00120 


1 r 

! .00233 1 • 


. 00363 


r ! 

! .00615 1 


j 

.00992 1 


1 .01260 


22202; AVI 70/ ZC 


! '. 00036 


1 ".00063 


1 ".00032 1 

L 1 


.00014 


1 .00097 1 

1 L 


.00209 1 


.00254 


PROBE U i 


uJ : FOR 


SAMPLE 31 


ZE3 N EQUAL 


. 70 3 


AND M EQUAL 70 27. 




7137 37A7. VALUE 


1 7 


i 10 


1 12 1 


15 


1 13 1 

' i 


23 1 


26 


EXAC7 C.D.7. 1 

i 


1 .00764 


1 .01650 
1 


1 .02512 1 


.04236 


1 .07734 1 

i t 


.12660 1 
1 


.17433 


J 

ERROR; NORMAL 1 


1 '.00265 


1 

1 '.00099 


r t 

1 .00072 1 

' ! 


.00333 


I \ 

1 .00374 1 

| ! 


l 

.01342 1 


.01630 


ERROR; NORM. U/CC I 

i 


'.00363 


1 '.00254 

i 


1 '.00133 1 

i i 


.00089 


[ I 

1 .00405 1 


.00665 1 

1 


.00331 


1 

ERROR; 7 DI37 1 

i 


'.00121 


1 

1 .00031 

i 


1 I 

1 .00173 1 

L 1 


.00414 


1 i 

1 .00741 1 

1 i 


.01026 1 

1 


.01253 


ERROR; 7 U/CC 1 


*.00219 


1 

1 *.00129 


r t 

1 *.00042 1 

i i 


.00097 


I 1 

1 .00250 1 

t ! 


1 

.00323 1 

L 


.00390 


ERROR; AVE T/Z 1 

i 


'.00133 


1 *.00034 

i 


1 .00122 1 

i | 


.00402 


1 t 

1 .00803 1 
: | 


I 

.01134 1 
1 


.01466 


i 

ERROR; AVE TC/ZC 1 

i 


'.00291 


1 

1 ".00192 


! '.00037 1 

! L_ 


.00093 


I .00327 1 

J L. 


r 

.00496 I 

L 


.00610 



As can be seen from the tables, the average 
of ZC and TC gives the most consistently accurate 
results. For sample sizes NxM > 80, nearly three 
decimal place accuracy is obtained in all cases. 
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Hypothesis Testing; 



3 . 

P-values are computed for three basic 
hypotheses comparing the medians or variances of the 
two populations as shown below. 



a. One-sided Tests 

( 1 ) H0 : Mx = My Versus H 1 : My < My or 



H0 : V X = Vy 


Versus 


HI : 


v x 


> Vy. 


The P-value 


equals 


Pr [U <_ u] , 


where u 


is the 


observed 


value of the 


test 


statistic . 
















(2) H0 : 


M X = 


My 


Versus 


HI : M X > 


My or 


H0 : V X = Vy 


Versus 


HI : 


v x 


< Vy. 


The P-value 


equals 


Pr [U >_ u] . 














b . 


Two-sided Test 










( 1 ) H0 : 


M X = 


My 


Versus 


HI : M x # 


My or 


H0 : V X = Vy 


Versus 


HI : 


v x 


# Vy. 


The P-value 


equals 



twice the smaller value of a( 1 ) or a(2), but not 
exceeding the value one. 

For sample sizes NxM greater than 80, the 
average of the normal and student’s T approximations, 
each with continuity correction, is used. Computations 
of the P-value for each alternative hypothesis are: 
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a. Hi : Mx < My or Vx > Vy 

Let Pzc s Pr [ Z ± ( u +.5 - Mu)/ a u3 and 



let P^c = P r 



T(n-2) i. 



|u - ^ u | - .5 

( N-f-M- 1 ) °u“ C|u-/i u !-.5] 2 



. 5 



N+M-2 



N+M-2 



where Z is standard normal, T( n _ 2 ) has a student’s T 
disoribuxion with (n-2) degrees of freedom, M u - NxM/2 
and a u 2 = ( N(M) ( N+M+1 ) )/1 2 . Then the P-value for the 
test is ( ?zc + O ~ ?tc))/ 2 for u less than Mu and 
(PZC + Ptc)/ 2, otherwise. The above formulas are 

obtained from those given by Iman [Ref. 7] after 
inclusion of the continuity correction. 

b. HI : My > My or Vy < Vy 

The P-value equals ((1 - Pzc) + Ptc)/ 2 if 
u is less than Mu and ((1 - Pzc) + (1 Ptc))/ 2 > 

otherwise. The computation of Pzc and- PTC is similar 

to the above except the sign of the continuity 
correction is changed. 

c . HI : Mx t My or Vx 4 Vy 

The P-value equals twice the smaller 
value of a or b above, but not exceeding the value one. 

4- . Confidence Interval Estimation 

Confidence intervals for the difference in 
medians, (My - Mx ) , are based on the ordered 
arrangement of the differences (Yj- j = 1,2,...,M; 
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i = 1 , 2 , . . . , N for all i and 3 . A 100(1 - a )% 
confidence interval is determined in the following 
manner . Let u be the number such that Pr [U u] 

( a /2). Then, the (u+1 )th and (m-u)th order 

statistics, where m = NxM or the total number of 
possible differences, constitute the end points of the 
confidence interval. 

For computing confidence intervals when 
sample sizes NxM are greater than 80, a normal 
approximation with continuity correction is used. 

3. KRUSKAL-WALLIS TEST 

The Kruskal -Wallis test is a nonpar ametr ic analog 
of the one-way classification analysis of variance test 
for equality of several population medians. Gibbons 
[Ref . 1 :pp . 199]. 

1 . Computation of the Test Statistic 

Calculations of the test statistic H center 
around the ordered arrangement of the combined samples 
from which the sum of ranks for each sample is derived. 
Let Xij , j= 1 , 2,..., ni and i= 1 , 2 ,...,k, be independent 

random samples from k populations. Let r^j = rank(Xj_j), 

nj_ k 

Rj_ = ^ r ji , an d N = 2 n i* Then , 

3=1 i = 1 

( k 

H = ( 12/(N(N+1 )) £ (RiS/m) 

(i = 1 



- 3(N+1 ) 
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If ties occur in the combined sample, they 
are resolved by assigning the average value of the 
ranks involved. A correction' based on the number of 
observations tied at a given rank and the number of 
ranks involved, is included in the calculations. A 
complete description of the correction factor is given 
in Gibbons [Ref. 2:pp. 178-179]. 

2 . The Null and Asymptotic distribution of H 

The null distribution of H is generated by 
enumeration. Each possible permutation of ranks is 
listed for the combined sample, and the corresponding H 
value computed. The frequency distribution of H.is the 
total number of occurrences of each distinct H value. 
The H values are arranged in increasing order while 
maintaining the frequency pairings. The null 
distribution is obtained by dividing the cumulative 
frequencies by n-j ! n 2 ! • • • ! /N ! . 

Due to computer limitations, generation of 
the exact distribution of H was only possible for k = 3 
populations with n = 4 observations in each, and 4 
populations with 3 observations in each. Most of the 
distributions were generated on the mainframe computer 
and saved in matrices for quick recall by the Xruskal- 
Wallis test program. 

Exact C.D.F. values were compared with the 
corresponding approximate values using the following 
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distributions: chi-square with (k-1) degrees of freedom 
(C), F distribution with (k-1 ) and (N-k) degrees of 
freedom (F), and F with (k-1 ) and (N-k-1 ) degrees of 
freedom (FI). The chi-square distribution uses the 
Kruskal-Waliis H statistic, while the F and FI 

distributions use a modified H statistic, HI = 
( (N-k)H)/(k-1 )( ( N— 1 )-H) ; see Iman and Davenport [Ref. 
3]. As can be seen in Table 4, FI gives the most 
consistently accurate estimates. 



TABLE <1. C.D.F. COMPARISONS FOR THE XRUSKAL- 
W ALL IS TEST 



PROBZH 5 hi; FOR .1 CROUP OF 3 SAMPLES CONSISTING OF , 4, AND 3 OBS. 



TEST ST AT. VALUE 


1 


7.1439 


1 5.7121 | 5.1313 


I 5.5935 


( 5.0536 


1 4.2121 


i 


3.5935 


EXACT C.D.F. 


i 


.10970 


1 .01905 I .02961 


| .04365 


( .07810 


1 .12918 


| 


.17784 


ERROR ; CRT SQUARE 


1 ".01340 


(". 01532 |". 01535 


(".01220 


(".00184 


I .00746 


i 


.01241 


ERROR ; F LIST 


1 


.00304 


I .00736 | .00836 


j .01113 


I .01820 


j .01696 


j 


.00990 


ERROR ; F Wt ~ 1 DF 


1 


.00084 


.00426 j .00403 


I .00544 


1 .01138 


I .00895 


i 


.00172 


PROBZH 2 hi i 


; FOR A GROUP OF 3 SAMPLES CONSISTING OF 4,4 


, AND 4 OBS . 




TEST ST AT. VALUE 


i 


7.6538 


j 6.9615 I 6.5000 


I 5.6923 


I 4.9615 


1 4.2692 


i 


3.5769 


EXACT C.D.F. 


i 


.00762 


I .01939 | .02996 


( .04366 


( .03000 


(’ .12190 


i 


.17299 


ERROR ; CEISQUARE 


r 


.01416 


(".01139 I". 00882 


(".00941 


(".00368 


I .00361 


i 


.00577 


ERROR ; F BIST 


i 


.00290 


j .00839 I .01204 


I .01100 


1 .01272 


j .01225 


i 


.00263 


ERROR ; F UTl DF 


i 


.00149 


I .00600 | .00890 


| .00647 


I .00708 


I .00592 


p 


.00384 


PROBZH Shi; 


FOR A GROUP OF 4 SAMPLES CONSISTING 


OF 3, 3, 


2, AND 2 


OBS . 


TEST ST AT. VALUE 


i 


7.6364 


1 7.1818 | 7.0000 


( 6.5273 


I 6.0132 


j 5.3818 


j 


4.8727 


EXACT C.D.F. 


i 


.01000 


1 .01921 I .02921 


J .04921 


( .07984 


( .12984 


i 


. 17952 


ERROR ; CEISQUARE 


r 


.04416 


I". 04712 |". 04269 


(".03939 


( ".03089 


(".01604 


p 


.00183 


ERROR; F DIST 


i 


.00234 


1 .00260 | .00730 


I .00330 


I .01034 


1 .01123 


i 


.00902 


ERROR; F U/~l DF 


1 


.00172 


~ . oougu I ~ . Q01U7 


I " .00307 


i“. 30375 


j 00570 


! " 


. 00843 


PROBZH 2 hi; 


FOR A GROUP OF 4 SAMPLES CONSISTING 


OF 3, 3, 


3 , AND 2 


OBS. 


TEST ST AT. VALUE 


i 


8.0152 


1 7.6364 j 7.1515 


I 6.7273 


1 6.1970 


I 5.4697 


j 


4.9697 


EXACT C.D.F. 


i 


.00961 


| .01831 I .02974 


j .04948 


( .07805 


I .12740 


i 


.17571 


ERROR ; CHI SQUARE 


r 


.03609 


(".03584 (".03748 


(".03164 


I " . 02436 


P . 01306 


i 


.00168 


ERROR ; F DIST 


i 


.00215 


j .00481 ( .00441 


I .00920 


( .01185 


I .01036 


i 


.01214 


ERROR-, F W/~ 1 DF 


r 


.00133 


|~. 00019 I " .00269 


1 .00030 


1 .00098 


1 " . 00230 


p 


.00091 
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A final accuracy comparison between the C and FI 
approximations was conducted by computer simulation for 



5 populations with 8 observations each. Initially, 
30,000 permutations of the 40 ranks were randomly 
generated (no tie ranks allowed) and the H statistic 
calculated for each permutation. Then the empirically 
determined percentiles Hp for selected values of p 
between .01 and .18 were compared with the 
approximations given by the C and FI distributions. 
The results are shown in Table 5. It can be seen that 
the FI approximation compares well with the simulated 
results, giving three decimal place accuracy, while 
the C approximation is less accurate. 



TABLE 5. C.D.F. COMPARISONS FOR THE KRUSKAL- 
W ALL IS TEST USING COMPUTER SIMULATION 



PROBiE 2 his BASED ON 10000 GENERATED B'S FOR 5 SAMPLES OF 8 OBS. EACH. 



TEST STAT. VALUE 


\ 12.229 


1 


11.065 


1 10.248 


j 


9.212 I 


8.129 


j 


7.030 


1 6.232 


C.D.F. 


VALUE 


1 .01000 


1 


.02000 


j .03000 


1 


.05000 | 


.08000 


i 


.13000 


j .18000 


ERROR ; 


C3ISQUARE 


(".00573 


r 


.00584 


|". 00646 


r 


.00601 j“ 


.00696 


r 


'.00432 


j". 00246 


ERROR; 


F WT 1 DF 


I .00081 


1 


.00231 


j .00259 


i 


.00350 j 


.00161 


i 


.00109 


(".00103 


PROBIE 2 h]; BASED ON 


20000 GENERATED E' 


S FOR 5 SAMPLES OF 


8 


OBS. EACH. 


TEST STAT. VALUE 


j 12.315 


1 


11.054 


| 10.184 


i 


9.163 | 


8 . 129 


j 


7.034 


1 6.220 


C.D.F. 


VALUE 


j .01000 


1 


.02000 


j .03000 


i 


.05000 j 


.08000 


i 


.13000 


j .18000 


ERROR ; 


CHISQUARE 


j~. 00516 


1" 


.00596 


I". 00745 


r 


.00716 |~ 


.00696 


i' 


’.00413 


j". 00334 


ERROR ; 


T V/~l DF 


1 .00126 


1 


.30221 


j .00186 


i 


.00235 j 


.00161 


i 


.00130 


j". 00199 


PRO. 


5ZH > h : ; 3ASED ON 


3 COCO GENERATED S'. 


S FOR 5 SAMPLES OF 


3 


OBS . EA 


.an. 


TEST STAT. VALUE 


j 12.305 


1 


10.976 


| 10.147 


j 


9.179 | 


8.168 


i 


7.072 


| 6.265 


C.D.F. 


VALUE 


1 .01000 


1 


.02000 


I .03000 


1 


.05000 I 


.08000 


i 


.13000 


| .18000 


ERROR ; 


CHISQUARE 


j". 00522 


r 


.00684 


(".00802 


i' 


.00677 j" 


.00563 


r 


.00213 


(".00020 


ERROR ; 


F WTl DF 


I .00121 


i 


.00143 


1 .00111 


i 


.00273 | 


.00301 


i 


.00344 


j .00142 
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3 . Hypothesis Testing 



p- 


values 


for 


the test 


H0 


: the popul 


at ion 


medians are 


: all 


equal 


versus 


HI : 


at least 


two 


population 


medians 


are 


not equal 


) 


are computed 


as : 


Pr [H _> h] , 


where 


h is 


the value 


of 


the observed 


test 



statistic . 

For three or more populations with at least 4 
observations in each, the FI approximation is used. 
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VI . TESTS FOR ASSOCIATION IN PAIRED-SAMPLES 



The tests described herein assume that the data 
consists of independent pairs of observations (Xj_,Yi) 
from a bivariate distribution. The general null 
hypothesis is that of no association between X and Y. 
Kendall’s B and Spearman’s R are considered. 

A . KENDALL ’ S 3 

1 . Computation of the Tesx Statistic 

The test statistic is computed by comparing 
each observation (Xj_,Yi ) with all other observations 
( X j ,Yj ) in the sample. If the changes in X and Y are 
of the sarnie sign, sgn(Xj - Xj.) = sgn(Yj - Yi), the pair 
(Xi,Yi) and (Xj.Yj) is "concordant" and a +1 is scored. 
If the signs are different, the pair is "discordant" 
and a -1 is scored. Any ties between either the X’s or 
the Y’s scores a zero for that pair. The sum of all 
scores divided by the total number of distinguishable 
pairs, (N(N-1))/2, gives B. If zeros are scored, the 
denominator is reduced by a correction factor which is 
based on *he number of observations tied ax a given 
rank and the number of ranks involved in each of the X 
and Y samples. A complete description of the 
correction for ties is given in Gibbons [Ref. 2:pp. 
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289] . The value of B ranges between 1 , indicating 
perfect concordance, and -1, for perfect discordance. 
Gibbons [Ref. 1:pp. 209-225]. 

2 . The Null and Asymptotic Distribution of B 

. The null distribution of B is derived from 
the following recursive formula given in Gibbons [Ref. 

1 : pp . 2 1 6 ] . 

u(N+1,P) = u(N,P) + u( N , P-1 ) + u( N , P-2 ) + ...+ u(N,P-N) 

where u(N,P) denotes the number of P concordant 
pairings of N ranks. This formula is used to generate 
the frequency with which the possible values of P 
occur. Division by N! results in the probability 
distribution of P. Since, B = ( 4P/ ( N( N- 1 ) ) ) - 1 , the 
null distribution of B is easily determined. 

Exact C.D.F. values were compared with those 
obtained using a normal approximation, with and without 
a continuity correction factor (CC = 6 /N(n 2-1), pro- 

posed by Pittman [Ref. 11] for the Spearman’s R test). 
The results for various sample sizes are provided in 
Table 6. As can be seen, for sample sizes greater than 
13, a normal approximation with continuity correction 
provides three decimal place accuracy. 
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TA3L2 5. 



C.D.F. COMPARISONS FOR KENDALL’S 3 





? 20 SC 3 2 


b]; 702 


SAMP LI 31 21 


EQUAL 


TO 


13. 






TEST 37AT. YAWS 


[ 

1 0.5123 


1 0.4315 
} 


i [ 

1 0.4103 1 

f ’ 


0.3530 


T" 

1 


; 

0.3333 1 


[ 

1 0.2564 ! 


! 0.2303 


EXACT C.O.F. 


1 .00743 
1 


1 .01524 
\ 


1 .02353 ! 


.04393 


I 

» 


.05443 1 

• i 


.12533 1 


! .15303 


222 OS; NORMAL 


1 .00014 
- 1 


1 .00121 


1 .00313 1 
[ « 


.00520 


1 

I 


.00303 1 
! 


.01473 1 
1 


! .01703 


ERROR: NORM. IJ/CC 


r 

1 '.00013 


1 .00073 


! .00233 1 

t 


.00437 


1 


.00538 : 
• 


.01223 1 
| 


.01415 




PR0BC3 2 


b] ; FOR 


SAMPLE SIZE 


EQUAL 


TO 


14. 






TIST S7AT. WLU5 


- T - 

1 0.4725 


■ r 

I 0.4233 


i i 

1 0.4053 1 


0.3525 


I 

l 

l 


I 

0.2367 1 


! 1 ‘1 
0.2527 1 


0.2033 


EXACT C .3.7. 


1 .00964 


1 .01773 


1 .02253 1 


.03373 


i 


.07353 i 

i 


.11535 I 


.16541 


ERROR; NORMAL 


i .00035 
1 


1 .00140 


; f 

1 .00213 1 

I » 


.00431 


i 

1 

i 


1 

.00333 1 


i 

.01256 1 


.01527 


ERROR; NORM. WCC 


! .00003 


! .00035 


1 .00152 ! 

i • 


.00345 


1 

1 


.00742 1 


.01057 1 


.01371 



3 . Hypothesis Testing 

P-values for tests of no association between 
X and Y are computed for three types of alternative 
hypotheses. Because the distribution of B is symmetric, 
all probabilities cam be taken from the upper tail 
using the absolute value of b, the observed value of 
the test statistic. Linear interpolation is used when 
b lies between tabulated values. The P-values are 
computed as follows. 

a. One-Sided Alternatives 

The one-sided alternative tested depends 
on the sign of b. A positive b will automatically test 
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for direct association or concordance, while a negative 
b will test for indirect association or discordance. 
The P-value equals Pr[B _> |b|]. 

b. Two-Sided Alternative 

The P-value equals tvice the probability 
computed for the one-sided hypothesis. 

For sample sizes greater than 12, a normal 
approximation with continuity correction is used. The 
approximate P-values are then: 



1 - Pr[ Z ± ( ( 1 b | 

is standard normal, CC is 
fi-fr = 0, and ^tq- = ( 4N + 
sided test and twice this 
test . 



- CC ) - M-fr)/ , where 'Z 

the continuity correction, 
10)/9N(N-1), for the one- 
P-value for the two-sided 



3. SPEARMAN’S R 

The Spearman’s R Test requires the added 
assumption that the underlying bivariate distribution 
is continuous. The test measures the degree of 
correspondence between rankings, instead of the actual 
variate values, and can be used as a measure of 
association between X and Y. Gibbons [Ref. 1:pp. 226]. 

1 . Computation of the Test Statistic 

The test statistic R is computed in the 
following manner. Let ri = rank(X^) and s^ = rank(Yi) 
and Dj_ = rj_ - Sj_. Then, 
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N 

6 V] Di^ 

R = I - i =J 

N(N2 - 1 ) 

where N is the size of the sample. If ties occur in X 
or Y, they are resolved by assigning the average value 
of the ranks involved. A correction factor, based on 
the number of observations tied at a given rank and the 
number of ranks involved, is included in the 
calculations. A complete description of the 

correction factor is given in Gibbons [Rex'. 2:pp. 279]. 
The value of R ranges between 1 , indicating perfect 
direct association, and -1, for perfect indirect 
association. Gibbons [Ref. 1:pp. 226-235]. 

2 . The Null and Asymptotic Distribution of R 

The null distribution of R for a given sample 
size N is generated by enumeration. The method, as 
presented in Kendall [Ref. 9], involves generation of 
an N by N array of all possible squared differences 
between any two paired ranks of X and Y. All N! 
permutations of N ranks are used to index values from 
the array. The sum of these indexed values for each 
permutation gives rise to N! sum of squared differences 
which are then converted to the R statistic. The 
frequency distribution of R is the total number of 
occurrences of each distinct value of R divided by N! . 
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Due to mainframe computer memory limitations 
in the APL environment, generation of the distribution 
of R was limited to sample sizes of 7 or less. 1 Using 
tables, provided by Gibbons [Ref. 2:pp. 417-418] to 
supplement computer computations, a numerical matrix, 
called PMATSP, was created .to store the cumulative 
distributions of R for sample sizes less than 1 1 . This 
matrix allows for quick recall of cumulative 
probabilities by the Spearman’s R Test program. 

Exact C.D.F. values were compared with those 
obtained using a student’s T approximation with (N-2) 
degrees of freedom (see Glasser and Winter [Ref. 10]). 
and a normal approximation. Both normal and T 
approximations were computed with and without a 
continutity correction factor, CC = 6/N(N2-i ) (Pittman 
[Ref. 11]). From the results presented in Table 7, the 
most consistently accurate approximation is given by 
the T distribution with a correction. 

3 . Hypothesis Testing 

P-values for tests of no association between 
X and Y can be computed for three types of alternative 
hypotheses. Because the distribution of R is symmetric, 
all probabilities are taken from the upper tail using 



1 The memory capacity of the mainframe computer in 
the APL environment is limited to 2.5 megabytes. 
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TABLE 7. 



C.D.F. COMPARISONS FOR SPEARMAN’S R 



• 


FS03C2 > 


rJ; FOR 


SAMPLE SIZE 


EQUAL TO 


9. 






TEST STAT. VALUE 


1 0.7333 
» 


i 

1 0.7167 


t 1 ~~i 

1 0.6667 1 

< L. 


0.6000 1 


0.5333 


1 0.4333 1 

i F 


0.3500 


EXACT C.3.F. 


1 .00361 


1 .01343 


1 .02944 1 


.04340 1 


.07275 


1 l 

! .12495 ! 

F 1 


.17929 


12202} NORMAL 


i ‘.00475 
1 


1 ‘.00290 

f 


i ‘.00023 1 


.00336 1 
1 


.00305 


1 .01479 1 

t i 


.01319 


1220 2: NORM. U/CC 


1 ‘.00553 


r 

1 ‘.00413 


} r 

1 ‘.00135 1 

f * 


.00123 1 


.00493 


{ . 
1 .01029 i 


.01226 


E2RQR; T 01 ST 


1 .00235 


1 .00352 


1 .00431 1 


.00459 1 


.00415 


i .00293 I 


.00133 


ERP.0R; T U/CC 


1 .00153 

F 


1 .00203 

i . 


1 .00251 1 

t t 


.00175 1 


.00042 


1 ‘.00212 1 

L _ L 


‘.00479 




P2G3C ?. 1 


rJ ; FOR I 


■AMPLE SIZE 


EQUAL TO 


10. 






TEST STAT. VALUE 


1 0 . 74 3 w 


1 0# o 727 


1 0.5264 i 


0 • ^626 i 


0.4909 1 


1 0.4061 1 

F 


0.3333 


EXACT 0 .1.7. 


1 .00370 


1 .01543 


1 .02722 1 


.04314 1 


.07741 1 


1 .12374 1 


.17437 


22202; NORMAL 


1 '. 00396 


1 ‘.00230 


1 ‘.00091 i 

i r 


.00271 1 


.00700 : 


.01215 i 

F 


.01572 


ERROR; NORM. U/CC 


1 ‘.00437 
1 


1 ‘.00327 
!. 


1 ‘.00210 1 
I - i 


.00095 1 
1 


.00451 1 

... ... . | 


.00367 1 
1 


.01123 


E220R; T DIST 


1 .00204 
» 


1 

1 .00236 

i 


1 1 

1 .00326 1 

i ’ 


f 

.00323 1 

!■ ■ 


t 

.00253 1 

■ - 1 


I 

.00159 1 

i 


.00107 


£2202; T U/CC 


1 .00144 

t 


i .00135 


t 1 

1 .00134 1 

j i 


r 

.00114 1 ' 

i 


t 

'.00035 1 

i 


I 

‘.00230 1 
» 


‘.00361 


the absolute 


value 


of r , 


the observed 


value 


of the 


test 



statistic. The P-values are computed as follows 
a. One-Sided Alternatives 

The one-sided alternative tested 
on the sign of r. A positive r will test for 
association, while negative r tests for 
association. The P-value equals Pr[R >_ |r|]. 



depends 
dir “^t 
indirect 
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b. Two-Sided Alternative 



The P-value equals twice the probability 
computed for the one-sided hypothesis. 

For sample sizes greater than 10, an 
approximation based on the student’s T distribution 
with (N-2) degrees of freedom and continuity, 

correction, is used. The P-values are: 

1 - Pr[ T( N _ 2 ) 1 ( ( | r | - CC) - u r )/ <7 r ] , 

where T(n_ 2) denotes the T distil ibut ion with (N-2) 
degrees of freedom, CC is the continuity correction, 
fj r = 0, and <7 r 2 = (1 - ( | r ! -CC )2 ) / ( n- 2 ) , for the one- 

sided test, and twice this P-value for the two-sided 
test. Gibbons [Ref. 1:pp. 218]. 
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VII.NONPARAMETRIC SIMPLE LINEAR REGRESSION 1 



Nonpar ametr ic Linear Regression assumes that the 
data consists of independent pairs of observations from 
a bivariate disxribution and that the regression of Y 
on X is linear. The program estimaxes linear 
regression parameters based on the data samples. It 
then allows the user to input X values to predict the Y 
values. Hypoxhesis Xesting and confidence interval 
estimation for the slope of the regression equation is 
offered. If the estimated slope lies outside the 
confidence interval, an alternate regression equation 
is offered with an opportunity to input X values to 
predict the corresponding Y values. 

A. COMPUTATION OF THE ESTIMATED REGRESSION EQUATION 

The least squares method is used to estimate A and 
B in .the regression equation Yi = A + BXi + ej_ 
( i= 1 , 2 , . . . N ) , where ei (unobservable errors) are 
assumed to be independent and identically distributed. 
A and B are computed from the following equations: 



1 Except 
information 
paraphrased 



for program design considerations, 
and concepts provided in the section 
from Conover [Ref. 12:pp. 263-271]. 



the 

are 
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N N N 



n E x i Y i - E x i E Y i 



3 



i = 1 i= 1 i = 1 





N N 



E Yi - B £ x i 



A 




N 



B. HYPOTHESIS TESTING 

P-values for testing hypotheses about the slope of 
the regression equation are based on the Spearman’s 
rank correlation coefficient R between the Xi and Uj_ = 



appropriate one-sided test of hypothesis, H 0 : B = 3 $ 

versus HI : B < Bgj or HI : B > B$ , is automatically 
chosen based on the sign of the computed test statistic 
r (positive r tests, HI : B > B @ ; negative r tests, HI : 
B < B0 ) . The P-value is computed as: Pr [R _> |r|]. 

P-values for two-sided tests, H 0 : B = Bq versus HI : B + 
B0 , are also presented. 

For sample sizes N greater than 10 , P-values are 
approximated using a T distribution with (N- 2 ) degrees 

of freedom and continuity correction. 

C. CONFIDENCE INTERVAL ESTIMATION 

100(1- a )% confidence bounds for the slope 
parameter B are determined as follows. The n possible 



3 0 X ^ , where 30 is the hypothesized slope. The 
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slopes, Si_i = ( Yj_-Yj )/( X^-Xh ) , are computed for all 
pairs of data (Xj_,Yj_) and (Xj,Yj) such that i < j and 
Xj[ t Xj and rearranged in increasing order to give s( 1 ) 
<_ s(2) ... <_ s( n ). Let w be the (1- a/2) percentile 
of the distribution of Kendall’s statistic with sample 
size n. 1 Let d be the largest integer less than or 
equal to (n-w)/2 and u the smallest integer greater 
than or equal to (n+w)/2 + 1 . Then s(6.) and s( u ) are 
the desired lower and upper confidence bounds, 
respectively . 

For sample sizes larger than 13, a normal 
approximation with continui.ty correction is used to 
estimate the confidence intervals. 

If the slope of the estimated regression equation 
does not lie within the computed confidence interval, 
the program automatically calculates a new regression 
equation where the slope is the median of the two-point 
slopes Sij and the intercept is the difference of the 
medians of the X and Y samples, My- Mx • 2 

^Kendall’s statistic is defined here as N c - N^, 
where N c is the number of concordant pairs of 
observations and Na is the number of discordant pairs. 
3onover [REF. 1 2 : pp . 256]. 

^This procedure is recommended by Conover [REF. 

1 2 : pp . 256] . 
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VI I I. AREAS FOR FURTHER WORK 



To create a more versatile and powerful software 
package, the NONPAR workspace could be expanded to 
include some or all of the following nonpararaetric 
xests : tests for randomness based on runs, Chisquare 
and Kolmogorov-Smirnov( K-S ) Goodness-of-f it tests, 
Chisquare and K-S general two sample distribution 
xests, Chisquare test for independence , and the 
Friedman test for association. 
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APPENDIX A 



DOCUMENTATION FOR THE MICROCOMPUTER WORKSPACE 



1 . General Information 

This appendix describes the organization and 
operation of the IBM- PC (or compatible) version of the 
workspace. Appendix C continues from where this 
appendix leaves off, to walk the user through each test 
by working practical examples. 

Before proceeding any further, the user should 
refer to section II (Workspace Design Issues) for 
general information about workspace requirements and 
assumptions regarding its use. 

To get started, enter the APL environment in the 
usual manner and load the NONPAR workspace. 

2 . Workspace Menus 

This workspace is designed around the use of 
menus. They guide the user through the selection 
process of choosing a nonpar ametr ic test and a test 
option. Three types of menus are used; the main menu, 
test menus, and help menus. 

a. The Main Menu 

Within moments of loading the NONPAR 
workspace, the main menu will appear. It is titled 
Nonpar ametric Statistical Tests. This menu presents 
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general information about the workspace. Its primary 
purpose is to list the choices of nonparametr ic tests 
available and provide an option which allows the user 
to exit the main menu into APL to copy data into the 
workspace or return to DOS. Each test choice is listed 
with some information about the test’s area of 
application. To make a selection from the menu, move 
the cursor (using the cursor keys) to highlight the 
desired choice. and press enter. As a reminder to the 
user. a footnote at the bottom of the screen describes 
the procedure for entering a choice. Once a test has 
been selected from the main menu, a sub-menu 
appropriate to the test appears. To exit from any menu 
back to the main menu, press the Escape key. 
b. Test Menus 

The title of the test menu is the name of the 
nonpar ametr ic test chosen. The text portion of the 
menu gives a general overview of the test, to include, 
the method used to compute the test statistic, and a 
description of the various options that may be 

exercised. The third section consists of the list of 

test options available. These options include 
returning to the main menu or choosing the help menu. 
Test menus may have options listed in single or 

multiple-paged formats. The comment in the final block 
of the menu lets the user know if a certain menu is 
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mul tiple-paged or not. To make a selection from a 
multiple-paged qjenu, use the page-up or page-down key 
to locate the desired option. Proceed with the scroll 
keys to highlight the choice. and press enter. Once a 
test option is entered, the user is prompted to input 
the data required to run the test. When the option for 
more information is selected, the help menu is 
displayed. 

c. Help Menus 

The title of the help menu usually begins 
with the words "More Information About..." followed by 
the title of the nonpar ametr ic test. The text portion 
of the menu explains the test and its options in 
greater detail. No choices are offered in the menu. 
To return to the test menu, press any key. 
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APPENDIX 3 



DOCUMENTATION FOR THE MAINFRAME COMPUTER WORKSPACE 
1 . General Information 

This appendix describes the organization and 
operation of the mainframe computer workspace. To load 
a copy of the NONPAR workspace from the APL library, 
enter the APL environment and type: )LOAD 9 NONPAR. 

Within a few moments the variables LIST and DESCRIBE 
are displayed on the screen. These variables provide a 
description of the workspace. 

2 . The NONPAR Workspace 

The NONPAR workspace consists of seven 

programs which call several subprograms during their 
execution. The exact syntax for each test and its 
corresponding nonpar ametric test name is given in the 
following format: 

SYNTAX: Nonpar ametr ic Test and Application. 

a. SIGN: Ordinary Sign Test for Location in Single 

and Paired-sample Data. 

b. WILCOX: Wilcoxon Signed-rank Test for Location 

in Single and Paired-sample Data. 

c. MANNWHIT: Mann-Whitney Test for Equal Medians or 

Variances in Two Independent Samples. 

d. KRUSKAL: Kruskal-Wallis Test for Equal Medians in 

K Independent Samples. 
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e. KENDALL: Kendall’s B; Measure of Association for 

Paired-sample Data. 

f. SPEARMAN: Spearman’s R; Measure of Association 

Between Rankings of Paired Data. 

g. NPSLR: Nonpar ametr ic Simple Linear Regression; 

Least Squares . 

The list presenxed above can be displayed at any 
time by typing: LIST. 

For each test program, there exists a HOW variable 
that gives a full description of the test and the 
various options that may be exercised. To display any 
of the HOW variables, just enter the test program’s 
name with the suffix HOW appended (i.e. SIGNHOW). 

A test is run by entering the program’s name. The 
user is immediately prompted to input data. Enter 
numerical data separated by spaces or as a variable to 
which the numbers have been previously assigned. 
Several of the tests require a considerable amount of 
prompting before all the necessary data has been 
entered . 
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APPENDIX C 



WORKSPACE FAMILIARIZATION THROUGH PRACTICAL EXAMPLES 
1 . General Information 

This appendix applies to both the mainframe and 
microcomputer workspaces. Its purpose is to acquaint 
the user with the organization of the programs and the 
type of prompts to be expected. 

Extensive error checking has been included in the 
programs to ensure that the data is of the proper form. 
Should a program become suspended, clear the state 
indicator by entering: )RESET , check over the data for 
errors, and restart the program. 330 kilobytes of 
computer memory are needed to load APL and the NONPAR 
workspace; to avoid filling up the remaining workspace 
area, the user should minimize data storage in the 
NONPAR workspace. To exit a program at any time, press 
the Control and Escape keys, simultaneously. 

2 . Practical Examples 
a. Sign Test 

(1 ) Description of Problem 1 . A Sinclair mine 
is manufactured to have a median explosive weight of 
not less than 16 ounces. The explosive weights of 15 
mines, randomly selected from the production line, were 
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recorded as follows: 



16.2 15.7 15.9 15.8 15.9 16 16.1 



15.3 15.9 16 16.1 15.7 15.8 15.9 15.8. 

(a) Is the manufacturing process packing 
enough explosives in the mines? 

(b) What range of values can be expected 
for the median of the explosive weights 90% of the 
t ime . 

(2) Solution . To see if the manufacturing 
process is meeting the specifications, we test the 
hypothesis Ho: M = 16 versus HI : M < 16. 

( 5 ) Workspace Decision Process . 

(a) Microcomputer: Choose the Sign Test 
from the main menu, and the option, Single Sample; 
Test H0 : M = Mo versus HI : M < Mo, from the test menu. 
Skip to the Program Interaction section below. 

(b) Mainframe: Enter SIGN at the keyboard 
and receive the prompt: 

DID YOU ENTER THIS PROGRAM FOR THE SOLE 
PURPOSE OF GENERATING CONFIDENCE INTERVALS FOR A 
SPECIFIED SAMPLE SIZE AND QUANTILE? (Y/N). 

Enter N (If Y is entered, the user will 
go directly to this last option of the test). The next 
prompt is: 

THE NULL HYPOTHESIS STATES - THE 
POPULATION MEDIAN (M) IS EQUAL TO THE HYPOTHESIZED 
MEDIAN (Mo); H0 : M = Mo. WHICH ALTERNATIVE DO YOU WISH 
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TO TEST? ENTER: 1 FOR HI: M < Mo; 2 FOR HI: M > Mo; 

3 FOR HI : M t Mo . 

Enter 1 . The next prompt is: 

ENTER: 1 FOR SINGLE-SAMPLE PROBLEM; 2 

FOR PAIRED-SAMPLE PROBLEM. 

Enter 1 . 

(4) Program Interaction . The prompt is: 

ENTER THE DATA (MORE THAN TWO OBSERVATIONS 
ARE REQUIRED). 

Enter the data separated by spaces or as a 
variable to which the data has been previously 
assigned. The next prompt is: 

ENTER THE HYPOTHESIZED MEDIAN. 

Enter 16. The following is dislayed. 
COMPUTATIONS ARE BASED ON A SAMPLE SIZE OF: 

1 3 . 

THE TOTAL NUMBER OF POSTIVE SIGNS IS: 3. 

THE P-VALUE FOR H0 : M = 16 Versus HI: 

M < 1 6 IS: . 0461 . 

Consider a significance level of .05. 
Since the P-value of .0461 is less than .05, we reject 
H0 : M = 16 in favor of HI : M < 16 and conclude that the 

manufacturing process is not packing enough explosives 
in the Sinclair mine. The next prompt is: 

WOULD YOU LIKE A CONFIDENCE INTERVAL FOR 
THE MEDIAN? (Y/N). 
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Enter Y ( If N is entered, the progam asks 
if confidence intervals for a quantile are desired). 
The next prompt is: 

ENTER THE DESIRED CONFIDENCE COEFFICIENT; 
FOR EXAMPLE: ENTER 95, FOR A 95% CONFIDENCE INTERVAL. 

Enter 90. The following is displayed. 

A 90% CONFIDENCE INTERVAL FOR THE MEDIAN OF 
THE POPULATION IS: ( 1 5 . 8 <_ MEDIAN £ 16 ). 

The next prompt is: 

WOULD YOU LIKE CONFIDENCE INTERVALS FOR A 
SPECIFIED QUANTILE? (Y/N). 

To see the form of the results, we generate 
confidence intervals for the 30th quantile. Sample 
size is automaticly set at the number of data points 
entered eariler . Enter Y (If N is entered, the 
mainframe program ends; or, the Sign test menu 
reappears). The next prompt is: 

ENTER DESIRED QUANTILE; FOR EXAMPLE: ENTER 
20, FOR THE 20TH QUANTILE. 

Enter 30. The following is displayed. 

ORDER STATISTICS _!_ COEFFICIENTS 
3 8 ! .823160 

2 9 ! .949490 

1 ‘0 , .991600 

***** THIS TABLE GIVES CONFIDENCE COEF- 
FICIENTS FOR VARIOUS INTERVALS WITH ORDER STATISTICS AS 
END POINTS FOR THE 30TH QUANTILE. 
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The menu- 



The mainframe program ends, 
driven microcomputer program pauses for input from the 
keyboard by prompting: 

PRESS ENTER WHEN READY. 

Press Enter and the Sign test menu 

reappears . 

b. Wilcoxon Signed-rank Test 

( 1 ) Description of Problem 2 . A special train- 
ing program is being considered to replace the regular 
training that Radio Telephone Operators receive. In 
order to evaluate the effectiveness of the new training 
program, proficiency tests were given during the third 
week of regular training. Twenty-four trainees were 
chosen at random and grouped into twelve pairs based on 
proficiency test scores. One member of each pair 
received specialized training while the other member 
received regular training. Upon graduation, the 
proficiency tests were given again with the following 
results . 

Specially Trained Group (X): 60 50 55 71 

45 59 64 49 61 54 47 70 

Regularly Trained Group (Y): 4.0 46 60 55 

49 57 51 55 4-5 59 40 55 

(a) Does the special training program 
ensure higher scores? 
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(b) By what range of values can the 
scores of the two groups be expected to differ 95 % of 
the time? 

(2) Solution . To test the hypothesis that the 
special training program raises profficiency scores, we 
test H0 : M(X-Y) = 0 versus HI: M(X-Y) > 0. 

( 3 ) Workspace Decision Process 

(a) Microcomputer: Choose the Wilcoxon 

Signed-rank Test from the main menu, and the option. 
Paired-sample; Test H0 : M = Mo versus HI : M > Mo, from 
the test menu. Skip to the Program Interaction section 
below. 

(b) Mainframe: Enter WILCOX at the key- 

board and receive the prompts: 

THE NULL HYPOTHESIS STATES - THE 
POPULATION MEDIAN (M) IS EQUAL TO THE HYPOTHESIZED 
MEDIAN (Mo); H0 : M = Mo. WHICH ALTERNATIVE DO YOU WISH 
TO TEST? ENTER: 1 FOR HI: M < Mo; 2 FOR HI: M > Mo; 
3 FOR HI : M * Mo . 

Enter 2. The next prompt is: 

ENTER: 1 FOR SINGLE-SAMPLE PROBLEM; 2 

FOR PAIRED-SAMPLE PROBLEM. 

Enter 2. 

(4) Program Interaction . The prompt is: 

ENTER X DATA (MORE THAN TWO OBSERVATIONS 
ARE REQUIRED). 
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Enter the X data separated by spaces. The 
next prompt is: 

ENTER Y DATA (NUMBER OF Y ENTRIES MUST 
EQUAL NUMBER OF X ENTRIES). 

Enter the Y data. The next prompt is: 

ENTER THE HYPOTHESIZED MEDIAN FOR THE 
DIFFERENCES OF THE PAIRED DATA. 

Enter 0. The following is dislayed. 

COMPUTATIONS ARE BASED ON A SAMPLE SIZE OF: 

1 2 . 

THE TOTAL SUM OF POSITIVE RANKS IS: 60.5. 

THE P-VALUE FOR COMPARING THE MEDIAN OF THE 
POPULATION OF DIFFERENCES TO THE HYPOTHESIZED MEDIAN, 
H0 : M(X-Y) = 0 Versus HI: M(X-Y) > 0, IS: .0505. 

Consider a significance level of .05. 
Since the P-value of .0505 is greater than .05, we do 
not reject the null hypothesis that the two training 
cources are equally effective. However, due to the 
closeness in values, the choice of rejecting or not 
rejecting the null hypothesis is strictly a judgement 
call. The next prompt is: 

WOULD YOU LIKE A CONFIDENCE INTERVAL FOR 
THE MEDIAN? (Y/N). 

Enter Y ( If N is entered, the mainframe 
progam ends; or, the Wilcoxon test menu reappears). 
The next prompt is: 
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ENTER THE DESIRED CONFIDENCE COEFFICIENT; 
FOR EXAMPLE: ENTER 95, FOR A 95^ CONFIDENCE INTERVAL. 

Enter 95. The following is displayed. 

A 95 % CONFIDENCE INTERVAL FOR THE MEDIAN OF 
THE POPULATION OF DIFFERENCES IS: 

( -1 <_ MEDIAN ( X-Y ) ± 16.5 ). 

The mainframe program ends. The menu- 
driven microcomputer program pauses for input from the 
keyboard by prompting: 

PRESS ENTER WHEN READY. 

Press Enter and the Wilcoxon test menu 

reappears . 

c. Mann-Whitney Test for Equality of Medians 

( 1 ) Description of Problem 5 . A group of Army 
and Navy officers were given the Defense Language 
Aptitude test. From the results, 14 Army and 17 Navy 
officers’ scores were randomly selected. These scores 
are listed below. 

Army (X): 35 30 55 51 28 25 16 63 60 44 20 

42 47 38. 

Navy (Y): 54 26 41 43 37 34 39 50 46 49 45 

33 29 36 38 42 34. 

(a) Is there sufficient evidence to claim 
that Navy officers score higher on this test than Army 
officers? 
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(b) By what range of values can the 
scores between the two groups be expected to differ 90% 
of the time. 





(2) 


Solution . 


To see 


if Navy 


officers 


score 


higher 


on 


the exam 


, we test 


H0 : Mx = My 


versus 


HI : Mx < 


My. 














(3) 


Workspace 


Decision Process. 










(a) Microcomputer: 


Choose 


the 


Mann- 


Whitney 


Test 


from the 


main menu. 


and the 


option , 


Test 



H0 : Mx = My versus HI : Mx < My, from the test menu. 

Skip to the Program Interaction section below. 

(b) Mainframe: Enter MANNWHIT at the 

keyboard and receive the prompts: 

DO YOU WISH TO COMPARE THE MEDIANS OR 
VARIANCES OF THE POPULATIONS? ENTER: 1 TO COMPARE 

MEDIANS; 2 TO COMPARE VARIANCES. 

Enter 1. The next prompt is: 

THE NULL HYPOTHESIS STATES - THE 
MEDIANS OF X AND Y ARE EQUAL; Mx = My. WHICH 
ALTERNATIVE DO YOU WISH TO TEST? ENTER: 

1 FOR HI : Mx < My; 2 FOR Mx > My ; 3 FOR Mx t My. 

Enter 1 . 

(4-) Program Interaction . The prompt is: 

ENTER X DATA (MORE THAN ONE OBSERVATION IS 

REQUIRED) . 
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Enter the X data separated by spaces. The 
next prompt is: 

ENTER Y DATA. 

Enter the Y data. The following is 

displayed. 

THE SUM OF THE X RANKS IS: 224. THE U 

STATISTIC EQUALS: 119. 

THE P- VALUE FOR H0 : Mx =My versus HI : 

Mx < My IS: .5078. 

We do not reject the hypothesis of equal 
population medians and conclude the median of all Army 
scores is equal to the Navy’s. The next prompt is: 

WOULD YOU LIKE A CONFIDENCE INTERVAL FOR 
THE SHIFT IN LOCATION (My - Mx)? (Y/N). 

Enter Y (If N is entered, the mainframe 
program ends; or, the Mann-Whitney test menu 
reappears). The next prompt is: 

ENTER THE DESIRED CONFIDENCE COEFFICIENT; 
FOR EXAMPLE: ENTER 95, FOR A 95% CONFIDENCE INTERVAL. 

Enter 95. The following is displayed. 

A 95% CONFIDENCE INTERVAL FOR THE SHIFT IN 
LOCATION BETWEEN POPULATIONS X AND Y IS: 

(-10 _£ My-Mx <_ 1 0 ) . 

The mainframe program ends. The microcom- 
puter program pauses for input from the keyboard by 
prompting : 
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PRESS ENTER WHEN READY. 



Press Enter and the Mann-Whitney test menu 



reappears . 

d. Mann-Whitney Test for Equality of Variances 

(1) Description of Problem .- Referring to 
problem 3 in section c(1). Is there sufficient 
evidence to claim that Army scores vary more than Navy 
scores? 



(2) Solution . To see if Army 
more, we test H0 : Vx = Vy versus HI : Vx > 

( 3 ) Workspace Decision Process , 

(a) Microcomputer: Choose 

Whitney Test from the main menu, and the 
H0 : Vx = Vy versus HI : Vx > Vy , from the 
receive the prompt: 



scores vary 
Vy. 

the Mann- 
option. Test 
test menu and 



ENTER THE DIFFERENCE OF THE MEANS OR 
MEDIANS (Mx - My). 



Because we believe the population 
medians to be approximately equal. We enter 0. Skip to 
the Program Interaction section below. 

(b) Mainframe: Enter MANNWHIT at the 
keyboard and receive the prompts: 

DO YOU WISH TO COMPARE THE MEDIANS OR 
VARIANCES OF THE POPULATIONS? ENTER: 1 TO COMPARE 
MEDIANS; 2 TO COMPARE VARIANCES. 

Enter 2. The next prompt is: 
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THE TEST TO COMPARE VARIANCES, REQUIRES 
THE TWO POPULATION MEANS OR MEDIANS TO BE EQUAL. IF 
THEY DIFFER BY A KNOWN AMOUNT, THE DATA CAN BE ADJUSTED 
BEFORE APPLYING THE TEST. ENTER THE DIFFERENCE OF 
MEDIANS (Mx -My) OR 900 TO QUIT. 

We enter 0. The next prompt is: 

THE NULL HYPOTHESIS STATES - THE 
VARIANCES OF X AND Y ARE EQUAL; Vx = Vy . WHICH 
ALTERNATIVE DO YOU WISH TO TEST? ENTER: 

1 FOR HI : Vx < Vy; 2 FOR Vx > Vy ; 3 FOR Vx t Vy . 

Enter 2. 

(4) Program Interaction . The prompt is: 

ENTER X DATA (MORE THAN ONE OBSERVATION IS 

REQUIRED) . 

Enter the X data separated by spaces. The 
next prompt is: 

ENTER Y DATA. 

Enter the Y data. The following is 

displayed. 

THE SUM OF THE X RANKS IS: 166. THE U 

STATISTIC EQUALS: 61. 

THE P-VALUE FOR H0 : Vx = Vy versus HI: 

Vx > Vy IS: .0112. 

Consider a significance level of .05. 
Since a P-value of .0112 is less than .05, we reject 
the null hypothesis of equal variances in favor of 
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Vx > Vy and conclude that Army scores do vary more than 
Navy scores. 

The mainframe program ends. The 

microcomputer program pauses for input from the key- 
board by prompting: 

PRESS ENTER WHEN READY. 

Press Enter and the Mann-Whitney test menu 

reappears . 

e. Kruskal-Wallis Test 

(1) Description of Problem 4 . During a recent 
Monster Mash involving four Navy SEAL Teams, one of the 
events consisted of the number of pushups a man could 
do in 2 minutes. Eight men were chosen randomly from 
each Team. The following scores were recorded. 

SEAL 1: 90 96 102 85 65 77 88 70. 

SEAL 2: 64 79 99 95 87 74 69 97. 

SEAL 5: 101 66 93 89 71 60 76 98. 

SEAL 4: 72 78 73 81 83 92 94 86. 

Are the different Seal Teams considered to 
be equally fit? 

(2) Solution . To see if the Seal Teams are 

equally fit, we test the hypothesis that all the 
population medians are equal . 
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( 3 ) Workspace Decision Process . 



(a) Microcomputer: Choose the Kruskal- 

Wallis Test from the main menu; and, once the test menu 
is displayed, press Enter. 

(b) Mainframe: Enter KRUSKAL at the 

keyboard. 

(4) Program Interaction . The prompt is: 

ENTER THE NUMBER OF POPULATIONS TO BE 
COMPARED (MUST BE GREATER THAN TWO). 

Enter 4. The next prompt is: 

ENTER YOUR FIRST SAMPLE. 

Enter the SEAL 1 data separated by spaces. 
The next prompt is: 

ENTER YOUR NEXT SAMPLE. 

Enter the SEAL 2 data. The next prompt is: 
ENTER YOUR NEXT SAMPLE. 

Enter the SEAL 3 data. The next prompt is: 
ENTER YOUR LAST SAMPLE. 

Enter the SEAL 4 data. The following is 

displayed. 

THE H STATISTIC EQUALS: .1335. 

THE P-VALUE FOR H0 : THE POPULATION MEDIANS 

ARE EQUAL versus HI : AT LEAST TWO POPULATION MEDIANS 

ARE NOT EQUAL IS: .98893. 



67 



We do not reject the null hypothesis that 
the population medians are equal and conclude that the 
SEAL Teams are equally fit. 

The mainframe program ends. The microcom- 
puter program pauses for input from the keyboard by 
prompting : 

PRESS ENTER WHEN READY. 

Press Enter and the Kruskal-Wallis test 
menu reappears. 

f. Kendall’s B 

(1 ) Description of Problem 5 . In order to 
determine if cold weather affects target marksmanship, 
Naval Special Warfare recorded small arms marksmanship 
scores and corresponding air temperatures for a period 
of one year. 20 men were chosen at random, and their 
scores averaged for different air temperatures. The 
average score for each air temperature is shown below. 

Air temperature (X): 50 55 20 50 65 55 30 

52 40 60. 

Average scores (Y): 210 200 165 165 260 

215 175 191 180 235. 

Can it be said that colder temperatures 
have an effect on marksmanship scores? Is that effect 
positive or negative? 
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(2) Solution . We test the null hypothesis 
that no association exists between cold temperatures 
and marksmanship . 

( 3 ) Workspace Decision Process . 

(a) Microcomputer: Choose Kendall’s B 

Test from the main menu; and, once the test menu is 
displayed, press Enter. 

(b) Mainframe: Enter KENDALL at the 

keyboard . 

(4) Prograjn Interaction . The prompt is: 

ENTER X DATA (MORE THAN TWO OBSERVATIONS 
ARE REQUIRED). 

Enter the X data separated by spaces. The 
next prompt is: 

ENTER Y DATA (NUMBER OF Y ENTRIES MUST 
EQUAL NUMBER OF X ENTRIES). 

Enter the Y data. The following is 

displayed. 

KENDALL’S B EQUALS: .7817. 

THE P-VALUE FOR H0 : NO ASSOCIATION EXISTS 

versus: HI: DIRECT ASSOCIATION EXISTS IS: .00045. 

THE P-VALUE FOR THE TWO-SIDED TEST OF 
HYPOTHESIS IS: .0009. 

Since the P-value for the one-sided test 
equals .00045, we reject the null hypothesis that no 
association exists between temperatures and 
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marksmanship in favor of direct association. We 
conclude that colder temperatures tend to cause lower 
marksmanship scores. 

The mainframe program ends. The microcom- 
puter program pauses for input from the keyboard by 
prompting: 

PRESS ENTER WHEN READY. 

Press Enter and Kendall’s B test menu 

reappears . 

g. Spearman’s R 

(1 ) Description of Problem 6 . When fitness 
reports are written, officers of the same grade are 
ranked against each other based upon their demonstrated 
level of performance. Last marking period, the 
Commanding and Executive Officers separately ranked 9 
Ensigns as shown below. 

Ensigns 

ABCDEFGHI 
CO (X): 64 1528379 

XO ( Y ) : 5634 1 9728 

Does any association exist between the two 
sets of rankings? 

(2) Solution . We test the null hypothesis 
that no association exists. 
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( 3 ) Workspace Decision Process . 



(a) Microcomputer: Choose Spearman’s R 

Test from the main menu; and, once the test menu is 
displayed, press Enter. 

(b) Mainframe: Enter SPEARMAN at the 

keyboard. 

(4) Program Interaction . The prompt is: 

ENTER X DATA (MORE THAN TWO OBSERVATIONS 
ARE REQUIRED). 

Enter the X data separated by spaces. The 
next prompt is: 

ENTER Y DATA (NUMBER OF Y ENTRIES MUST 
EQUAL NUMBER OF X ENTRIES). 

Enter the Y data. The following is 

displayed. 

SPEARMAN’S R EQUALS: .5500. 

THE P-VALUE FOR H0 : NO ASSOCIATION EXISTS 

versus: HI: DIRECT ASSOCIATION EXISTS IS: .0664. 

THE P-VALUE FOR THE TWO-SIDED TEST OF 
HYPOTHESIS IS: • .1328. 

Consider a significance level of .05. 
Since a P-value of .0664 exceeds .05, we do not reject 
the null hypothesis that no correspondence exists 
between the two sets of rankings. 
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The mainframe program ends. 



The microcom- 



puter program pauses for input from the keyboard by 
prompting : 



reappears . 

h. Nonparametr ic Simple Linear Regression; Least 
Squares 

( 1 ) Description of Problem 7 . Battery-powered 
Swimmer Proplusion Units are sometimes used to aide 
swimmers during long underwater swims. Recent tests 
have shown that a nearly linear relationship exists 
between water temperature and battery life for these 
units. The following 17 data points were randomly 
selected from the test results. 



PRESS ENTER WHEN READY. 



Press Enter and the Spearman’s R test menu 



Water temperature (X) 



Battery life (Y) 



70 

65 

50 

40 

60 

55 

52 

50 

43 

40 

72 

55 

48 

35 

70 

68 

57 



3 



2.75 
1 .8 
1 .2 
2.4 
1 .9 
1 .75 

1 .7 
1 .6 
1 . 1 
2.75 

2 



1 -5 
.9 

3.3 

3 



2.3 



(a) Find the fitted regression equation. 
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(b) For the following water temperatures, 

predict the battery life of the units: 61 52 46 56. . 

(c) Can we determine with any certainty 
if the slope of the regression line equals -.05. 

(d) What range of values could be used 
as the slope of the estimated equation line 9<Z>% of the 
time? 

(2) Solution ♦ To determine the estimated 
regression equation, we use nonparametr ic linear 
regression . 

( 3 ) Workspace decision process . 

(a) Microcomputer: Choose Nonpar ametr ic 

Simple Linear Regression from the main menu; and, once 
the test menu is displayed, press Enter. 

(b) Mainframe: Enter NPSLR at the 

keyboard . 

(4) Program Interaction . The prompt is: 

ENTER X DATA (MORE THAN TWO OBSERVATIONS 

ARE REQUIRED). 

Enter the X data separated by spaces. The 
next prompt is: 

ENTER Y DATA (NUMBER OF Y ENTRIES MUST 
EQUAL NUMBER OF X ENTRIES). 

Enter the Y data. The following is 

displayed . 
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THE LEAST SQUARES ESTIMATED 



REGRESSION 



EQUATION IS: 

Y = -1 .263 + . 060668X . 

The next prompt is: 

DO YOU WISH TO ENTER SOME X VALUES TO GET 
THE PREDICTED Y ’ S? (Y/N). 

Enter Y ( If N is entered, the program skips 
to hypothesis testing for the slope). 

ENTER X VALUES. 

Enter 61 52 46 36. The next prompt is: 

THE PREDICTED Y VALUES ARE: 2.44 1.89 

1.53 .92. 

WOULD YOU LIKE TO RUN SOME MORE X VALUES? 
Enter N. The next prompt is: 

WOULD YOU LIKE TO TEST HYPOTHESIS ON B, THE 
SLOPE OF THE EQUATION? (Y/N). 

Enter Y ( If N is entered, the program skips 
to confidence interval estimation). The next prompt 
is : 

ENTER THE HYPOTHESIZED SLOPE. 

Enter .05. The following is displayed. 
SPEARMAN’S R EQUALS: .5756. 

THE P-VALUE FOR H0 : B = .05 versus HI : B > 

.05 IS: .0079. 

THE P-VALUE FOR THE TWO-SIDED TEST OF 
HYPOTHESIS IS: .0158. 
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Consider a significance level of .05. Since 
a P-value of. 0079 is less than .05, we reject the null 
hypothesis that B = .05 in favor of B > .05, and 
conclude the slope of the regression line is greater 
than .05. The next prompt is: 

WOULD YOU LIKE A CONFIDENCE INTERVAL FOR 
THE SLOPE? (Y/N). 

Enter Y ( If N is entered, the mainframe 
program ends; or, the Nonpar ametr ic Regression test 
menu reappears). The next prompt is: 

ENTER THE DESIRED CONFIDENCE COEFFICIENT; 
FOR EXAMPLE: ENTER 95, FOR A 95£ CONFIDENCE INTERVAL. 

Enter 90. The following is displayed. 

A 90^ CONFIDENCE INTERVAL FOR B, THE SLOPE 
OF THE ESTIMATED REGRESSION LINE, IS: 

( .05333 <_ B <_ .07 ). 

If the estimated slope, does not lie within 
the confidence interval, the following would he 
displayed . 

THE LEAST SQUARES ESTIMATOR OF B LIES 
OUTSIDE THE CONFIDENCE INTERVAL. DISCARD THE LEAST 
SQUARES EQUATION AND USE: 

Y = -1.4-458 + . 060833X . 

THIS EQUATION IS BASED ON THE MEDIANS OF 
THE X AND Y DATA, AND THE MEDIAN OF THE TWO- POINT 
SLOPES CALCULATED FOR THE CONFIDENCE INTERVAL ON B. 
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The next prompt is: 

DO YOU WISH TO ENTER SOME X VALUES TO GET 
THE PREDICTED Y’S FROM THE NEW EQUATION? (Y/N). 



To 


compare results, 1 


et us 


input 


the 


temperatures 


in the new equation. 


Enter Y ( If N 


is 


entered, the 


mainframe program 


ends ; 


or , 


the 


Nonpar ametric 


Regression test menu 


reappears ) . 


The 



next prompt is: 

ENTER X VALUES. 

Enter 61 52 46 36. The following is 

displayed. 

THE PREDICTED Y VALUES ARE: 2.265 1.72 

1.35 .74. 

The next prompt is: 

WOULD YOU LIKE TO RUN SOME MORE X VALUES? 
Enter N. The next prompt is: 

WOULD YOU LIKE TO TEST HYPOTHESIS ON B, THE 
SLOPE OF THE EQUATION? (Y/N). 

To compare results once again, we enter Y 
(If N is entered, the mainframe program ends; or, the 
Nonparametr ic Regression test menu reappears). The 
next prompt is: 

ENTER THE HYPOTHESIZED SLOPE. 

Enter .05. The following is displayed. 
SPEARMAN’S R EQUALS: .8287. 
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THE P-VALUE FOR H0 : B 



. 05 versus HI : B > 



.05 IS: .0000. 



THE 


P-VALUE FOR THE TWO-SIDED TEST OF 


HYPOTHESIS IS: 


. 0000 . 


The 


mainframe program ends. The microcom- 


puter program 


pauses for input from the keyboard by 


prompting : 





PRESS ENTER WHEN READY. 

Press Enter and the Nonparametr ic 
regression test menu reappears. 
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APPENDIX D 



MAIN PROGRAM LISTINGS FOR MICROCOMPUTER WORKSPACE 



p 

L9 j 

C 10 * 

111 , 

112, 

[13 

14 

r Id 
;i7 J 
Cia, 

t Q 1 

■r 

ci» 

I? 

:li 

30 
1 3 1 

32 

33 

34 

35 

36 
,37 

:ii 

:52 

:5I 

,44 

:5i 

:48 

>9 
■ $ 

:s4 

.55 

56 

:§? 

[58 

[59 

[60 

,61, 

[62j 

1 6 3 j 
£643 

Ess; 

66 

67 

.68 

.69, 

.70 

.71 

.72 

l73. 



V KENzA :AA ;B;BX:BY; CzCX ;CY:D ;DD; DX ;DY zDXY ; S ; POS ; NEC ; XX ;YY;N; DEN ; NN ; NUM ; P 
PVAL • SU • SV • T • U • V • AT • Z • X • * ' WW • CHA •E*PV*Q*R 
p THIS FtiNCtldN' COMPOTES TEE kENDJuL B* STATISTIC WE ICE IS A MEASURE 
p OF ASSOCIATION BETWEEN SAMPLES. P- VALUES ARE GIVEN FOR TESTING ONE 
A AND TWO-SIDED EIPOTEESIS FOR NO ASSOCIATION VERSUS ASSOCIATION . 
p SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE : TIES , TIESK , KEN D ALP , 
p INTERP , INPUT AND NORMCDF . 

a DISPLAY TEST MENU AND INPUT DATA 

NlzE+MENU KEN03J 
->(E=1)/B1 
MENU MAINOBJ 

Bl:R+INPUT 2 
0-<-l f R 

3 + 1 ) *R 

ORDER Y IN INCREASING ORDER OF X 



A+IZ&X1 

*3+XZ&Xl 

C+&&A 



ORDER X IN INCREASING ORDER 
COMPUTE CURRENT RANKING OF I 
NOW ORDER I RANKS IN INCREASING ORDER 



D+AZU1 

p IF TIES EXIST IN EITHER X OR I RANKED VECTOR USE MID -RANK METHOD 
DD + l TIES D 
XX+1 TIES 3 

p FIND ORIGINAL RANKING OF I WITH TIES RESOLVED 

YY+DDZCl 
N+oX 

p COMPUTE NUMBER OF DISTINGUISHABLE PAIRS V 

NN+(N*(N- l))+2 
S«-p 0 
AA+ 0 

P POSITIVE ONES COME FROM A RUNS UP CONDITION z NEGATIVE 1 FROM RUNS DOWN 
p ZERO IS SCORED FOR TIES. MULTIPLY THE RESULTS FOR EVERY ELEMENT AND SUM 
LllAA+AA+l 
BX+CXXLAA}>(AA*XX)) 

CX+IXXLAA1<UA*XX)UC1) 

DX+BX+CX 

BY+ ( YY ZAA I > CAA +YY ) ) _ 

c y+ IyyiaaUIaa+yy)U(i) 

DY+BY+CY 

DXY+DXxDY 

posAdxy> 0 ) N 

NEgADXY< O)x(-I) 

S+S t POS,NEG 

fl SUM FINAL VECTOR TO DETERMINE S 

S++/S 

p OBTAIN THE NUMBER OF TIES IN EACH VECTOR USING THE TIESK FUNCTION 
U+TIESK B 
V+TIESK D 
SU++/C21U) 

SV++/ (217) 

p CALCULATE THE B STATISTIC INCLUDING THE CORRECTION FOR TIES 
T+S* ( ( NN-SU )x ( NN-SV ) )*0 . 5 
AT+ T 

+ (N> 13 ) /NORM 

p CALL KENDALP TO CALCULATE THE RIGHT TAIL OF THE CDF OF B 

P+KENDALP N 

p CALL INTERP TO CALCULATE P-VALUE BY INTERPOLATION 
PVAL+AT INTERP ? 

+-(PVAL*~1)/L 3 
PVAL+Q . 5 
+L 3 

fl \ CALCULATE P VALUE USING NORMAL APPROX . 

NORM:NUM+(3*AT)x ( (2*NN)*Q . 5 ) 

DEN+(2* ( C2xiV) + 5 ) )*0 , 5 
Z+NUM+DEN 

PVAL+1- (NORMCDF Z) 

P ^POSITIVE PRINT OUT DIRECT ASSOCIATION . 

CHA+ 1 INDIRECT 1 
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£5 :C3 A* 'DIRECT' 

L7 :PV*2*PVAL 
*(PVS1)/L3 
PV* 1 

LB :> KENDALL' ' S S EQUALS: ' .(W) .OTCNL 
'THE P -VALUE FOR HO: NO ASSOCIATION EXISTS VERSUS' 

< HI: ',(*CSi) ' ASSOCIATION EXISTS IS: ' 

•TEE P -VALUE FOR THE TWO-SIDED TEST OF HYPOTHESIS IS: ' 

< PRESS ENTER WHEN READY . ' 

WW+Q 

-tfl 

7 



• , (nvPVAL) .OTCNL 
, (u sPV),OTCNL 



7 KRWLiNUM:DENOM:A:C:H:B:K:AA:BB:DD:E:F:N:OF:P:PVAL:R:SOFR:SR:TSOR:CBA:B 

a this function compOtZs the kIjskal-wallis test statistic h which is 

, 2 ' a A MEASURE OF THE EQUALITY OF K INDEPENDENT SAMPLES. 

3 ’ a SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: TIES, TIESK . INDEXPLS, 

* 4 j a FDISTN , INTER P AND THE VARIABLES ? MAT KW 20 , PMAIKW 31, PMATKW33 , PMATKW 34 
a PMATKW 41, PMATKW 42, AND RMATKW 43. 

|7] a MENU CHOICES AND ROUTE TO PROPER STATEMET FOR ACTION. 

8 ] N 1 : B+MENU KRWL03J 

9] *(B=1)/31 

10] MENU MAI NOB J 
' - * -*-0 — 

■ r -5 ] 2P°-*-4 

III] 31 :* ENTER TEE NUMBER OF POPULATIONS TO BE COMPARED (.MUST 3E GREATER THAN T 
U0). ' 

;i5] i((H<3)v(fl|iO=0))/31 

16 -((p£)>l)/31 

17] a INITIALIZE VECTORS E AND F AND VARIABLE C 

18, E~F~SOFR+ pO 
, 19 ] 00 

20] a THIS LOOP FACILITATES ENTERING THE SAMPLE VECTORS AND STORING THEM 
21, OH A- ' FIRST ' 

22 : £l:C<-Ol 

' ENTER YOUR • , (.9CHA ) , ' SAMPLE. < 

+$poD)*0)/NEXT 
D* lpD 

a CONCATENATE SAMPLES AS THEY ARE ENTERED AND STORE THEM IN VECTOR E 
NEXT * E+E D 

ft * ’RECORD THE LENGTHS OF THE SAMPLES AS THE X ARE ENTERED 
F+F,pD 
CHA+'NEXT' 

+ (C<U(- 1 ) )/Ll 
CHA+'LAST' 

+{C<K)/L 1 

p RECORD SIZE OF ALL SAMPLES WHEN COMBINED 

N++/F 

p ORDER SAMPLE SIZES LARGEST TO SMALLEST 
0F+F11F^ 

p 4 ORDER COMBINED SAMPLE VECTOR TO BE USED BY TIES FUNCTION 
D+ElkEl 

p CALL INDEXPLS TO INCREMENT INDEXES WHEN TIES OCCUR WITHIN ONE SAMPLE 
AA+F INDEXPLS E 

p CALL TIES TO BREAK TIES BY MIDRANK METHOD 

B3+1 TIES D 
C + 0 

p THIS LOOP CALCULATES THE H STATISTIC 

L2 : C+C+l 

« EACH SAMPLE IS CALCULATED 

SR++/BBZ (FlCl ♦ (AA [C: ] ) )] 

* CALCULATE SUM OF RANKS SQUARED DIVIDED BY THE INDIVIDUAL SAMPLE SIZE 
SR+ v SR * 2 j LC J 

a„„„„ STORE EACH CALCULATION 

SOFR+SOFR , SR 
*(C<K)/L 2 

« „ , SUM ACROSS ALL SAMPLES 

TSOR+-/SOFR 

a CALCULATE FINAL 3 STATISTIC 

3* (TSORx ( 12+ (Nx(N±l J ) ) )-(3x(/7+l ) ) 

a RECALCULATE H WITH CORRECTION FOR TIES 
A*TIESK E 

NUM+(*/(A+3 ))-( + AU 
DENOM+Nx((N* 2)-l) 

H+H* (1- (NUM*DENOM ) ) 

a SYSTEM OF LOCICAL STATEMENTS ENSURE PROPER PROB . IS ACCESSED 
(OFClD <2 MOUTPUT 
(X>4 )/F APPROX 



*(K=3)/IF 
*{ (a /(0F= 



.3 3 3 3 ) )v (03C1] >3 ) ) /FAPPROX 
a/ (OF- .2 1 1 1 ) )/ OUTPUT 
Cl] =2 ) 5/342 



-*-((H=4)a(0F[1] : 
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[71] 

[ 72 ] 
73 

[74 » 

w 7 5 
.76 

77 

78 
[79 
[80 
[31 
[82 
[83 

84 

85 
[86 
[87 
[88 
[89 
[SO 
[91 
[92 

till 
31 

[99, 
LlCO] 

[ici 5 

[1023 

Ll03 j 

[104] 
[1053 
l 136 j 
? 10 7 ] 
L108J 



((a /(.OF- 3 1 1 1) )v(a/( OF = 3 2 1 1 ))v(a/(OF= 3 3 1 1))v(a /(OF- 3 3 2 133 

h 



1 ) ) ) /OUTPUT 

)v(a/?0F= _ 3 3 1))v(a/(OF= 3 3 2))v(a/( 0F= 3 3 3333/P33 
)v(0FCl3=4))/P3i 



3/P41 

-(F=4)/P43 

JF: ■►(OFClD >4 )/ FAPPROX 
♦((A /(OF- 2 1 1 ))v(a/(OP= 3 1 
♦ ( (0FL13 =4 ) a (0FC33 =1 ) 3/P31 
+ a /\0F- 3 2 1 5 ) v ( a/ (0F= 

♦((a/(0F= 3 2 2 ))v(OP[i3=4))/P34 
r CASE APPROPRIATE VARIABLE ACCESS CDF 

P23:P+PMATKW2Ql(N-*)\ ;3 
♦PM 

P3 1 : P+PMATKW3 1 C (M- 5 3 : : 3 
♦PM 

P33 :P*PMATKU33 C (M-5 ) ; ; 3 
♦PM 

P34:P-ePMA2’KM34[(M-63; -.3 
♦PM 

P4 1 : P*PMATKW^ 1 C (M ~ 5 3 • ; 3 
♦PM 

P4 2 :P*PMTFW42 C (M- 5 ) ; ; 3 
♦PM 

P43 :P*PMATKWU3 C (/V-7 3 ; ; 3 

r M£I INTER? TO CALCULATE P-VALUE BE INTERPOLATION 
PM:PVAL~3 INTER P ? 

*(PVAL= 1 )/ OUTPUT 
♦E 5 

R CALCULATE P-VALUE USING THE F DIST W/ONE LESS D .F .IN DENOM APPROX . 
FAPPROX : F* ( (N-K)xR 3* ( (X-l )x ( (M-l )-H 3 3 
PWiE+l-(((X-13,((M-X3-l33 FDISTN F3 
♦E 5 

OUTPUT :PVAL * 1 GREATER T3AN .25' 

L3:'TBE H STATISTIC EQUALS: ' . (4®S ) .UTCNL 

'TEE P-VALUE FOR HQ: TEE POPULATION MEDIANS ARE EQUAL VERSUS '.UTCNL 
' 31: AT LEAST TWO POPULATION MEDIANS ARE NOT EQUAL 1ST ’,(4 s?VAL), 

UTCNL 

1 PRESS ENTER WREN READ I . ' 

AA+a 

♦/VI 

El: ’ERROR: ENTER A SINGLE INTEGER VALUE GREATER TRAN 2 : TRI AGAIN .UTCNL 
♦SI 
V 



I 

i 

.11. 

. 12 , 

13 

>: 

31 : 

:io; 

. 21 . 

II: 

;is: 

327 ; 

*30j 

[ 3 2 j 
^33, 
>34, 
.35 
36, 
.37, 
.38. 
H 39. 

40, 

41, 

42, 

43, 
44 J 






7 MANW:N;M:PV2:A;B;C;G\MM:NN;RX;U;NMUP:NU;PVAL;NM;NUMZ;NUMZ1:DEN:DENC;DE 
NCI : TC xTC 1 ; tfCW : Z ; Z 1 ; jILPM ;CDF; INDEX : IPX: Cl ; ; S5 ; CC ; U1 ; U2 ; PV ; NN1 ; NN 

7 ♦ PVT • DIFF • AX • AY • A A • GC* AX\ • PVM • PV7 • D*Q*R 

p to is' FUNchioh uhEs * THfc sum of raMs procedure to calculate tee mann- 

p WHITNEY U STATISTIC WHICH IS USED IN COMPUTING THE P-VALUE FOR THE TEST 
P OF LOCATION AND SCALE. THE C.I. FOR (MY-MX ) , THE SHIFT IN LOCATION IS 
p ALSO COMPUTED. SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE : TIES , TIES 2 
p INDEXPLS , MANWP, INPUT , CONFMW , NORMCDF, AND NORMPTH . 

*N 3 

Ml : MENU MANURE LP 

r MENU CR0ICES~2ND ROUTE TO PROPER STATEMENTS FOR ACTION. 

N3:D*CR0ICEM PACEDMENU MANWOBJ 
DIFF * 0 

♦(0=2,3,4.53/51 

♦(0=6,7.83/53 

+ ( 0=1 J/ifrl 

MSM7 MAINOBJ 
-►0 

B3: J £:iVrFi? THF DIFFERENCE OF THE MEANS OR MEDIANS (MX - MY).' 

DIFF + □ 

+((pDJFF)>l)/E3 

p ENTER DATA VECTORS 

Bl : 'ENTER X DATA ( MORE THAN ONE OBSERVATION IS REQUIRED).' 

N+U 

+(TqoN)=Q)/E 1 
' ENTER Y DATA . ' 

M+ □ 

p IF CALCULATIONS INVOLVE VARIANCES ADJUST X BY THE DIFFERENCE IN MEANS 
N+N-DIFF 

P CONCATENATE X AND Y SAMPLE VECTORS 

A+N ,M 

p DETERMINE SIZE OF X AND Y VECTORS AND ASSIGN TO NN AND MM 

NN+oN 
MM+oM 

r COMPUTE SIZE LIMIT OF LEFT TAIL OF NULL DISTRIBUTION 

NM* (.NNxMM 3 + 2 

mi+imt 

r , ORDER A AND ASSICN TO B 

B+A.IHA1 , 

OUVM.MM 3 INDEXPLS A 

a CALL TIES FUNCTION TO BREAK TIES USING MIDRANK METBOD 

G* 1 TIES B 

r IF FALSE CALCULATE TEST FOR VARIANCES 
♦(0=2,3,4,53/55 

r CALL VARMU TO GENERATE RANKS REQUIRED FOR VARIANCE TEST 

GG+VARMW ( NN+MM 3 

R CALL TIES TO RECORD TIES IN TRE DATA AND BREAK TIES IN CG 



80 



33 



G*GG TIES 3 

a CALCULATE SUM OF X RANKS 

35 :RX*+/ (GZ (CCl; (\NN)1 )] ) „„ 

a CONVERT TO MANNWBIT U STATISTIC 

U*RX-((NNx(NN+ 1))*2) 

IF SIZE OF X TIMES SIZE OF I > 80; GO TO NORMAL APPROX 
*( ( NMx2 )>8 0 )/L2 
AWl+UJW.W? 

a MANWPFUficTIO N CALCULATES LEFT TAIL CUMULATIVE PROBS. OF U STATISTIC 
P*NN2 MANWP NN1 

LOGICAL STATEMENT ENSURES ON LI LEFT SIDE OF NULL DIST IS USED 
*(U£NM 1 )/L 3 

a CONVERT U STAT WREN GREATER TRAN LEFT TAIL VALUES 

U*(NNxMM)-U r „„ 

a IF U 1 IS A FRACTIONAL, INTERPOLATE P VALUE 
L3:*((l\U)=0)/N0N 
U2 *' U 
~(U2>Q )/?l 
PV+1-UPZU2+U )*2) 

*?3 

zm . vy+i - ( (^r,72 ] 2 )* 2 ) 

P3 :~PVI*(PZ‘U2*ll *PZU2+2] )*2 
*CRECK 

N0N:*(U>0)/GO 
PV* 1 
♦P2 

00:?7~1-?01 
P2 :P T /r-*-PC (£7*1 )1 
03303 : - ( U1SNM1 ) /L<* 

■Dyo + D T/ 

>T/l?yr 

? 7 /J^?72 

-Z4 

p I? CONFIDENCE INTERVAL DESIRED CO TO L 10 
£2:-(£=5 )/ZlO 

p COMPUTE IRE NORMAL APPROXIMATION W /CORRECTION FACTOR 
NUMZ+(U+0.o)-NM 
NUMZ1+(U-Q.5)-NM 
DEN+ ( (^Vx/VNx(OT+iViV^i ) ) + i 2 )*G . 5 
Z+NUMZ+DEN 
Zl+NUMZl+DEN 
NUM+\ (U-NM) 

DENC+l ( ( ( NN+MM-1 )x (DEN*2 ) )*(NN+MM-2 ) )- ( ( (NUM-Q . 5 )*2 )+ (NN+MM-2 ) 1 )*0 . 5 
DENCl+ll ( (NN+MM-1 )x (DEN*2 ) )* ( NN+MM-2 ) )- ( ( (M7M+0 . 5 )*2 )* (NN+MM-2 ) ) )*0 . 5 
TC+(NUM-0 .S)*DENC 
TCl+(NUM+0 . 5 )*DENC1 
+(OZNM)/ SECOND 

PVI+UNORMCDF Z)+((NN+MM-2) TDISTN TO) *2 
PV+((l-(NORMCDF Zl))+(1- ((NN+MM-2) TDISTN TC l)))+2 
+L4 

SECOND : PVI+( (NORMCDF Z ) + (l - ( (NN+MM-2 ) TDISTN TC)))* 2 



PV+{ (l-(NORMCDF Zl))+((NN+MM-2) TDISTN TC 1))*2 

. — - ’)) 

_ J1 )/iV5 
P73+1 



L4:P^3^2x<;l/(P7,P7I] 
-►(P73£l)/tf5 



[1051 
[1061 
L107 J 

[1081 

[109] 

1110 j 

[111] 

[ 112 , 

[113, 

[114] 

[115] 



iVS:P7M-*-(3 ,1 )p(P7J,P7.P73 ) 

' TRE SUM OF TRE X RANKS IS: ',(vRX),<. TEE U STATISTIC EQUALS: ’,(vUl) 
OTCNL 

a LOGICAL STATEMENT FOR VARIANCE OUTPUT 

♦0=6,7, 8 )/VAR ' 

1 TRE P-VALUE FOR BO: MX = MI VERSUS B 1: MX ' , (vLOCICZD-l ; 11 ) , ' MI IS: 1 
(6 u 9PVM CP- 1 ; 1 ] ) , DTCNZ 
*L 8 

VAR : PVM * ( 3 . 1 ) p (PV , PVI , P73 ) 

1 TRE P-VAtUE FOR 30: VX = VI VERSUS Hi: VX ' , (vLOCICZD-S ;1] ) , 1 VI IS: ' 

(6 u SP7MC3-5; 1 ] ) ,QTCNZ 

1 PRESS ENTER WREN RE ADI . 1 
BB+a 
*N3 

L 8 : 1 WOULD I0U LIKE A CONFIDENCE INTERVAL FOR TEE SEIFT IN LOCATION (.MI - 
MX)? (I/N). ' .OTCNL 
B3-a 

— ( 33-' IV )/N 3 
310 : CC*INPUT 3 
ALPS A * ( 1 0 0 -00 ) * 2 00 

a ROUTE TO NORMAL APPROX. FOR CONF . JiVr. OP LARGER SAMPLE SIZES 
+((NMx 2)>80 )/L 5 

a COMPUTING CONFIDENCE INTERVALS BI EXACT P-VALUE 

CDF*P 

a INDEX POSITION OF VALUE IN CDF <. ALPRA 

INDEX * ( +/ ( CDFZALPRA ) ) 

*(INDEX> 0 )/£6 
INDEX* 1 

♦i 6 

fl COMPUTING CONFIDENCE INTERVALS USING NORMAL APPROX. W/C.F. 

L5:DEN*((MMxNNx(MM+NN+l ) )*12)*0. 5 
C/.4LPHA-*- (DEN* (NORMPTH ALPRA ) ) +NM- 0 . 5 
« 30 0W0 U ALPRA DOWN AND INCREMENT BI ONE 
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; 1 3 2 ] INDEX* L U ALPHA + 1 
.133] L 6 : IPX*NN .INDEX 
‘ " CI+IPX cdNFMW A 

» a » , (wco , 1 ssffidffsf liviFFVAi fsf the shift in ' 

'LSSAIISW BETWEEN POPULATIONS X AND I IS:',QTCNL 
;i37] * ( ' (5CIC1]),» S MY - MX Z »,(*SI C2]), ! )',0ISFL 

13 81 1 PRESS ENTER WHEN READY . ' 

139J SF«-Q * 

1401 -WV3 ° 

141 1 FI : ' ERROR : SAMFLF CONTAINS LESS THAN TWO ENTRIES ; M AGAIN ,UTCNL 
.142J -►SI 

14 3 3 F3 :'FFFSF: YSF FAFF ENTERED MORE THAN ONE VALUE. TRY AGAIN .' t UTCNL 

.144] *3 3 

7 



ASSICN THE SIZE OF X (AND Y) TO N 
COMPUTE THE SUM OF X' S AiVD Y'S 

COMPUTE THE MEAN OF X AND Y 

K COMPUTE THE SUM OF THE X' S SQUARED 
COMPUTE THE SUM OF X TIMES Y 



7 NPLR : N ; SFMX : SFMY ; XBAR : YFAF : SUMX2 ; SFMXY ; B ; A ; JW ; XX ; BB ; a : D ; ALPHA ;P;CC ; CDF ; 
T ALPHA ; NN ; Cl ; SLOPES ; RR : SR : 77 : DENOM : 7/VDFX ;FF:X:Y:Q;R; CHA ;E;PV 
a PFSCFAM CONDUCTS NONPARAMETRIC LINEAR REGRESSION. THE LEAST SQUARES 
A ESTIMATED REGRESSION LINE IS COMPUTED WITH HYPOTHESIS TESTING AND 
a CONFIDENCE INTERVAL AVAILABLE FOR THE SLOPE 3. IF 3 DOES NOT LIE IN THE 
a C.I. AN ALTERNATE REGRESSION LINE IS PROPOSED. SUBPROGRAMS CALLED ARE 
a SPMANP , KENDALP , NORMPTH , INPUT, AND CONFLR. 

DPP* 5 

a DISPLAY MFJVF AFD 7FPFI DAIA • 

Nl'.E+MENU NPLROBJ 
+(E=1)/N2 
MENU MAIN03J 
-*•0 

N2 \R* INPUT 2 
Q+l+R 

Z*l±( Q + D+R 
Y*(Q + 1 HP 
a 

N*pX 

*SUMX*+/X 
SUMY++/Y 

XBAR+SUMX*N 
YBAR+SUMY+N 

SUMX2*+/(X*2Y 

SUMXY*+/(X*Y ) 

A COMPUTE 'S' IFF SLSPF SF IFF FSIIMAIFD LFASI SCFAFFS REGRESSION LINE 
B+( (NxSUMXY)- (SUMXxSUMY) H ( (NxSUMX 2 )-(SUMX* 2 > ) 

A COMPUTE ' A ' , IFF Y- INTERCEPT 

A+YBAR- (BxXBAR ) 

FF«-»JV' 

'IFF LFASI SQ FAFFS FSIIMAIFD REGRESSION EQUATION IS:',UTCNL 
'Y - ',(5*A),» + ' (5*F). 'X. ' DISFL 

'DO YSF J/ISF IS FiVIFF SSMF X /ALFFS IS SFI IFF PREDICTED Y"S? (Y/N).' 

WW+ □ 

+ (WW='N' )/L 1 
12: 'ENTER X VALUES. ' 

XX*U 

A CALCULATE PREDICTED Y'S 

YY*A+BxXX 

' IFF PPFDISIFD Y 7ALFFS AFF : ' . ( 9 YY ) , UTCNL 
'WOULD YOU LIKE TO RUN SOME MORE X VALUES? (Y/N).' 

WW + Q 

*(WW-'Y' )/L2 

Li: 'WOULD YOU LIKE TO TEST HYPOTHESIS ON B , IFF SLOPE OF THE EQUATION? (Y / 
N ) . ' 

ww + a 

-►(JW='F» )/£3 

'ENTER THE HYPOTHESIZED SLOPE. ' 

BB*U 

a COMPUTE UI'S 

U*Y-(BBxX) 

CHA*'> ' 

a CALL SPMANP TO COMPUTE R 
D*X SPMANP U 
~(DZll>0)/Lll 
CHA— ( < 1 
Lll:PV*2xDZ2l 
+(PV£1 )/Ll 9 
P7«-l 

£19: 'SPEARMAN' 'S R FQFALS: 



AND ASSSSIAIFD P-7ALFFS. 



1 , (4®DC1 ] ) .UTCNL 



_ rtUkJ. | \ / • LJ jl UJH 

' TEE P -VALUE FOR SO: B = ',(955),' 7F.RS5S' 

' 51: 5 ' , (.bCBA [1 2]7Tl9BB).' IS: 



' . (495C23 ) .aicwi 



' TBE P-VALUE FOR TBE TWO-SIDED TEST bF HIPblBESIS Is: ' , (49 ( 2 xDC 2 ] ) ) ,UTCN 

L 

a IF USING TBE NEW REGRESSION EQUATION BASED ON MEDIANS, EXIT BE RE. 
L3:+(FF='N' )/L 18 
' PRESS ENTER WEEN READ!.' 

WW+0 
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[6 9] 
[7o; 

C 71 

72 

I 73 * 

.7o, 

76 

77 
[78< 

1 79 j 

,82; 

[83 

84 

35 

[ 86 , 

[87 

L 8 8 , 

[ 89 j 
C50l 

C 9 1 

,32. 

i93 

94 

las; 

96 
[97, 
[9 8, 
[39; 
.130^ 

1.1013 

il02' 

l!03] 

Clou 

[105 

[106 

[107 

[108 

[109 

,110 

11*? 

[H2 

[113 

[ 1 1 4 

Cll5 
[116 
117 
,118 
.119 
[120 
[121 
[122 
123 
[124 
.125 
L126 
[127 
,128 
.129 
[ 13 0 
.131 
[132 
[133 

[134] 

[135 

[136. 

[137 

[138, 

[139 

[140] 

L 141 



-►2V 1 

p COMP UTS CONFIDENCE INTERVALS ON S x 

LIS : ‘ WOULD IOU LIKE A CONFIDENCE INTERVAL FOR TEE' SLOPE? (J/2V).' 

ww * a 

-►(W= , 2V» )/N 1 
£10 ICC* INPUT 5 

o CHANGE ENTERED VALUE TO ALPEA 

ALPHA+{\QQ-CC)*2Q0 

a ROUTE TO NORMAL APPROX . FOR CONF . INT. OF LARGER SAMPLE SIZES 
* (2V>13 )/£5 

p COMPUTING CONFIDENCE INTERVALS 8X EXACT P -VALUE 

P+KSNDALP N 
CD F+P 12 :] 

p iiVIFX POSITION OF VALUE IN CDF S ALPEA 

INDEX * ( +/ {CDFSALPEA ) ) 

-(JiV£FX>0)/£6 
INDEX* 1 
-►£6 

p COMPUTING CONFIDENCE INTERVALS USING NORMAL APPROX . W/C.F . 

L 5 : DENOM* ( (2Vx (2V-1 )x ( (2xtf ) + 5 ) )*18 )*0 . 5 
r^IFS^OFiVOMx ( | ( NORMPTE ALPEA ) ) 

£6^£FF;i«-?[3 ;I2V£FX] 

3 i C± CON P +jR I 

NN+1+CI 

SLOPES+l+CI 

RR+l ( (NN-TALPHA )*2 ) 

RR 

*(RR-0 )/L2 0 

£20 (l+( (2V2V-T<;£PF<1 )*2 ) ) 

SR 

* SRsiz SLOPES) )/L21 
SR^-z SLOPES 

L2l:'A ' («CC),» CONFIDENCE INTERVAL FOF 5. IFF S£0PF OF f 
' IFF FSIIMIFP REGRESSION LINE , IS: ' DrC7V£ % T „ 

' ( ' (5 SSLOPZSZRR2),' 2 3s * , (5*S£0PSSCS3] ) , ' K'.DX’Ctfi 

'PRESS ENTER WHEN READ X. 1 
wv *-& 

o I? 3 OUTSIDE TEE C.I. CALCULATE NEW EQUATION BASED ON MEDIANS 
-(( SZSLOPESiRR 2 )a (3ZSL0PESZ3R1 ) )/Nl 
p ORDER X AND I 

x-<-x[<wn 

i+jCaj] 

p CHECK TO SEE IF THE SIZE OF SS IS EVEN OR ODD FOR FINDING MEDIANS 
♦( C 2 T JViV ) = 0 )/Sl 

o COMPUTE MEDIAN FOR ODD CASE 

B+SLOPESZ(.(NN+l)*2)l 
*S 2 

p COMPUTE MEDIAN FOR EVEN CASE 

SI : B* ( SLOPES C (AW* 2 ) ] + SLOPES C ( ( JVJV+2 ) * 2 ) ] ) + 2 
p 30 2*33 S^MB 303 2*33 X A3D X VECTORS 

S2:*((2I3)=0)/S3 
XBAB+XC ((3+1)t2)] 

XB4B-*-X[ ( CJV+1 5 + 2 )] 

*OUT 

S3 : YBAR+(.Yl(N*2 )] +XC ( (3+2 ) + 2 )] ) + 2 
XB«l3-*-(X[ (3*2 )] +XC ( (3+2 )*2 )] )*2 
p COMPUTE NEW INTERCEPT 'A' 

OUT A+1BAR- (,B*XBAR ) 

'THE LEAST SQUARES ESTIMATOR OF B LIES OUTSIDE THE CONFIDENCE INTERVAL.' 
'DISCARD THE LEAST SQUARES EQUATION AND USE: < .OTCNL 
' X = ',(®4J ' + ' (®B) ,'X OTCNL 

' THIS EQUATION IS BASED ON THE MEDIANS OF THE X AND X DAT A, AND TH 

E * 

1 MEDIAN OF THE TWO -POINT SLOPES CALCULATED FOR THE CONFIDENCE INTERVAL 0 
N B * 0 TCNL 

p "aLlOW USER TO DO SOME ANALYSIS ON NEW EQUATION 
' DO YOU WISH TO ENTER SOME X VALUES TO GET PREDICTED Y"S ' 

'FROM THE NEW EQUATION ? (.Y/N).' 

FF-a 

*► ( FF = 1 7 1 )/L 2 
FF+* I 1 
-►£1 
7 



[ 10 ] 



7 SIGN ; A ;C : B : D :PV AL :X :MO ; N : CDF ; ALPHA : Cl : I : A A ; BB ; CC ; DD ; PV ; PVI ; NNN ; KPOS ; ORD 
D ; II : 0kDk : kklPRA : i : il :QUh :WW :PVM :Pfr3 ik ; 6 

THIS FUNCTION USES THE ORDINARI SIGN TEST TO CALCULATE THE K 
STATISTIC . P -VALUE. AND CONFIDENCE INTERVAL AS A TEST FOR MEDIANS . 

THE LAST OPTION WILL DISPLAI A TABLE OF CONFIDENCE INTERVALS OF ORDERED 



STATISTICS WITH CONFIDENCE COEFFICIENTS . 

SUBPROGRAMS CALLED BI THIS FUNCTION INCLUDE X BINOM , NORMCDF , NORMPTH , 
INPUT, AND QUANC. 

-JV4 

NUMENO SIGNHELP 

p MENU^ CHOICES AND ROUTE FOR PROPER ACTIONS 

2V4 : C+CHOICES PAGEDMENU SIGNQBi 
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PAIRED SAMPLE CASE 



COMPRESS D TO REMOVE ZEROS 
RECORD LENGTH OF A AND ASSIGN TO N 
KEEPING TRACK OF POSITIVE SIGNS 



*(C=2 ,3,4,5 )/£8 
*lC= 7,8,9,10 )/£ 9 
-*-(C=l)/JVl 
*(C= 11 )/L 20 
MENU MAINOBJ 

a*° INPUT DATA FOR SINGLE SAMPLE CASE 

L8:AA*1 
X* INPUT 1 
NNN*oX 

-*• (C= 5 )/£16 
MO* INPUT 3 
D*X-MO 
*L 11 

L 9 :AA*2 
R+INPUT 2 
Q*liR 

X*l + (Q+l )+U 
Y*(Q+1)+R 
DD*X-Y 
NNN*oDD 

->-(C=10 )/£16 
MO* INPUT u. 

D*(X-Y)-MO 

Lll:A+(Dxo)/D 

p 

N*pA 

EPOS**/ ZA> 0 ) 

*(N>25)/NORM 
PVAL*BINOM N 
*(KPOS>OWPl 
PVI-1 
*P2 

PI : PVI+1 -PVAL ZKPOS] 

P2:PV*PVALZ(XPOS+l)l 
*L 6 

p IF N IS GREATER THAN 30 USE NORMAL APPROX W/ CONTINUITI CORRECTION 
NORM: Zl*( (KPOS-O.S)-(O.SxN) )* (C . 5* (N*0 . 5 ) ) 

Z*( (KPOS*C . 5 ) - (0 . SxJV) ) + (0 . 5x (iV*o . 5 ) ) 

PV+NORMCD? Z 
PVI*1- (NORMCDF Zl) 

p IF PAIRED SAMPLE TEST GO TO L 17 FOR OUTPUT STATEMENT 

Lo:PV3+2x(l/ (PV ,PVI)) 

*(PV3£1)/NS 
PV 3*1 

NS :PVM*( 3,1 )o ( PV.PVI.PV3 ) 

'COMPUTATIONS ARE BASED ON A SAMPLE SIZE OF: < . (*N) , ' . < .DTCNL 
'THE TOTAL NUMBER OF POSITIVE SIGNS IS: < , (.9KP0S ) , DTCNL 
-*-(C=7,8,9)/£17 

'THE P-VALUE FOR HO: M = ' .(*MO).< VERSUS HI: M ' . (9LOGICL (C-l ) : 1] ) . ' 

9M0),< IS: ' , (5®PraCC-l;l ] ), DTCNL 

-*■£18 

Lll'.'THE P-VALUE FOR COMPARING THE MEDIAN OF THE POPULATION OF ' 
'DIFFERENCES TO THE HYPOTHESIZED MEDIAN. '.DTCNL 

<H 0: M(X-Y) = '.(9MO),' VERSUS Hi: M(X-Y ) 1 , (vLOGICC (C-6 ) ;1] ) , ' ' ' 

, IS: ' ,(S9PVMZlC-6);l] )7DTCNZ 

L18: 'MOULD YOU LIKE A CONFIDENCE INTERVAL FOR THE MEDIAN? (Y/N). ' .DTCNL 

BB*D 

*(BB='Y ' )/£ 16 
+QUANT 

a INPUT SIZE OF CONFIDENCE INTERVAL 

£16 :CC*INPUT 5 
ALPHA* ( 100-CO + 20 0 
*(NNN>2S)/N0RM1 

p COMPUTING CONFIDENCE INTERVALS BY EXACT P-VALUE 

CDF*BINOM NNN 

P INDEX POSITION OF CDF FOR ALPHA * 2 

B*+/ (CDFiALPBA ) 

* (B>0) / SKIP 
B* 1 
*SKIP 

p COMPUTING CONFIDENCE INTERVALS BY NORMAL APPROX. 

NORM1 : KALPRA* ( ( 0 . 5 x (NNN* 0 . 5 ) )x ( NORMPTH ALPHA ) ) + ( 0 . SxNNN ) - 0 . 5 
5 ROUND KALPHA DOWN TO NEAREST INTEGER AND INCREMENT 3Y ONE 
3*>. RALPH A + l 

a IF SINGLE SAMPLE CASE CO TO LI 

SKIP:*(C=2 ,3,4,5 )/£7 

p CALCULATE AND PRINT OUT CONF. INT . FOR PAIRED SAMPLE CASE 

L5:ORDD*DDZ&DDl 

YY*f>ORDD 

CI*0RDDZB1 .ORDDZ (JJ-(fl-l))] 

' A ' ,(9CC)' CONFIDENCE INTERVAL FOR THE MEDIAN OF THE < 

< POPULATION OF DIFFERENCES IS:' , DTCNL 

' ( ' , (®CT[1] ), 1 Z MEDIAN (X-Y) Z ',<®CT C2] ) , ' )' .DTCNL 

*QUANT 

p CALCULATE AND PRINT OUT CONF. INT. FOR ONE SAMPLE CASE 

L7:ORDX*XZ£in 
YY*oORDX 
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'1001 CI+ORDXZB] .ORDXZiYY- (5-1 ) ) 1 

!l 01 j 'yl ' (see),' CONFIDENCE INTERVAL FOR TEE MEDIAN OF TEE POPULATION IS: ' 

0 TCNL 

■ 102 ] f ( '.(serci]),' S MEDIAN Z ' (seJC2]),' )', □ TCNL 

'103' QUANT: 'WOULD YOU LIKE CONFIDENCE INTERVALS FOR 1 SPECIFIED QUANTILE? (.Y/N 
) . 1 , QTCNL 
WW+Q 

)/B 1 

*N 4 

L20:' ENTER TEE SIZE OF TEE SAMPLE . ' 

NNN+C 

31: 'ENTER DESIRED QUANTILE-, FOR EXAMPLE: ENTER 20, FOR THE 20 TE QUANTILE. 
QUA* □ 

-*-( (<3M<l0)v(iJM>100))/51 

NNN QUANC QUA 

t I 

» ***** THIS TABLE GIVES CONFIDENCE COEFFICIENTS FOR VARIOUS INTERVALS ' 
'WITH ORDER STATISTICS AS END POINTS FOR THE ',(*QUA),'TH QUANTILE. ' ,DTC 
NL 

' PRESS ENTER WHEN READY.' 

33+iL 

WV4 

,.. 3 , Eli' ERROR: TEE QUANTILE VALUE MUST LIE BETWEEN 0 AND 10 0; TRY AGAIN.' 

Cl20] -31 



7 SPMA N : X : Y ; A : 0 : R : CH A : 3 3 ; 3 ; P V 

a THIS FUNCTION COMPUTES TEE SPEARMAN R STATISTIC WHICH MEASURES 

a TEE DEGREE OF CORRESPONDENCE BETWEEN RANKINGS OF TWO SAMPLES. THE ?- 
a VALUE IS GIVEN FOR TESTING ONE AND TWO-SIDED HYPOTHESIS OF ASSOCIATION 
a SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: TIES , IIESK, SPEARP , 
a INPUT , SPAPROX, INTERP , AND THE VARIABLE PMATSP . 

a DISPLAY MENU AND INPUT DATA. 

N1 : 3+MENU SPMANOBJ 
-(£=1)/S1 
MENU MAI NOB I 
-*0 

Bl:R+INPUT 2 
y«- (§+i )+r 

o CALL SPMANP TO CALCULATE THE STATISTIC AND ASSOCIATED P -VALUES 
A+X SPMANP Y 
+ (Al 1]>0 ) /LI 
CHA+' INDIRECT ' 

+L2 

L 1 ; CHA + 1 DIRECT ' 

L2:PV+2*AZ2l 
+(PVZ1)/L3 
PV + 1 

L3:' SPEARMAN "S R EQUALS: ' , (**A [1 ] ) , UTCNL 
'THE P-VALUE FOR HO: NO ASSOCIATION EXISTS VERSUS' 

' Hi: '.(*CHA)' ASSOCIATION~EXISTS IS: 1 ,(4®A[2l) 

'THE P-VALUE FOR THE TWO-SIDED TEST OF HYPOTHESIS IS: *,(4*(2 *A 
L 

' PRESS ENTER WHEN READY. ' 

B3 + Q 
+N 1 
V 



\1^ T ) C ,utcn 



V WISIG;AiB:D:E:F:PV2;Zl ;Z:DEN;NUMZ:NUMZl :PVAL\X:MO \N xTPLUS: CDF : T ALPHA :AL 
PH A ; HyCI ; Y : A A ; 33 ; CC:NN:Du ; PV ; POS ; TP OS ; NM ; PVI ; TPOS 1 ; NNN ; C ; PVM ;PV3;R;Q;NU 



M;TC:Tc 1 ]rkAP\ DENT :bENtl 

this Function uses the wil 



a thIs Function Uses The wilcoxon signed rank test to calculate the tplus 
a statistic, p-value, and confidence interval as a test for medians. 

a SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: TIES , WILP , NORMCDF , 
a NORMPTH. CONFW , INPUT AND THE VARIABLE ? MATRIX . 

NUMENU WILHELP 

a 'MENU CHOICES AND ROUTE FOR PROPER ACTIONS. 

NV.C+CHOICEW PAGEDMENU WIL03J 
-► (C=2 ,3,u,5 )/L3 
-(C=6.7 ,3,9)/£9 
+(C=l)/Nl 
MENU MAINOBJ 
+0 

LS-.AA+l 
X+ INPUT 1 
NNN+oX 
+(C=5)/L1S 
MO+INPUT 3 



INPUT DATA FOR SINGLE SAMPLE CASE 
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PAIRED SAMPLE CASE 



[20 

,21 

[22 

,23 

,24 

25 

26 
27 

.28 

"29 

,30 

I- 31 

[32 

[33 

[34 

[35 

.36 

[37„ 

if 5 

IL'J 

J 43J 
44 

45 
46 
47 
48, 
L 4 9 

tf ? 4 

kli; 

C 54 , 

t 55 

[56 

[57 

Css 

[59 

[60 

[61 

.62 

.63 

[64 

65 

66 

.67 

.68 

[69 

70 

71 

72 

.73 

.74 

.75 

76 

[77 

[78 

79 

[80 

[81 

[82 

,83 

[84 

[85 

[86 

[87 

.88 

89 

.90 

91 

[92 

[93 

,94 

:-22^ 

[98 

.99 

.100 

.101 

.102 

1 103 

[104 
f 10 5 
[ 106 
1107 



D+X-MO 
+L 11 

L9 :AA + 2 
R+INPUT 2 

X+l+CSt 1 ^ 

J«-(Q+1 )+£ 

£D«-X-I 

NNN+oDD 

-(C=9)/L16 

MO+INPUT 4 

D«-(X-I)-M0 

Lll:A+O>sO)/0 

p 

fl ?05^U>0) 



COMPRESS D TO REMOVE ZEROS 
RECORD LENGTH OF A AND ASSIGN TO N 
KEEPING TRACK OF POSITIVE SIGNS 



p TAKE THE A3 SOLUTE VALUE OF A; ASSIGN TO 3 AND ORDER 3 

B+\A 
3+BI&B1 

o REORDER POSITIVE SIGNS TO COINCIDE WITH PROPER POSITIONS IN 3 

POS+POSUC\Ajj 

p CALL FUNCTION TO BREAK TIES 

E+l TIES B 

a CALCULATE TPOS BY ADDING ACROSS ALL POSITIVE VALUES OF E 

TPOS-*-*/ (POSxE) 

TPOSl*TPOS 

a GIVES SIZE OF LEFT TAIL OF PROBABILITY DISTRIBUTION 

NM+(U(*/\N)*2 ) )*1 

a GO TO STATEMENTS 3ASED ON LENGTB OF VECTOR 3 

~(N><3)/L 3 

a CENERATE NULL DISTRIBUTION FOR TPLUS 

F*WILP N 

a IF TPOS FALLS IN LEFT HALF OF PROB DIST CALCULATE PVALUE AS NORMAL 
LX: + (TPOS<, ( NM-X))/TP.EG 

a OTRERWISE USE TEE NEGATIVE T STATISTIC 

TPOS-*- ( + / uV ) -TPOS 

a IF TPOS IS FRACTIONAL USE BOTE TEE INTEGER ABOVE AND BELOW AS TPLUS 
TREG TPOS } = 0 ) /NON 
TPLUS+ITPOS 

Pl:PK+l-( (FtTPLUS)+F ZTPLUS+ll )*2) 

PVI-*-(FlTPLUS*l)*FZTPLUS*2l )*2 
r v 

NON:+(TPOS>0)/CO 
PV* 1 
-*-P2 

go:pv*x-fz(tpos)\ 

P2 : PVI*F C (TPOS +1)3 
CHECK : ♦ ( TPOS1 Z (NM-1))/L6 
PV2+PV 
PV-*-PVI 
PVI+-PV2 
*L 6 

a COMPUTE NORMAL APPROX . W /CONTINUITY CORRECTION FACTOR 
L2 :TRAP*(Nx (N+l j ) + 4 
NUMZ*(TPOS+0 . 5 )-TRAP 
NUMZ 1 * ( TPOS -0.5) -TRAP 
DEN*( (Nx (N*l)}*( (2xN ) + l)) + 24)*0.5 
Z*NUMZ*DEN 
ZX*NUMZX*DEN 

a COMPUTE STUDENT T APPROXIMATION WITS CONTINUITY CORRECTION FACTOR 
NUM* 1 (TPOS -TRAP) 

DENT*(^Nx(DEN*2))*(N-l))-((,(,NUM-0.S)*2)*(N-l)))*0.S 



... D)-((CiVra-0.5)*2)+(iV-l)))v _ 

. . (.NX (DEN*2))* (N-l ))-( ( (NUM* 0 . 5 )*2 )* (N-l ) ) )*0 . 5 

TC*(NUM-0 . 5 )*DENT 



DENT1- 



TC1*(NUM*0 . 5 )*DENT 1 
A COMPUTE AVERAGE OF TC AND ZC 
*(TPOS£ ((*/ \N)*2) )/ SECOND 

PV*((l-(NORMCDF Zl))+(l-( (N-l ) TDISTN TC 1)))*2 

pvi*Unormcdf z)*X(n- 1 ) - 

*L 6 



TDISTN TC ))* 2 



SECOND :PV*( (1- (NQRMCDF Z1))*((N-X) TDISTN TC1))*2 
~?VI*((NORMCDF Z)*(l-((N-1) TDIS — 



'<6 : PV 3-2 x ( L / (PV , PVI ) ) 
:x)/ns 



TDISTN TC)))*2 



?V2*1 

N 5 : PVM* ( 3 , 1 ) p ( PVI , P V . PV3 ) 

1 COMPUTATIONS ARE BASED ON A SAMPLE SIZE OF ' .(fN) ,< .< .DTCNL 
'TEE TOTAL SUM OF POSITIVE RANKS IS: ' . (vTPOSl ) , ' . ’ . DTCNL 

a / IF PAIRED SAMPLE TEST GO TO L 17 FOR OUTPUT STATEMENT 

*(AA=2)/L17 

'THE P-VALUE FOR E0: M = '.(sMO),' VERSUS EX: M < , (nLOGIClC-1 ; 1] ) , ' ' , (u> 
MO)' IS:' . (4*praCC-l;l) ) ,i)TCNL 
*L 18 

L 17:' THE P-VALUE FOR COMPARING TEE MEDIAN OF TEE POPULATION OF ' 

' DIFFERENCES TO THE HYPOTHESIZED MEDIAN. ' , DTCNL 

H1: MU ~ h '-^WCICLCS-.U).' ',( *MO).'. 
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[103] 

[109] 

110 

'111 

1112, 

[ll3, 

,114, 

.115. 

116. 

[117 

.ns; 

,113, 
]l20 
121 
L 12 2 , 
[123, 
124, 

r 1 2 5 " 



L128J 
[ 129 1 

■13 0' 

:iai3 

|l32] 

133] 

,13U 

135, 

:i3§; 

*133 1 

[1231 

1^0 

i'IUI 

L1423 

riu 3 i 
r * uu* 
1145] 
L146] 



Ll 8: 'WOULD YOU LIKE' A * CONFIDENCE INTERVAL FOR TEE MEDIAN? ( Z/N ).' 
+(BB='Z' )/L16 

-*Nu> c ' 

LIS : CC*-INPUT 5 
ALPHA+ (10 Q-CC )t200 

p ROUTE TO NORMAL APPROX. FOR CONF INT OF LARGE SAMPLE SIZE 

+(NNN>9 )/£4 
CDF+WILP NNN 

p INDEX POSITION OF CDF FOR ALPHA + 2 

CDF+(CDF*0)/CDF 
T ALPHA*- ( + / ( CPPSJLPffA ) ) 

+(TALPHA*0)/JUMP 

T ALPHA *- 1 

+JUMP 

p COMPUTING CONFIDENCE INTERVALS BZ NORMAL APPROX. W/C.F. 

Lu : TP;LP+ ( MWx ( AWiV+ 1 ) ) + 4 
DEN+( (NNN* (NNN+1 )x ( (2*NNN )+l ) )*24 )*0 . 5 
TALPHA*- (DEN* ( NORMPTH ALPHA ) )+MP- 0 . 5 
P^PPPa 

p ROUND TALPHA DOWN TO INTEGER VALUE AND INCREMENT 31 ONE 

TALPHA+lTALPHA+1 

a I? ONE SAMPLE CASE CO TO LI 

JUMP: + (C=2 ,2,4,5 )/L7 

fl CALCULATE AND PRINT OUT CONF. INT. FOR PAIRED SAMPLE CASE 

L5 :CI+T ALPHA CONFW DD 

■ A '.(ICC) » CONFIDENCE INTERVAL FOR THE MEDIAN OF THE * 

'PCPCLAriW CP DIFFERENCES IS:' ,QTCNL 

' ( ' , (sCJCl] ), » Z MEDIAN (X-Z) Z » , (*CT C2]), r ) » , Drew 

' PPPSS ENTER WHEN READZ . 1 
E5+Z 
-*N 4 

p juto ppjivr out conf . 

Ll iCI+T ALPHA CONFW X 
'A ' (®CC) f r CONFIDENCE INTERVAL FOR THE MEDIAN OF THE POPULATION IS:' 
OTCNL 

( » <®CTCl2), r Z MEDIAN C, »,(®CJC 2]),' ) * , CrCtfL 
'PRESS ENTER WHEN READY. ' 

33+L 
-*N 4 
7 



:/VP. PCP C/VP SAJ*PPP CAPP 
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APPENDIX E 



MAIN PROGRAM LISTINGS FOR MAINFRAME COMPUTER WORKSPACE 



[l, 
[ 2 . 
"3; 

71 

!8J 

*9] 

[10; 

.li. 

[12, 

[i3, 

!S= 

L -L / , 

C 18 

L19, 
[20 
L 21 
[22 
[23 

'2U 

[25 

[26 

27 

28 
.29 

30 

31 

32 

33 
.34 
.35 

36 

37 
.38 
.39 
.40 
.41 
42 
.43 
.44 
.45 
[46 

5S 

kl 

51 

£52: 

,53, 

:ls; 

hi 

.ol 

[62 

[63 

[64 

.65 

66 

.67 

68 

.69 

.70 

71 

[ft 



7 KENDALL ;A:AAiB;BX;BYiC:CXiCY;D;DD;DX;DY:DXY;S;POS;NEG}XX;YY\N;DEN;NN;NU 
M:P :PV:PvAl :&U ; SV : T : (j 5 1 ; At ; Z : X : X : CM s Q ; A 
A fois PUNCTtON COkPUTES THE kEtiD&LL B STATISTIC WHICH IS A MEASURE 
a OF ASSOCIATION BETWEEN SAMPLES . P-V ALVES ARE GIVEN FOR TESTING ONE 
a AND TWO -SIDED HYPOTHESIS FOR NO ASSOCIATION VERSUS ASSOCIATION . 
q SUBPROGRAMS CALLED 31 THIS FUNCTION INCLUDE : TIES , IIESK, KEN D ALP , 
q INTERP INPUT AND NORMCDF. 

R+INPUT 2 



3<-i+c 



.. _/(Q+l)*J? 

x«-(q+i)+p 

VxCAX] 

P+XCAX3 

* 



OMM X INCREASING ORDER OF X 
ORDER X JJV INCREASING ORDER 
COMPUTE CURRENT RANKING OF I 
NOW ORDER I RANKS IN INCREASING ORDER 



D+Al&Al 

A IF TIES EXIST IN EITHER X OR I RANKED VECTOR USE MID -RANK METHOD 

DD + 1 TIES D 
XX+1 TI r S 3 

A FIND ORIGINAL RANKING OF I WITH TIES RESOLVED 

YY+DDCCl 
N + cX 

q COMPUTE NUMBER OF DISTINGUISHABLE PAIRS 
MMtfx(tf-l) ) + 2 
S-«-pO 
AA + 0 

A POSIXJ7P CtfPS CCMP FROM A RUNS UP CONDITION : NEGATIVE 1 FROM RUNS DOWN 

A ZERO IS SCORED FOR TIES . MULTIPLY THE RESULTS FOR EVERY ELEMENT AND SUM 

LliAA+AA+1 

bx+(,xxZaai>(aa+xx}) 

CX+(XXtAAUUA+XX)U( 1) 

DX+BX+CX 



BY+(YYlAAl>(AA+YY)) 

aaUU' ' 

CX 

DXY+DXxDY 



CX+(XX_„. 
PX+PX+CX 



(M +XX ) )x (“1 ) 



POS+(PXX>0) 

neg+(dxy< 0 )x( 1 ) 

S+S t POS,NEG 

+(A&<(N-1))/Ll 

fi SUM FINAL VECTOR TO DETERMINE S 

S++/S 

A OBTAIN TEE NUMBER OF TIES IN EACE VECTOR USING TEE TIESK FUNCTION 
U*TIESK B 
V+TIESK D 

SU*+/(2'.U) 

SV*+/(21V) 

A CALCULATE TEE B STATISTIC INCLUDING TEE CORRECTION FOR TIES 
T-*-S*((NN-SU)x(NN-SV))*O.S 
AT-*- 1 T 

*(N>12)/NORM 

fi CALL KENDALP TO CALCULATE TEE RIGET TAIL OF TEE CDF OF B 

P+KENDALP N 

a CALL INTERP TO CALCULATE P-VALUE BI INTERPOLATION 
PVAL+AT INTERP P 
+ (PVAL* D/LB 
PVAL+ 0 . 5 

*£3 

a CALCULATE P VALUE USINC NORMAL APPROX. 

NORM : MUM* (. 3 xAT) * ( ( 2 x NN ) * 0 . 5 ) 

DEN+(2x((2*N)+S ) )*0. 5 
Z*NUM*DEN 

PVAL+l- (NORMCDF Z) 

fi IF B IS POSITIVE PRINT OUT DIRECT ASSOCIATION. 

LB:*(T>Q)/LS 
CEA-*- ' INDIRECT' 

*L7 

L 5 : CEA* ' DIRECT 1 
L7 :PV*2xPVAL 
■+(PV£1 )/L8 
PV* 1 

LB:' KENDALL" S B EQUALS ',(4 vT) 

I I 
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'751 

.76] 

77] 

73j 

[79J 



'TEE P-VALOE FOR E C: NO ASSOCIATION EXISTS VERSUS ' 

| El: ' , (*CEA ) , 1 ASSOCIATION EXISTS IS: ' , (4 *PVAL ) 

'TEE P-VALUE FOR TEE TWO-SIDED TEST OF EIPOTEESIS JS:‘,(4*P7) 

t t 



CJ3 



: 



__ j1 

[29 

[3C 

[31 

:ii 

is 

,36 

.37 

.38 

[39 

40 

.41 

[42 

>3 

.44 

.45 

.46 

[47 

[48 

51 

52 
[53 

54 
! 55 . 
56 
.57 

58 

59 
[60 
:ei: 
*62] 



P°°4 

[67. 

:ii; 

N : 

■74; 
1 75 
.76^ 
. 77 j 
.78] 



7 KRUSKAL : WSM : CE.VOM : £ ; C : S : S : F : /M : S3 : OS ; S ; F ; S ; OF ; P : PVAL : R:SOFR ;SR; TSOR ; C5A 
a ISIS FUNCTION COMMUTES T’JE KbSKAL -WALLIS TEST STATISTIC 3 WHICH IS 
a i? MEASURE OF THE EQUALITI OF K INDEPENDENT SAMPLES. 

a SUBPROCRAMS CALLED BY THIS FUNCTION INCLUDE : TIES. TIESK, INDEXPLS, 
a FDISTN, INTERP AND THE VARIABLES PMATKW20 , P MAT KM 21 , PMATKW 33, 
a PMATKW3 4, PMATKM 41, PMATKM *2 , AND P MAT KM ^ 2 . 

QPP-*-5 

Si : 'ENTER THE NUMBER OF POPULATIONS TO BE COMPARED ( MUST BE CHEATER THAN I 
WO i . ' 

*((F<3)v((l| F)=0))/S1 
-((pF)>l)/Sl 

a INITIALIZE VECTORS E AND F AND VARIABLE C 

E-*F*-SOFR+o 0 
C*“ c 

arsis LOOP FACILITATES ENTERING THE SAMPLE VECTORS AND STORING THEM 
CH A* 'FIRST' 

£i:OC+l • 

1 ENTER IOUR ' , ( vCHA ) , ' SAMPLE . ' 

D+I 

(poD ) = C )/NEXT 
D*’ 1 c 0 

a COSC^ISV^rS SAMPLES AS THEY ARE ENTERED AND STORE THEM IN VECTOR E 
NEXT:E*S,D 

a RECORD THE LENGTHS OF THE SAMPLES AS THEY ARE ENTERED 

F*F,oD 
CEA*' NEXT' 

*(C<(F-1))/I1 

CHA-'LAST' 

*{C<K)/L 1 

a RECORD SIZE OF ALL SAMPLES WHEN COMBINED 

N+-/F 

a ORDER SAMPLE SIZES LARGEST TO SMALLEST 
OF+FZVP1 

a ORDER COMBINED SAMPLE VECTOR TO BE USED BI TIES FUNCTION 

Df-El&El 

a CALL INDEXPLS TO INCREMENT INDEXES WHEN TIES OCCUR WITHIN ONE SAMPLE 
AA+F INDEXPLS E 

a CALL TIES TO BREAK TIES BI MIDRANK METHOD 

BB+ 1 TIES D 
C* 0 

a THIS LOOP CALCULATES THE H STATISTIC 

L2:C+C+1 

a SUM OF RANKS FOR EACH SAMPLE IS CALCULATED 

SR++/BBZ (FCC] + (AA CC:] ) n 

a CALCULATE SUM OF kANKS SQUARED DIVIDED BI THE INDIVIDUAL SAMPLE SIZE 

SR+(SR*2)*FIC1 

a STORE EACH CALCULATION 

SOFR+SOFR ,SR 
*(C<K)/L2 

a , SUM ACROSS ALL SAMPLES 

TSOR++/SOFR 

H* (TSOR* ( 1 2 + ( Nx (N+l ) ) ) ) - ( 3 x (.N+l ) ) 
a RECALCULATE H WITH CORRECTION FOR TIES 
A+TIESK E 

NUM+(+/(A*3))-(+/A ) 

DEN0M+Nx((N*2)-ll 
H+H* ( 1 - (NUM*DEN0M ) ) 

a SISTEM OF LOGICAL STATEMENTS ENSURE PROPER PROB . IS ACCESSED 
0F[1]<2 ) /OUTPUT 
f F> 4 ) / F APPROX 
[K= 3 )/IF 

,(a/(0F= 3 3 3 3 ) )v (OFCl] >3 ) ) / FAPPROX 
,a/ {OF- 211 1 )) /OUTPUT 

2 1 1 ) )v (a/ (OF- 3 3 1 1))v(a/(0F= 3 3 2 1 ) ) 



j/Dli'l 

1 ) )v(a/(0F= 



(K= 4 )A(0F 
. (a/(0F= 3 
)/?41 

’(a=4)/?U3 

IF:*(0FCl]>4 )/F4PPS0X 
-►((a /( 0F= 2 1 1 ) ) v ( a/ (0F = 3 
-*-((0FC1]=4)a(0FC3]=1))/P31 

(a/ (0F= 321 

*((a/(0F= 322.. 

a CASS APPROPRIATE VARIABLE ACCESS CDF 

P23iP+PMATKW2QZ(.N-n)-, ;] 

■*PM 

P31 :P+PMATKW3lt (N-5 ) ; ; 3 
■*PM 

P 3 3 : P+PMATKW3 3 C (N- 5 ) ; ; ] 

■*PM 



1)))/ OUTPUT 

1) )v(A?fOF=-3 3 1 ) ) v ( a/ {OF- 3 

2) )v(0SCl]=4))/P34 



3 2))v(a/(0F= 3 3 3 ) ) )/P33 
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[7 9] P34:P^PMrJW34C(iV-6)5 ;] 

[80 +PM 

*81 P4l:P+PMAPPV4l[(tf-5); ;] 
r 82 +PM 

[83! PU 2 : P+PMATKW^ 2 C {N - 5 ) ; ; ] 

*84 ♦PM 

85 P43:P+PMAPPJ/43C(Ar-7); :] 

"86 a CALL INTERP TO CALCULATE P-VALUE BY INTERPOLATION 
*87] PM.-P7AL+3 INTERP P 
"88] +(PVAL= 1)/ OUTPUT 
* 89 , +L 5 

90 p CALCULATE P-VALUE USING TEE F DIST W/ONE LESS D.F.IN DENOM APPROX. 

91 FAPPROX:F+UN-K)*H)*{{K-l )x( {N-l )-H) ) 

92] P7AL+1-((U-1), ((#-£) -1)) FDISTN F) 

93 +L5 

94 OUTPUT :PVAL+' GREATER TEAN .25 1 

95 LS:'THE E STATISTIC EQUALS: ' , (4*5) 

96, 1 1 

97 'TEE P-VALUE FOR E 0: TEE POPULATION MEDIANS ARE EQUAL VERSUS 1 
“g 8 j 1 1 

C9 9J 1 PI : AT LEAST TWO POPULATION MEDIANS ARE NOT EQUAL IS: ' ,U?VAL) 



ty y_ 

Jioo] 
[ 101 ] r ^- 
.102 j ? 

Ci03] 
L104] 



t t 
*0 

k: 

t f 

+5 1 



'ERROR: YOU MUST ENTER A SINGLE INTEGER VALUE GREATER TEAN 2; TRY AGAI 



V MANNWHIT;N;M;PV2:A:3:C:G;MM:NN:RX;U:NM1;P;NU;PVAL:NM;NUMZ:NUMZ1;DENC1 :Z 
;Zl ; ALPHA : CDF : INb EX : IPX : Cl ; U ALPHA ; 33 ; CC ; Ul ; U2 ; PV ; MV1 ; LW2 : P7P ; LPPP ; aa ; GG 
;PVM;PV3:D;Q:R; DEN :DENC ;TC : PCI ; jVCM 
A PPis FbNCTiON usks TEk SlJM OF RANKS PROCEDURE TO CALCULATE TEE 
A MANN-WEITNEY U STATISTIC USED IN COMPUTING TEE P-VALUE FOR TEE TEST 
A OF LOCATION AND SCALE . P5P C.P. 505 M(I)-M(X) f P5P S5IPP PiV 
A LOCATION, IS ALSO COMPUTED . SUBPROGRAMS CALLED INCLUDE: TIES, TIES 2 
A INDEXPLS, VARMW, MANWP , INPUT, CONFMW , NORMCDF , AND NORMPTE . 

DIFF+-0 

i r 

l L0 IOC VJS5 PC COMPARE TEE MEDIANS OR VARIANCES OF TEE POPULATIONS? ' 

t ? 

B 2: f ENTER: 1 P0 CCMPA55 MPLJAtfS; 2 P0 C0MPA55 7A5JACCPS. 1 

AiH-0 

*((.AAx1)a(.AAx2))/E2 
+ UA = 1)/N1 

' THE TEST TO COMPARE VARIANCES REQUIRES THE TWO POPULATION MEANS' 

'OR MEDIANS TO BE EQUAL. IP THE I DIFFER BY A KNOWN AMOUNT , ' 

'THE DATA CAN BE ADJUSTED BEFORE APPLYING THE TEST. ' 
t r 

S3:' ENTER THE DIFFERENCE OF MEDIANS (M(X) - M(Y)) OR 900 TO QUIT).' 
DIFF + □ 

-*•( (oDIFF)>l )/E3 
■*lDIFF=900)/0 

'THE NULL HYPOTHESIS STATES - TEE POPULATION VARIANCES ARE EQUAL ; V(X) = 
7(7). ' 

■ i 
t 

i i 

N 2: 'ENTER: 
x 7(7). ' 



WHICH ALTERNATIVE DO YOU WISH TO TEST?' 

1 FOR Hi: 7(X) < 7(1); 2 FOR Hi: V(X) > 7(1); 



3 FOR Hi: V(X) 



D 'Rl 



'( (0*1 )a(L*2 )a(L*3 ) )/54 
+51 

Nl : ' TEE NULL HYPOTHESIS STATES - THE MEDIANS OF X AND Y ARE EQUAL; MU) 
= MU).' 

! t 

f WHICH ALTERNATIVE DO YOU WISH TO TEST? ' 

i r 

56 : ' ENTER : 1 FOR 51: M{X) < MU) ; 2 FOR Hi: MU) > M(Y); 

* MU). ' 

p+g 

+ ((L*i)a(Lx2)a(Lx3 >)/E u 



3 FOR Hi: MU) 



31 : ' ENTER 



DATA 



N+Q 

+ ( :oq_N)-0)/El 

1 ENT 

M+Q 



ENTER DATA VECTORS 
{MORE THAN ONE OBSERVATION IS REQUIRED).' 



Y DATA . ' 



A IF CALCULATIONS INVOLVE VARIANCES ADJUST X BY THE DIFFERENCE IN MEANS 
N+N-DIFF 

A CONCATENATE X AND Y SAMPLE VECTORS 

A+N ,M 

A DETERMINE SIZE OF X AND Y VECTORS AND ASSIGN TO NN AND MM 

NN+pN 
MM+pM 

A COMPUTE SIZE LIMIT OF LEFT TAIL OF NULL DISTRIBUTION 

NM+{NN*MM)*2 
NM1+INM 



90 



ORDER A AND ASSIGN TO B 



■4§i 



B+At&Al 
C*(NN , MM) INDEXPLS A 

A CALL TIES FUNCTION TO BREAK TIES USING MIDRANK METHOD 

G * 1 TIES B 

a IF FALSE CALCULATE TEST FOR VARIANCES 
*(AA=1 )/B 5 

o CALL VARMW TO GENERATE RANKS REQUIRED FOR VARIANCE TEST 

GG*VARMW (NN+MM ) 

a CALL TIES TO RECORD TIES IN THE DATA AND BREAK TIES IN GG 
G*GG TIES 3 

a , , CALCULATE SUM OF X RANKS 

3z:RX*+/(GL(CLl;(\NN)) )] ) 

a ,, , CONVERT TO MANNWRIT U STATISTIC 

U+RX- ( (NNx (NN+ l))+2) 

Ul+U 

a IF SIZE OF X TIMES SIZE OF 1 > 80; GO TO NORMAL APPROX 

*( (NMx 2 )>80 )/L 2 
NN1*\/NN,MM 
NN2*< /NN, MM 

aMANWP FUNCTION CALCULATES LEFT TAIL CUMULATIVE PROBS . OF U STATISTIC 
P+NN 2 MANWP NN1 

a LOGICAL STATEMENT ENSURES ONLI LEFT SIDE OF NULL DIST IS USED 

*(UZNM1 )/L 3 

a , CONVERT U STAT WREN GREATER TRAN LEFT TAIL VALUES 
U*(NNxMM)-U 

a IF U 1 IS A FRACTIONAL, INTERPOLATE P VALUE 
L3:*((1\U)=0)/NON 
U2*IU 
*(U2>0)/P1 
?V*l-((pr 02+12 ) + 2) 

~?3 

?1:?7«-1-((?C02)+P!;02 + 13 ) + 2) 

P3 :PVI*(P[U2+Ij +PCU2+2J ) + 2 
*CRECK 

NON:*(U>0)/GO 
PV* 1 
-?2 

GO:PV*l-?CU2 

P2:PVi+pr{u+i)} 

CHECK : + ( U1SNM1 } /Lu 
PV2+PV 
PV+PVI 
PVI*PV 2 
♦LU 

a COMPUTE TRE ■ NORMAL APPROXIMATION W /CORRECTION FACTOR 
L2:NUMZ*(U+0.5)-NM 
NUMZ1+U l-0.5)-NM 
DEN*((MMxNNx(MM+NN+l))+12)*0.S 
Z*NUMZ+DEN 

Z1*NUMZ1+DEN 
NUM* I ( U-NM ) 

DENC+l ( ( (NN+MM-1 )x (DEN*2 ))+ (NN+MM- 2 ) )-( ( (JVOM-0 . 5 )*2)+(NN+MM-2 ) ) )*0 . 5 
DENC1*( ( ( (NN+MM-l)x(DEN*2 ))+ (NN+MM- 2 ) )- ( ( (SOM+O . S )*2 )+ (NN+MM-2 ) ) )*0 . 5 
TC+(NUM- 0.5 )*DENC 
TCl*(NUM+0 . 5 )*DENC 1 
*(USNM)/ SECOND 

PVI*( (NORMCDF Z) + ( (NN+MM-2 ) TDISTN TC)) +2 

PV* ((1 - (NORMCDF Z1 ) ) + (l- ( (NN+MM-2 ) TDISTN TCl))) + 2 

*L 4 

SECOND :PVI*( (NORMCDF Z )+(!-( (NN+MM-2) TDISTN TC)))*2 



PV*(( 1 - (NORMCDF _Z 1 + ( (NN+MM-2 ) TDISTN TCI) ) + 2 ' 

.. . . J1 3/05 
P73+1 



£4:PV3-t-2x(L/(P7,P7I 
•(P73Sl)/iV5 



S5:P7M-*-(3,l)p(P7.T,P7.P73) 

jraF SOM OP TRE X RANKS IS: ',(vRX),’. TRE U STATISTIC EQUALS: ' , ($01) 



;i22, 

:i233 



a . . LOGICAL STATEMENT FOR VARIANCE OUTPUT 

' _TRE P-V ALUE_FOR _ RO : M(X) = M(X) VERSUS SI: MOO ' , (sLOC-TCCO; 1] ) , ' M(7) 



■W-4 = 2)/7;lS 
TSS P -VALUE FOR SC 
IS: ' , (4®P7MC0:13 ) 



t i 
-►18 



7/lS:P7M-(3 ,l)p(P7,P7I,P73) 

■ISS P-74S0P 70S SO: 7(.v) = 7(1) VERSUS SI: 7(X) ' , (vLOGICZD ; 1 ] ) , ' 7(7) 
PS: ' , (usP7MC£:i: ) 



LS: ’WOULD 70U LIKE A CONFIDENCE INTERVAL FOR TRE SHIFT IN LOCATION (M (I ) - 
M(X ))? (I/N).' 

BB*Z 

*(BB= <N' )/0 
110 : CC+INPUT 5 
ALPHA* (100 -CO + 200 

A ROUTE TO NORMAL APPROX. FOR CONF . INT . OF LARGER SAMPLE SIZES 
+((NMx2)>S0)/LS 

a COMPUTING CONFIDENCE INTERVALS BI EXACT P-VALUE 

CDF*P 

A INDEX POSITION OF VALUE IN CDF £ ALPHA 

INDEX* (+/( CDFSALPRA) ) 

*(INDEX>0 )/L6 

INDEX*1 




,166 

,167, 

El68J 



*L 6 

o COMPUTING CONFIDENCE INTERVALS USING NORMAL APPROX. W/C.F. 

L5 : UALPRA+ (.DENOMZx (.NORMPTH ALPHA ) )+NM-Q . 5 
fl ROUND U ALPHA DOWN AND INCREMENT BT ONE 

INDEX+lUALPHA+1 
L6:IPX*NN, INDEX 
CI+IPX CONFMW A 

> A ' , («CC), ' PERCENT CONFIDENCE INTERVAL FOR IRE SHIFT IN ' 

' LOCATION BETWEEN POPULATIONS X AND I IS:' 

( ' , UCTCl] ) , ' £ M(7) - M(X) £ ' ,(»CIC2] ),' )■ 
ERROR: THE SIZE OF TOUR SAMPLE IS LESS TRAN TWO ; TR7 AGAIN . ' 
ERROR: IOU RAVE ENTERED MORE TRAN ONE VALUE. TRI AGAIN.' 

ERROR: IOU RAVE NOT ENTERED A VALUE OF 1 OR 2; TRI AGAIN.' 

*32 

£4 : | ERROR : IOU RAVE NOT ENTERED A VALUE OF 1, 2, OR 3; TRI AGAIN.' 

*(AA=2)/N2 

+B 6 
V 



.23 



7 NPSLR: N ; SUMX ;SUM2:XBAR : 13 AR ; SUMX2 : SUMX2 : 3 ;A : WW ; XX ; 33 ; U : D ; ALPHA ;P:CC\ CDF 
; TALPkA ;NN ;CI : SLOPES : RR ; 5/?;77 xDENOM : INDEX ;FF:X;2 :Q; R: CHA ;PV 
A PROGRAM CONDUCTS NONPARAMETRIC LINEAR REGRESSION. THE LEAST SQUARES 
A REGRESSION LINE IS COMPUTED WITH HYPOTHESIS TESTING AND 

A CONFIDENCE INTERVAL AVAILABLE FOR THE SLOPE 3 . IF 3 DOES NOT LIE IN 
A THE C.I. AN ALTERNATE REGRESSION LINE IS PROPOSED . SUBPROGRAMS CALLED 
A ARE: SPMANP , KEN D ALP , NORMPTH , INPUT, AND CONFLR. 

UP P+5 

a livptfr mm 

R+INPUT 2 
Q+l *R 

X+1±(Q+1 )+R 
2+(Q+l)+R 



N+pX 

* SUMX++/X 
SUM2++/2 

XBAR+SUMX+N 

2BAR+SUM2+N 

SUMX 2++/(X*2) 



ASSIGN THE SIZE OF X ( AND 7) TO N 
COMPUTE THE SUM OF X' S AND 2'S 

COMPUTE THE MEAN OF X AND 1 

COMPUTE THE SUM OF THE X' S SQUARED 
COMPUTE THE SUM OF X TIMES 2 



SUMX2++/ (XxT) 

A COMPUTE • 5 • , THE SLOPE OF THE ESTIMATED LEAST SQUARES REGRESSION LINE 
B+{ (NxSUMXI)- {SUMXxSUMD )+ ( (NxSUMX2 )- {SUMX*2 ) ) 
a COMPUTE 'A 1 9 THE 2-INTERCEPT 

A+2BAR- (B*XBAR ) 

FF+ ' N ' 

| THE LEAST SQUARES ESTIMATED REGRESSION EQUATION IS: ' 

; 7= ',( *4),' + 

*00 700 70 ENTER SOME X VALUES TO GET THE PREDICTED 2"S? (2/N).' 

WW+U 

+(WW-'N' )/Ll 
L2\ 'ENTER X VALUES . ' 

XX+U 

a CALCULATE PREDICTED 2'S 

22+A+BxXX 

; THE PREDICTED 2 VALUES ARE: ',(*77) 

• WOULD 20U LIKE TO RUN SOME MORE X VALUES? ( 2/N 

ww+u 

-(WW=' 7' )/L2 

Li:' WOULD IOU LIKE TO TEST H2P0THESIS ON 3 , THE SLOPE OF THE EQUATION? (7/ 
;V ) . ' 

WW+U 

+(WW='N' )/L 3 

1 ENTER THE H2P0THESIZED SLOPE.' 

BB+Q 

a , % COMPUTE UI'S 

U+7-(BB*X) 

CHA + ’ > 1 

A CALL SPMANP TO COMPUTE RHO AND ASSOCIATED P-VALUES . 

D+X SPMANP U 
+(D111>0)/L11 
CHA+ ' < ' 

Lll:PV+2xDl2l 
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♦<P7S1)/£19 

PV+1 

L 19: 1 SPPAPMAtf' ! S 5 5QC/A5S: 



* , (4*5[l] ) 

= » , ($55), 1 VERSUS Hi: 3 



A*C3AL1 22 ) , ($55 ) , ' IS: 



63 

6U 

65 

66 
,57 

:i! 

72 

73 

74 

:p: 

bzlj 
”80 * 
C8 1 ! 
;s2 

V 83, 
84 
,35 
•36] 
*37] 
i.83 j 
8 3] 

M ZL Z 1 

_ 90 ^ 

C 9 1 1 

C§2. 

[9 3, 
;94j 
C 9 5 j 
96? 
;973 
[93] 
99] 
100 , 
, 101 , 
102 , 
,103 , 

,104 

105 

106 
.107, 
,108 1 
.109 
HO, 
. 111 . 
. 112 . 
.113, 
,114, 
[115 
.116, 
.117, 

118 
1 119, 
,120] 
.121] 
.122, 
.123 
.124, 
.125, 
.126, 
[127 
L 12 3 
.129, 
.130, 
11313 



'THE ? -VALUE FOR HO: 
f 1 f (4$5[2] ) 

» THE P -VALUE FOR THE TWO SIDED TEST OF HYPOTHESIS IS: *,(4$P7) 

1 1 

a I? USING TEE NEW REGRESSION EQUATION BASED ON MEDIANS, EXIT EERE . 
L3:*(.FF='I< )/0 

a COMPUTE CONFIDENCE INTERVALS ON 3 

> WOULD IOU LIKE A CONFIDENCE INTERVAL FOR TEE SLOPE? ( I/N ).' 

WW*S 

*(.WW='N' )/0 
LlO: CC*INPUT 5 

a CEANGE ENTERED VALUE TO ALPEA 

ALPEA* ( 1 0 0 -CC) * 2 0 0 

a ROUTE TO NORMAL APPROX. FOR CONF. INT . OF LARGER SAMPLE SIZES 
*(N>!2)/Lz 

a COMPUTING CONFIDENCE INTERVALS 31 EXACT P -VALUE 

P+XENDALP N 
CDF*?: 2:2 

a INDEX POSITION OF VALUE IN CDF <. ALPEA 

INDEX* ( + / ( CDFSALPEA ) ) 

*(INDEX>0)/LS 

INDEX*! 

*L 6 

a COMPUTING CONFIDENCE INTERVALS USING NORMAL APPROX. W/C.F. 

Lz :DENOM*( (.Nx(N-! )x< (2x/V) + 5 ) )*I8 )*0.5 
TALES A*DENCMx ( ( NORMPTE ALPEA ) ) 

*L 9 

Lo :T ALPS A*? :z-. INDEX! 

Ls :CI*X CONFLR I 
NN*!*CI 
SLOPES*!* Cl 
RR*i ( (NN-TALPEA )*2 ) 

*(RRs 0 )/L 20 
RR* 1 

L 20: SR* '(! + ( ( NN+TALPEA ) + 2 ) ) 

*(SEZ(oSLOPES))/L2! 

SR*o SLOPES 

L2!:< A ' (sCC),’ PERCENT CONFIDENCE INTERVAL FOR B, TEE SLOPE OP < 

1 TEE ESTIMATED REGRESSION LINE, IS:' 

i i 

J t ( ' , (.vSLOPEStRRl ). ' < B < ' , (vSLOPESZSR) ) , ' 

a IF 3 OUTSIDE TEE C.I. CALCULATE NEW EQUATION BASED ON MEDIANS 
a*( ( BiSLOPESLRRl )a (.BSSLOPESLSR2 ) )/0 
a ORDER X AND I 



). 



x*xztx: 

~ 7CA*] 



i*n . _ 

a CHECK TO SEE IF TEE SIZE OF SS IS EVEN OR ODD FOR FINDING MEDIANS 
*U2\NN) = Q)/S! 
a COMPUTE MEDIAN FOR ODD CASE 

B*SL0PESZUNN+!)*2)1 
*S2 

a , COMPUTE MEDIAN FOR EVEN CASE 

S! :B*(.SLOPES C (.NN* 2 )] +SLOPESZ ( (.NN+2 ) + 2 )] ) + 2 
a DO TEE SAME FOR TEE X AND I VECTORS 

S2 :*( (2\N)=0 )/S 3 
IBAR*lt((N+l}*2)l 
XBAR*XL ( (W+l ;+2 )J 
*OUT 

S3 : IBAR*(IZ(N*2 )] + JT ( (.N+2 ) + 2 53 )*2 
XBAR*(XZZN*2)2+XZ(.In+2)*2)2 )*2 
a COMPUTE NEW INTERCEPT 'A 1 

OUT : A*IBAR - ( BxXBAR ) 

'TEE LEAST SQUARES ESTIMATOR OF B LIES OUTSIDE TEE CONFIDENCE INTERVAL . ' 
J DISCARD TEE LEAST SQUARES EQUATION AND USE:' 

; I = + < ,(.vB),'X ' 

' THIS EQUATION IS BASED ON TEE MEDIANS OF TEE X AND I DATA AND ' 

'TSE MEDIAN OF TEE TWO -POINT SLOPES CALCULATED FOR TEE CONFIDENCE INTERV 
AL ON 3.' 



r !32 

1:2 3 1 a ALLOW USER TO DO SOME AN ALT SIS ON NEW EQUATION 
Cl 3 43 1 DO IOU WISE TO ENTER SOME X VALUES TO GET PREDICTED I"S FROM TEE NEW E 

QUATION? (.I/N).' 

FF* 0 

*(FF='I< )/L2 
FF* ' I ' 

+L1 

1 
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[29 

l30 

[31 

[32 

[33 

[34 

if 

ci2 

;4o 

.42 

;44 

:5i 

kl 

;5o’ 

il 

[55 

.56 

57 

[58 

:io 

:il: 

63 

:!3: 

[66 J 

[67] 

68 

leg; 

tZ?j 

l 72 j 

:?4; 

;76; 

77 

ill i 

,78, 

.79 

80 

:ai: 

. 32 . 

C 83 3 

[8 4 j 
,8 5 j 
: 36 1 



7 SIGN;A:C;3:D:PVAL;X:M0;N:CDF;ALPHA:CI:Y;AA;BB;CC;DD;PV;PVI;NNN;KP0S;0RD 
DiXXiOhDkiKALPaAiixZliQUAxmtiPVMiPVZiRiQ „ m „ 

a inis Function (jses the okDiNARY sign test to calculate the k 

a STATISTIC . P -VALUE. AND CONFIDENCE INTERVAL AS A TEST FOR MEDIANS . 
a THE LAST OPTION WILL DISPLAY A TABLE OF CONFIDENCE INTERVALS OF 
p ORDERED STATISTICS WITH CONFIDENCE COEFFICIENTS. 

a SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE : BINOM, NORMCDF , NORMPTH , 
p INPUT, AND QUANC. 

i i 

i DID YOU ENTER THIS PROGRAM FOR THE SOLE PURPOSE OF GENERATING CON 

FIDENCE INTERVALS FOR A SPECIFIED SAMPLE SIZE AND QUANTILE? ( Y/N 

ww*-g 

-+(WW= ' Y ' )/B4 

i THE NULL HYPOTHESIS STATES - THE POPULATION MEDIAN (M) IS EQUAL T 

0 THE HYPOTHESIZED MEDIAN (MO); HO: M = MO. 1 



WHICH ALTERNATIVE DO YOU WISH TO TEST?' 
1 FOR Hi: M < MO; 2 FOR HI : M > MO; 



3 FOR Hi: M * MO . 1 



2 FOR PAIRED-SAMPLE PROBLEM. 



COMPRESS D TO REMOVE ZEROS 
RECORD LENGTH OF A AND ASSIGN TO N 
KEEPING TRACK OF POSITIVE SIGNS 



i 

i i 

S3 : ' ENTER : 

C + 0 

-*< (Cxi )a(Cx2)a(Cx3 ) )/E3 
32: 'ENTER: 1 FOR SINGLE-SAMPLE PROBLEM; 

AA*-Z\ 

+((AA*1)a(AA=2))/E2 
-KiU = 2)/£9 

p INPUT DATA FOR SINGLE SAMPLE CASE 

X+INPUT 1 
NNN+oX 
MO+INPUT 3 
D+X-MO 
-£ 11 

p PAIRED SAMPLE CASE 

S3 : R+INPUT 2 
Q*-l +R 
X«-l + (Q + l 
Y*-(Q+l)+R 
DD+X-Y 
NNN+oDD 
MO+INPUT 4 
D+(X-Y)-MO 

Lll:A*-(D*Q )/D 

p 

N+pA 

*KPOS*-+/(A> 0) 

-►0^30 )/NORM 
PVAL+BINQM N 
+ (KPOS>0)/P1 
PVI+1 
+P2 

PI : PVI+l -PVAL ZKPOSl 
P 2 : PV+PVAL [ (KPOS+1 ) ] 

-►£6 

p IF N IS GREATER THAN 30 USE NORMAL APPROX W/ CONTINUITY CORRECTION 
NORM:Zl+( (KPOS-O.5 )- (0 . $xN ) )+ ( 0 . 5x (N*0 . 5 ) ) 

Z*-( (FPCS+O •5;-(0.5xjV);*(0*5x \N *0 . 5 ) ) 

PV+NORMCDF Z 
PVI+1- (NORMCDF Zl) 

p , IF P4IPFS S4MP£F TEST GO TO L 17 FCF OUTPUT STATEMENT 

L6:PV3+2*(l/(PV,PVI ) ) 

+ (P73<;i)/iV5 
P73+1 , % , 

tf5 :P7M«-(3 ,1 )p (P7.P7J.PV3 ) 

| COMPUTAT IONS ARE BASED ON A SAMPLE SIZE OF: ' , (vN) 

| THE TOTAL NUMBER OF POSITIVE SIGNS IS: ' , (vKPOS ) 

+(AA=2)/L17 

'THE P-VALUE FOR HO: M = ',(*M0),' VERSUS Hi: M ' , (vLOGIClC ; 1] ) , * ' ,(*M0) 
f 1 f J5: 1 , (4$P7M[C; 13 ) 

-►£18 

£17: 'THE P-VALUE FOR COMPARING THE MEDIAN OF THE POPULATION OF ' 

' DIFFERENCES TO THE HYPOTHESIZED MEDIAN , ' 



'HO: M(X-Y ) = * , (*M0 ) , 
: ' , (4«P7MCC;13 ) 

i I 



VERSUS Hi: M(X-I ) » , (zLOGICZC; 1] ) . 1 ',(*M0), 



L18: 'WOULD YOU LIKE A CONFIDENCE INTERVAL FOR THE MEDIAN? (Y/N).' 
BB *- Q 

-►(££ = 'Y' )/£16 
+QUANT 

p INPUT SIZE OF CONFIDENCE INTERVAL 

£16 : CC+INPUT 5 
ALPHA*- (100 -CC ) + 20 0 
+ (NNN'ZZQ)/NORMl 

p COMPUTING CONFIDENCE INTERVALS BY EXACT P-VALUE 

CDF+BINOM NNN 

p , x INDEX POSITION OF CDF FOR ALPHA + 2 

B*-+/ (CDFZALPHA ) 

-+(B>0) /SKIP 
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3*1 

*SKIP 

b COMPUTING CONFIDENCE INTERVALS BY NORMAL APPROX. 

NORM 1 : XALPHA* ((0.5* (NNN*0 . 5 ) ) x (.NORMPTS ALPS A ) )+ ( 0 . 5 xJVAW ) - 0 . 5 
b ROUND RALPH A DOWN TO NEAREST INTEGER AND INCREMENT 31 ONE 
B*IKALPHA+1 

b , IF SINCLE SAMPLE CASE GO TO LI 

SKIP:*(AA=1)/L7 

fl CALCULATE AND PRINT OUT CONF. INT . FOR PAIRED SAMPLE CASE 

LS :0RDD*DDL£xDDl 
YY*oORDD 

CI*0RDD:31,0RDDt(,YY-(3-l'))l „„ , 

'A ',(5CC) ' PERCENT CONFIDENCE INTERVAL FOR TEE MEDIAN OF THE ' 



( ' , ), 1 S MEDIAN (X-Y ) 5 ' , ($CT[2] ) , ' )' 



"I • V / « a uilvL'irX UL/iTr 1 L/Dif Uu 

'POPULATION OF DIFFERENCES IS: 1 

r f 
I 

I I 

+QUANT 

q CALCULATE AND PRINT OUT CONF. INT. FOR ONE SAMPLE CASE 

LI :ORDX+XZAX3 
77+0 ORDX 



HO, 

.111, 

.112, 

113J 



CI+ORDXL31 ,ORDXt(YY-(B-l))l 

!A-L>(*CC), ' PERCENT CONFIDENCE INTERVAL FOR TEE MEDIAN OF TEE POPULATIO 
( 1 , (*CT[1] ), ' £ MEDIAN £ 1 ,(sCJ[2]),' )' 



N IS: 1 

i i 
i 

i i 

ANT: 'WOULD 10U LIKE CONFIDENCE INTERVALS FOR A SPECIFIED QUANTILE? (I /N 



QUAN 
WW+T 

+ CWW= '7' )/B 1 
-0 

3^: ' ENTER DESIRED SAMPLE SIZE ( SINGLE INTEGER VALUE) 
NNN^L 



+ ^pHNNJ>l)(E± 



, . 1\NNN)*0)/E 4 

31: 1 Ex TER DESIRED QUANTILE; FOR EXAMPLE : ENTER 20, FOR TEE 20 TH QUANTILE. 
QUA+H 

+( (QUAS3)v(gUA>100) )/El 
NNN QUANC QUA 

t t 

* ***** this table gives confidence coefficients for various intervals w 

ITE 1 

1 0RDER STATISTICS AS TEE END POINTS FOR TEE ' $ (*QUA) ,'TE QUANTILE. 1 

t i 
+0 

El: 'ERROR: TEE QUANTILE VALUE MUST LIE BETWEEN 0 AND 100; TRY AGAIN. ' 

-Si 

El: 'ERROR: YOU EAVE NOT ENTERED A VALUE OF 1 OR 2; TRY AGAIN.' 

i t 

-S2 

E3:' ERROR: YOU EAVE NOT ENTERED A VALUE OF 1 , 2 , OR 3 ; TRY AGAIN. ' 

i t 

-S3 

34: 'ERROR: YOU EAVE NOT ENTERED A SINGLE , INTEGER VALUE ; TRY AGAIN. ' 

i t 

-S4 

7 



7 SPEARMAN *X • T* A • Q *R • CEA *PV 

q TEIS* FbNCTloh COMPOTES TEE SPEARMAN R STATISTIC WEICE MEASURES 

q TEE DEGREE OF CORRESPONDENCE BETWEEN RANKINGS OF TWO SAMPLES. THE P- 
q VALUE IS GIVEN FOR TESTING ONE AND TWO-SIDED EYPOTEESIS OF ASSOC. 
q SUBPROGRAMS CALLED BY TEIS FUNCTION INCLUDE: TIES , TIESK , SPEARP , 
q INPUT , SPAPROX, INTERP , AND TEE VARIABLE PMATSP. 

R+INPUt 2 
0-1 fS 

X-l+(0+l)+£ 

r-($+i)+s 
o CALL SPMANP 
A+X SPMANP I 
-<~Ci:>0 )/Ll 
CEA+ ' INDIRECT ' 

+L2 

L 1 :CEA+' DIRECT' 

L2:PV+2*AZ2l 
+(PVZ1)/L3 
PV + 1 

L3 : 1 SPEARMAN ' 'S 

i 



TO CALCULATE THE STATISTIC AND ASSOCIATED P -VALUES 



R EQUALS ' , (4»ACl] ) 



THE P-VALUE FOR HO: 

Hi: 

i 



NO ASSOCIATION EXISTS VERSUS ' 

' , (»OT), ' ASSOCIATION EXISTS IS : ' , (4®A [2] ) 



THE P-VALUE FOR THE TWO-SIDED TEST OF HYPOTHESIS IS : ' , (4«P7) 
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$ 

"u 1 

it] 

Bj 

£ 9 ] 

[ 10 , 

Cll 

[ 12 , 

[13] 

[14] 

[is, 

£ 16 

[ 17 , 

Lis] 

[19] 

, 20 ] 

[ 21 ] 

[ 22 , 

£ 23 3 

[24, 

in ] 

tih 

:Io; 

.31, 

.32, 

33, 

34, 
35] 
36, 

[37 

til 

;22 

■21 

!2S 

:47 

48 

,49 

ji 

Js 

:s7. 

i ! 3 

Ei 5 i 

[62: 

,63 

.64; 

.65 

.66. 

.67, 

[ 68 , 

[69, 

[70, 

[71, 

[72 

[76 

77 

C79 

:|f 

.83 

,84 

85 

*86, 

.87 



V WILCOX ; A ;B } D :E ; F:PV 2 ; Z1 ; Z;DEN ; NUMZ ; NUMZ1 \PVAL ; X : MO ;N: TPLUS : CPF : MLPPA ; A 
PPFA:F:CI: I ;AA ;3B ;CC ;NW ;DD ;PV ;P0S ; IPOS ;M ;PVI ;IP0S1 ; NNN ;C ;PVM;p73 ;R;Q;T 
RAP ; : D&NT1 : TC : TCI 

p raJS PFCCFAW CSPS rap WILCOXON SIGNED RANK TEST TO CALCULATE THE T 
p STATISTIC , P-7APC/P, AM) CONFJPF/VCP J/vraFFAL AS A TEST FOR MEDIANS . 
p SP5PFCCHAWS CALLED BY TRIS FUNCTION INCLUDE : TIPS, WPP, NORMCDF , 
p^SPWPra, CONFW , AtfP INPUT. 

1 TEE NULL HYPOTHESIS STATES - OTP POPULATION MEDIAN (M) IS EQUAL T 

HYPOTHESIZED MEDIAN (MO); HO: M = WO . 1 



0 rap 

1 t 



> 1 

S3:' ENTER: 

c (?c 



WHICH ALTERNATIVE DO 100 WISH TO TEST?' 
1 FOR HI: M < MO i 2 FOR HI: M > MO; 



3 FOR Hi: M x MO . ' 



2 FOR PAIRED -SAMPLE PROBLEM . ' 



PAIRED SAMPLE CASE 



COMPRESS D TO REMOVE ZEROS 
RECORD LENGTH OF A AND ASSICN TO N 
KEEPING TRACK OF POSITIVE SICNS 



(.(Cx l ) a (Cxi )a(C*3 ))/E 3 
52:' ENTER: 1 505 SINGLE-SAMPLE PROBLEM ; 

AA*0 

*((AAxi)\(AAx2))/E2 
+(AAx2)/L$ 

a INPOT DATA FOR SINGLE SAMPLE CASE 

X*INPOT 1 
NNN*oX 
MO*INPOT 3 
D*X-MO 
*L 11 

£9 :R*INPOT 2 
Q+l+B 

X+l + (<3+l HP 
r+(«+i)*5 
DD*X-Y 
NNN*oDD 
MO+INPOT u 
D*(X-Y)-MO 

Lll:A*(DxO )/D 
A 

tf+p;i ' 

POS*(A>0) 

a , 5455 THE ABSOLOTE 7AL0E OF A: ASSICN TO 3 AND ORDER 3 

5+ 14 
5+5 CAB] 

A REORDER POSITIVE SIGNS TO COINCIDE WITH PROPER POSITIONS IN B 

POS*POSZl(\A)l 

a CALL FONCTION TO BREAK TIES 

5+1 TIES B 

a CALCOLATE TPOS BY ADDING ACROSS ALL POSITIVE VALOES OF E 

TPOS++/ (POS*E) 

TPOSl*TPOS 

A GIVES SIZE OF LEFT TAIL OF PROBABILITY DISTRIBOTION 

MJ +( l ((+/ i 5 )* 2 ))+l 

A CO TO STATEMENTS BASED ON LENGTH OF VECTOR E 

+(5>16)/£3 

a GENERATE NOLL DISTRIBOTION FOR TPLOS 

F*WILP N 

a IF TPOS FALLS IN LEFT HALF OF PROB DIST CALCOLATE PVALOE AS NORMAL 
* (TPOSS, (NM - 1 ) ) /TREC 

A OTHERWISE OSE THE NEGATIVE T STATISTIC 

TPOS* (+ / :N ) -TPOS 

A IP TPOS IS FRACTIONAL OSE BOTH THE INTEGER ABOVE AND BELOW AS TPLOS 

TREC :+( (1 JTPOS ) x 0 ) /NON 

TPLOS+LTPOS 

Pi :PV*l-( (FITPLOS1 +FZTPLOS+11 )+2) 

pvi*(fztplos+h+fztplos+21 )+2 

+CHECK 

NON :* (TP0S>0 ) /GO 
PV* 1 
+ P 2 

co:PV+i-pr(TPOS)i 
P2:PVI*FZ(TPOS+l )] 

CHECK:* (TPOSl<.(NM-l ) )/L& 

PV2*PV 
PV+PVI 
PVI*PV 2 
"Co 

3 COMPOTE NORMAL APPROX. 

L2 :TRAP*(Nx(N+ 1 ) )44 
NOMZ*(TPOS+Q . 5 )-TRAP 
NOMZl* (TPOS-Q . S )-TRAP 
DEN*((Nx(N+ 1 )x ( (2xtf)+l ) )+24)*0. 5 
Z*NOMZ*DEN 
Z1*N0MZ1*DEN 

a COMPOTE STODENT T APPROXIMATION WITH CONTINOITY CORRECTION FACTOR 



W/CONTINOITY CORRECTION FACTOR 



NOM * 1 (TPOS- TRAP) 

DENT* l ( (Nx (DEN* 2 ) )*(N-1 ) )r ( ( (NOM-O . S)*2)+(iV-l ) ) )*0 . 5 
DENT 1*(( (Nx (DEN* 2 ) )*(N- 1 ) ) - ( ( (NOM+O . 5 )*2 ) + (N- 1 ) ) )*0 . 5 



TC*(NUM-0 . 5 )*DENT 
TC1*(NOM+0 . 5 )*DENT 1 
A COMPOTE AVERAGE OF TC AND ZC 
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+{TPQSS((+/\N)*2 ) )/ SECOND 

PV+( Cl ‘(NORMCDF Zl ) )+ (1- ( (tf-1 ) TDISTN TC l)))+2 
PVI+{(NORMCDF Z ) + ( (iV- 1 ) TDISTN TC)) *2 
+L 6 

SECOND :PV+ C Cl- CNORMCDF Z1))+CCN-1) TDISTN TC1))*2 
PVI+CCNORXCDF Z)+Cl-CCN-l) TDISTN TC)))* 2 
LZ>:PV3+2 *C'l/CPV t PVI) ) 

+CP73S1 )/W5 
P73 + 1 



//5 :P7M«-(3 , 1 )o CP7I ,PV % PV3 ) 

1 COMPUTATIONS ARE 3ASED ON 



145 

146 
147] 
148 

.1451 

150] 

,151j 

,152J 

,153J 

154] 



SAMPLE SIZE OF: ' f C*N) 

, (*rposi) 



[10 5] 
[106] 
'107* 

Eios] 



Cm, 
tll2 j 
,113] 

■m- 



'TEE TOTAL SOM OF POSITIVE RANKS IS 

t t 

a IF PAIRED SAMPLE TEST CO TO L17 FOR OOTPOT STATEMENT 

-(AA=2)/L17 

' THE P-VALOE FOR R0: M = ' , VERSUS Hi: M < , (.vLOGICLC ;1] ) , ' <sM0 
), ; IS:' , (4*pracc;l] ) 

—519 

LIT. 'THE ? -VALUE FOR COMPARING THE MEDIAN OF TEE POPULATION OF ' 

; DIFFERENCES TO TEE HYPOTHESIZED MEDIAN,' 

•EC: MCA-2 ) = ' .(sMO),' 7ERSOS 31 : Mf.X-1) ' , izLOGICZC; 1Z ) , ' '.CsMO),', * 
S: J , (U9?raCC;l] ) 

LIS: 'WOULD IOU LIKE A CONFIDENCE INTERVAL FOR THE MEDIAN? (7 /N).< 

B3-G 

*(3B='1< )/L 16 
♦0 

L--3:CC*INPUT 5 
^£P5a*UC0-C^) + 2C0 

a ROUTE TO NORMAL APPROX. FOR CONF INT OF LARGE SAMPLE SIZE 

-(NNN>lo )/L- 

a COMPUTING CONFIDENCE INTERVALS 31 EXACT P-VALUE 

CDF+WILP NNN 

a INDEX POSITION OF CDF FOR ALPEA * 2 

T ALPHA * ( + / < CDFSALPEA ) ) 

*(TALP3A> 0 ) /JUMP 
T ALPHA* 1 
-JUMP 

a COMPUTING CONFIDENCE INTERVALS 31 NORMAL APPROX. W/C.F. 

m:DENOMZ+<.<.2xNNNx(NNN-H)*(.(.2xNNN)+l ) ) + 3 )*0. 5 
T ALPS A* (( (DENOMZx (NORMPTE ALPEA ) ) + ( NNNx (.NNN-t-1 ) )-2 ) + 4 ) 

« ROUND TALPEA DOWN TO INTEGER VALUE AND INCREMENT HI ONE 

TALPBA*\.T ALPHA + 1 

a IF ONE SAMPLE CASE CO TO LI ■ 

JUMP:*(.AA-1 )/L 1 

a CALCULATE AND PRINT OUT CONF. INT. FOR PAIRED SAMPLE CASE 

L5:CI*T ALPHA CONFW DD 

'A >,(9CC).' PERCENT CONFIDENCE INTERVAL FOR TEE MEDIAN OF TEE ' 
'POPULATION OF DIFFERENCES IS:' 

I I 

J t ( ' , (sCI[l] ), ' £ MEDIANCX-I) ' ,(*cr[2]),» )' 

♦0 

o CALCULATE AND PRINT OUT CONF . INT . FOR ONE SAMPLE CASE 

L7-.CI+T ALPHA CONFW X 

' A 1 ,( QCC),' PERCENT CONFIDENCE INTERVAL FOR THE MEDIAN OF THE POPULATIO 
N IS:' 



( »,(*CI[l]) i ' £ 
ERROR: IOU HAVE NOT ENTERED A 
ERROR: IOU HAVE NOT ENTERED A 



MEDIAN £ 1 ,(*CI[ 2]), 1 )» 

VALUE OF 1 OR 2; TRI AGAIN . 1 
VALUE OF 1 , 2 , OR 3 ; TRI AGAIN.' 
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APPENDIX F 



LISTINGS OF SUBPROGRAMS BASIC TO BOTH WORKSPACES 



7 CBN+BINOM N-.N-.PiXiK-.CDF 

a TBIS FUNCTION iS A SUBPROGRAM OF TBE SIGN TEST (SIGN) . IT CALCULATES 
A TBE CDF OF TBE BINOMAL WBEN PROBABILITI = .5. N=S AMPLE SIZE. 

P* 0.5 
K* 0 , iN 

CBN*+\ (K’N)x (P*K)x ( (l-?)*N-K) 

7 



LSj 



7 U+L CBIN R 



A TBIS FUNCTION IS A SUBPROGRAM OF CONFIDENCE INTERVAL GENERATOR FOR 
A TBE QTB QUANTILE. IT RETURNS TBE VALUE OF TBE BINOMIAL CDF AT R, 
a WITH N ,P-L WBERE N-SAMPLE SIZE AND P=PROBABILITI . 



£M+\(( l + il+£ClJ )!£Cl3 )*C£l2]» 

1J 

7 



L-rtl+£ClJ )x ((1-LZ22 )*LZl2-~l + il+LZll))LR+ 



10 ] 

, 11 , 

12 , 

13 

;i 4 ! 

:iS: 

III: 

19] 

. 20 , 

. 21 ] 

, 22 , 

.23, 

.24, 

[25] 



V CON+XX CONFLR 71 ;3B ;SS ;AA:A-,XR:YR;3 :C;S 
A , TBIS FUNCTION IS A SUBPROGRAM OF NONPARAMETRIC LINEAR REGRESSION 
A (NPLR ) . IT CALCULATES THE TWO -POINT SLOPE FOR EACH PAIR OF POINTS 
a (XI .11) AND iXJ.YJ.) ALL KJ AND XI=XJ . ALL SLOPES ARE ORDERED AND 
a USED TO FIND TBE CONFIDENCE* INTERVAL FOR 8 t THE SLOPE OF ESTIMATED 
A EQUATION . XX = X DATA SAMPLE AND 11= I DATA SAMPLE. 

A RECORD TBE SIZE OF XX AND INITIALIZE VARIABLES 
BB+pXX 

SS + pO * 

AA+ 0 

a TBIS LOOP COMPRESSES XX AND YY DOWN TO WBERE TBE XX < ALL OTBER XX l S 
L 2 :AA+AA + 1 
A+ZXXZAA1<XX) 

XR+A/XX 
YR+A/Yl 
B+oXR 
"► (S = 0 )/L3 
C+ 0 

a TBIS LOOP CALCULATES TBE SLOPE OF EACH PAIR OF PAIRED DATA. 

LlzC+C+1 

S+illZAAl-XRZCl )+{XXZAAl-XRZCl ) 

SS+SS.S 
+(C<B)/L1 x 
L3:+(AA<BB)/L2 
CON+(pSS) t ( SSZ&SS1 ) 

7 



7 CONFM+AA CONFMW BB:A ;B;C:D:E;F:G;B 
A TBIS FUNCTION IS A SUBPROGRAM OF TEE MANN-WBITNEY TEST (MANW ) « 

A IT COMPUTES CONFIDENCE INTERVAL ENPOINTS FOR THETA , TBE SHIFT IN 
A LOCATION BETWEEN X AND Y. AA- INDEX POSITION OF C.I. ENPOINT. 

A BB=COMBINED DATA SAMPLES . 

a ASSIGN SIZE OF X VECTOR TO A; INDEX POSITION FOR CONF INT TO B. 

A+AAZll 



3+AAl 2j 

*C+A+3B 

D+A+BB 

OCCACJ 

D+DlbD) 



ASSIGN X VECTOR TO C; Y VECTOR TO D 

REORDER X AND Y VECTOR VALUES TO ASCENDING ORDER 

INITIALIZE VECTOR B AND VARIABLE F 



H+ p 0 
F + 1 

& INNER AND OUTER LOOPS CALCULATE ALL POSSIBLE DIFFERENCES : EVERY Y 
a ELEMENT MINUS EVERY X ELEMENT. B VECTOR STORES THESE DIFFERENCES. 
L2:G+1 

Ll:E+DZFl-CtGl 

B+H,E 
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i r 223 C+c+i 
23 J -(C£(pC))/Ll 
2UJ F+F+ 1 
23] +(FZ(pD))/L 2 

26 J a REORDER B VECTOR VALUES TO ASCENDING ORDER 

27] B+Rr&Bl 

23] fl IWC2X APPROPRIATE CONF INT VALUES FROM B 

29] C 0 NFM+BZB 1 ,B[ ( (ptf)-(S-l ))] 

V 



!ii] 

[ill 

[jU 

?| 

13 

L 20j 

E 21 5 

L 22 J 
L23J 

L2U] 



V CONF+A CONFW 3 ; C ; D ; E :G ; H :N ;F ; A :B 

o TEJ5 FUNCTION IS A SUBPROGRAM OF THE WILCOXON SIGNED RANK TEST 

0 (WISIG). IT PROVIDES CONFIDENCE INTERVAL END POINTS BASED ON THE 

0 AVERAGES OF ALL PAIRS OF DATA SUCH THAT ALL XI S XK. A= INDEX 

o POSITION OF C. I. END POINT AND B= INPUT DAT AD 

C+BZ&31 

p START H ACCUMULATION VECTOR OFF WITH ORGINAL VECTOR VALUES 

3+0 0 
B+B , C 
D+0 

a OUTSIDE LOOP INCREMENTS D AND RESETS E TO D 
L 2 • D+D + 1 
E+D 

P INSIDE LOOP GENERATES NEXT SET OF AVERAGES AND CONCATENATES TO ORIGINAL 
VECTOR 
LHE+B+1 
F+(CZD1+CZE1 ) + 2 
3+3,? 

* CONTINUE INNER LOOP UNTIL E EQUALS THE SIZE OF C 

+(E<(pC))/Ll 

4 CONTINUE OUTER LOOP UNTIL D EQUALS THE SIZE OF C LESS ONE 

+iD<((?C)-l))/L2 

p ORDER FINAL ACCUMULATED VECTOR H 

B+BUB1 

p INDEX CQNF INT VALUES OUT OF H 

conf+bza: , sc (( pff )- a - i ))] 

V 



[l] 

*3 j 
[ 4 * 

w 

■A ■ 

.10 

111, 

[12 

-13* 

14 

15 
.16, 

17 

18J 

,20 

.21 

.22 

.23 

,24 
25 
,26 ^ 
27 
.23 
.29 
.30 

C 3 — 
22 
* 2 3 1 

C 3 5 

£36 

37 

I 38 

39 

[40 

,41 

.42 

43 

.44 

[45 

.46 

[47 

L43 



7 P+DF FDISTN X;A;M;N :RM:RN ;LM;LN :SM: SN :M2 ;N2 
p THIS FUNCTION IS SUBPROGRAM OF THE KRUSKAL -WALLIS TEST ( KRWL ) AND 
p THE STUDENT T DIST. (TDISTN ) . IT CALCULATES APPROX . CUMULATIVE PROBS . 
p OF X USING TBE F DISTRIBUTION W/ DF= (M ,N ) DECREES OF FREEDOM . 

R TAKES VECTOR ARG. 

A+ ,M*X*N+XxM+l+DF t N+l+DF 
+L*\ (M>2 )vW> 2 

R TREAT TBE 1ST . FOUR SPECIAL CASES . 

+*'L' . 1 0 35 M.N 
Lll:P+{ lo4*o.5)+o+2 
*►0 

LI 2 :P«-4*0 . 5 

L21:P«-1-(1-A)*0.5 

*►0 

L22 :P+A 
+ 0 

p BEGIN THE GENERAL CASES. INITIALIZE TBE QUANTITIES . 

LlRM+il \M*2) ,0/RN+l \N*2 
SM+(RM= 0 )xRMSLM+LM2+0 .5*M-2 

- — - - . 5 xN-2 

FIRST. 



SN+(RN= 0 )*RNSLN+IN2+Q . 

P TREAT SN.CASE N-EVEN 
+N0DDx\l=2)N 
+MEVEN*\Q=N2 

SN+ l++/x\(l-4 )°.xl+( 1+0. 5xM)*\N2 
+MEVEN 

NODD:+MEVENx\l=N 

SN+(l-A)*0.5 

P TREAT THE PORTION OF SN THAT DOESN'T CHANGE W/ M ODD OR EVEN. 

** ( 2 + 1 + ULC ) x 1 0 = LN 

SN+SN*l++/*\ (l-.4)o t xl+2xM2 + l+2xiLiV 
p TREAT THE M-EVEN SUBCASE 
+{ (M-2 ) ♦ 1=2 \ M )/MEVEN , 3+ltCLC 
SN~SNxx/l^*2xiM2 
+MEVEN 

p TREAT TBE M-ODD SUBCASE 
SN+SN* 2+ol 
+ENDx\l=M _ 

SN+SNxx / 1 + + l+2x X 0.5+M2 
MEVEN:+M0DDx\l=2\M 
+ENDx\2=M 

SM«*l ++/x\A° , xi - + 2x \M2 
+END 

p NEXT TREAT TBE SPECIAL CASES FOR ODD-M. 

M0DD:SM+( QxM-1 )+( ( A*0 . 5 ) + o + 2)x3 =M 
+END*\(l=M)v3=M 

SM+(A*0 . 5 )x ( + 0+2 )xl + + /x\4«» . xl-+l+2xiLAf 
END :P+ (SN*A*M+2 ) + (4 xRNxRN-RM)+( (RMxRNx loA*Q<.S ) + o + 8 ) -2xRNxSMx (1 -4 )*0. 5 
7 
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7 XX*N INDEXPLS 3 :AA :3B ; CC:DD ;X1 :NN:NC ; N : 3 : C :D : S ;F 
A THIS FUNCTION lb A SUE PRO CRAM OF WRUSKAL -WALLIS (KRWL ) , AND MANN- 

p WHITNEY ( MANW ). IT COMPUTES A MATRIX OP RANKS OF THE ORIGINAL DATA 
A WHERE THE FIRST ROW OF THE MATRIX IS THE FIRST SAMPLE. IF TIES OCCUR 
a IN THE RANKINC WITHIN A SAMPLE, SUCCEEDING TIES ARE INCREMENTED BY 
A ONE . N- SIZES OF THE SAMPLES BEING PASSED IN B. B- ALL DATA SAMPLES 
A COMBINED. 



AA*pN 

*3B+r/N 



RECORD SIZE OF N TO DETERMINE NUMBER OF SAMPLES IN B 
LOCATE LARGEST SAMPLE SIZE 
ORDER B SMALLEST TO LARGEST 



DD+BL&B1 

a SET UP MATRIX OF SIZE REQUIRED TO STORE SAMPLE RANKING. 

XX*(AA ,BB )pO 

a CONCATENATE 0 AND N 

NN*0 ,N 

a FIND CUMULATIVE SUMS OF SAMPLE SIZES. 

NC**\NN 
CC* 0 

a THIS LOOP INDEXES OUT ORGINAL SAMPLE VALUES FOR FURTHER CALCULATIONS. 

LI * CC*CC+1 

A ' LOCATE CC'TH ORIGINAL SAMPLE VALUES IN S 

Xl*BL (NC(CC) +\NtCC] )J 

a RECORD POSITIONS OF ELEMENTS OF X VECTOR IN B AND ASSIGN TO C 

C*DD\X 1 
D* 1 

a THIS LOOP DETECTS TIED INDEXED POSITIONS AND INCREMENTS THE INDEXING OF 
EACH SUCCESSIVE TIED POSITION BY ONE 

a TEST DTH ELEMENT OF C AGAINST REST OF C FOR TIES 

L2:E*(CID1=D*C) 

A SET F EQUAL TO THE APPROPRIATE SIZE (.DTH SIZE ) ZERO VECTOR 

F*DoO 

a CONCATENATE F AND VECTOR OF O' S AND 1 'S (1 ' S APPEAR WHEN TIE OCCURED ) 
F+F E 

A* ’ ADD RESULTANT F VECTOR TO C 

C*C+F 
D*D + 1 

a % CONTINUE FOR ENTIRE C VECTOR 

*(Di(pC))/L2 
XX[CC:]*C. ( BB-NLCC ] )pO 
*(CC<AA )/Ll 
7 



F13 

i 

‘ 8 

;gj 

B io; 

.13, 

■is 

¥ 

,20 * 
.21 
,22< 
.23 
.241 
25 

[26* 
.27 
[28 
. 29 
,30 
.31 
32 
.33 

34 

35 
.36, 
.37, 
.38 
.39 
[40 

L 4 1 

l 42 



7 IN+INPUT A; B;X;Y 

A THIS FUNCTION DOES MOST OF THE INPUT PROMPTING AND ERROR CHECKING . 
ft A=THE TYPE OF PROMPT DESIRED . 

A IT IS A SUBPROGRAM OF SIGN , WISIG, MANN, KEN . SPMAN , AND NPLR. 
+(L1,L2,L3 ,L4.L5 ) LAI 

LI : * ENTER THE DATA ( MORE THAN TWO OBSERVATIONS ARE REQUIRED). ' 

IN+Q 






*o 



:2)/El 



L2:< ENTER X DATA (.MORE THAN TWO OBSERVATIONS ARE REQUIRED).' 

*((ppX)=0)/£l 

-*-((ox)=2)m 

1 ENTER Y DATA (NUMBER OF Y ENTRIES MUST EQUAL NUMBER OF X ENTRIES).' 

)/£2 

IN+(pX\,X,Y 

L3:' ENTER THE HYPOTHESIZED MEDIAN.' 

IN+U 

+ UpIN)>l)/E3 

L*:' ENTER THE HYPOTHESIZED MEDIAN FOR THE DIFFERENCES OF THE PAIRED DATA , 
IN+U 

+ UpIN)>l)/E3 
— 0 

Lz : 1 ENTER THE DESIRED CONFIDENCE COEFFICIENT ; ! 

i i 

* FOR EXAMPLE: ENTER 9S, FOR A 95 PERCENT CONFIDENCE INTERVAL.' 

IN* 0 



(JW>100 ) )/E 4 



♦0 



) /E5 



El: 'ERROR: THE SIZE OF YOUR SAMPLE IS LESS THAN THREE ; TRY AGAIN.' 

*(A=1)/L1 

E2:< ERROR: SAMPLE SIZES ARE NOT EQUAL ; WANT TO TRY AGAIN? (Y/N). ' 
B* 0 

*(B = ' Y< )/L2 

'ENTER RIGHT ARROW * TO QUIT.' 



1 00 



;uuj 

, 4 5j 

C 46 

C47 

.43, 

,49* 

[5C 

L51] 

i§2i 

f 53 J 



5-*C 

F3 : J ERROR : XFF HIPOTHESIZSD MEDIAN MOST BE A SINGLE VALUE ; TRI AGAIN.' 

+(A=3)/L3 

*►£4 

E'-j : j ERROR : THIS VALUE MUST LIE BETWEEN 0 AND 100; TRI AGAIN.' 

*L5 

E5 : j ERROR : T3IS VALUE MUST BE AN INTEGER ; Ti?J AGAIN. ' 

*L 5 
7 



7 INTER* A INTERP B ; C : D \E ; F :G : FF ;GG ; PL 
a TfflS FUNCTION IS A SUBPROGRAM OF THE KRUSKAL-W ALLIS (KRVL) , 

a tfSMMLL’S S (.KEN), AND SPEARMAN'S R (SPMAN AND SPMAN1 ) . f/HE’W 
a P4SSP0 rap rasr STATISTIC as tee left argument , it calculates the 

a ASSOCIATED P-VALUE 31 INTERPOLATION OF THE TABLE VALUES OF STATS . 
a AND CDF'S VETCH ARE PASSED AS TEE RIGHT ARGUMENT. 

a SEPARATE TEE CDF TABLE AND STATS . INTO SINGLE VARIABLES. 

GG*Si2 - j 

a FIND VEERS A FIRST EXCEEDS OR EQUALS ONE OF TEE TABLE VALUES. 
C*<\ (AZFF ) 

a IF TEE STAT. DOES NOT EQUAL ANI TABLE VALUES SET INTER - 1. 

* ((+ / C )-Q) / L2 

a INDEX LOCATION OF FIRST OCCURENCE OF MATCE 

a"'"''* IF INDEXED POSITION EQUALS ONE INDEX P-VALUE OUT OF GG. 
*(D> 1 )/Ll 
INTER*GGZD ] 

a"" ’ OTEERWISE CONDUCT INTERPOLATION TO GET PROPER P-VALUE. 

Ll:E*FFZD-lj-FFZD1 
F*FFZD-ll-A 
G*GGZD1-GGZD-ll 
?L*(F*C)*E 

INTER*GGZD-ll*PL . 

-*0 

L2: INTERS 1 
7 



V RENP+RENDALP NiA:BiC:DiE;NN:X:XX:AAiF\TiG;P:TPL 
P THIS FUNCTION IS A SUBPROGRAM OF KENDALL'S B (KEN) AND NON - 

p PARAMETRIC LINEAR REGRESSION (NPLR). IT CALCULATES TEE CUMULATIVE Q 

p DISTRIBUTION OF B FOR A SAMPLE SIZE N. 

X+ 1 1 

p INITIALIZE FREQUENCIES FOR X FOR SAMPLE SIZE NN . 

NN+2 

P DETERMINE SIZE OF RIGHT PROB . TAIL VECTOR 

<7«-(L((JYx(JV-1))+4))+i 

p OUTER LOOP INCREMENTS THROUGH THE N SAMPLE SIZES TILL THE DESIRED ONE IS 
GENERATED 
Ll\D+pQ 
O 0 
F+oX 
NN+NN+1 

B+((NN*(NN- 1))*2)+1 

p INNER LOOP GENERATES JW+1 FREQUENCIES FROM THE VECTOR OF NN FREQUENCIES 
L2:A+A-1 
C+C+X C (B-A )] 

D+D.C 
AA + 0 



p IF SIZE OF 0 EQUALS NN AND INDEXES OF X STILL REMAIN GO TO L 4 

+(((pD)=NN)a((B-A)<F ) )/£4 

P OTHERWISE CONTINUE TO INCREMENT THRU OLD FREQS 

+ ((B-A)<F)/L2 

p WHEN LEFT HALF OF NN+ 1 FFCrOP J5 COMPLETE GO TO L 5 

p TFXS Z00P ALLOWS ONLI NN TERMS TO 3E USED TO GENERATE NEW VECTOR 

Z- iA*-A “1 
AA’*~AA+‘l 

C+C+XZ(B-A)1-XZAA1 

D+D,C 

*((U(B-l)*2)+l)>(B-A))/LU 
a INVERT VECTOR D AND ASSIGN TO E 

L 5 : E*eD 

a FOR NN OF APPROPRIATE SIZE EITHER DROP FIRST VALUE OFF E OR NOT 
* (NN=3 ,6, 7, 10, 11, 14, 15, 13, 19) /L3 

E*l+E 



a COMPLETE NEW VECTOR OF FREQS X BJ CONCATENATING D WITH E 
L3:X*D ,E 

a n CONTINUE UNTIL SIZE OF SAMPLE N IS REACHED 

*(NN<N)/L1 

a GENERATE VECTOR OF CORRESPONDING P STATS OF PROPER SIZE C 



101 



[43 

44 

45 

48 

49] 

503 



,Q ? CALCULATE B STATS FROM THE P VECTOR 
T+\ (((4xP)+(iyx(i7-i)))-i) 

p TAKE ONLY G ENTRIES FROM TEE FREQUENCY VECTOR 
TPL+T*(Nx(N-l) )t2 
XX+G+X 

p CHANGE FREQS TO CDF VALUES AND OUTPUT B STATS W/ APPRO. 
KENP+ (3 ,G)pT , (T+\XX)*('.N)),TPL 



CDF VALUES 



V MAN+N MANWP M\?\Q \P\S\T\ U : V:B xNN:N;M : UU; MM 
[1] p THIS FUNCTION IS A SUBPROGRAM OF THE MANN -WHITNEY TEST (MANN). 

“2] P IT GENERATES THE CUMULATIVE DIST. FOR THE U STATISTIC. N = SIZE OF 
;3] p LARGER SAMPLE ; M = SIZE OF OTHER SAMPLE . 

* 5 J ^COMPUTE NUMBER OF TERMS TO BE INCLUDED IN LEFT TAIL DISTRIBUTION LESS 1 
[61 MM+(l ( OVxM)*2 ) ) + l 

[7] p SET F VECTOR EQUAL TO 1 CONCATENATED WITH MM ZERO'S 

[83 F+l,MMpQ 

[9] p SET P EQUAL TO THE MINIMUM OF N+M OR MM 

[10] P+L/((N+M) ,MM) 

[ill p SET Q EQUAL TO THE MINIMUM OF M OR MM 

[12] Q+l/MjMM , 

[13] P cb TO LINE DENOM IF MM IS LESS THAN N+USIZE OF X+l') 

[14] +DENOM*\ (MM<N+1 ) 

[15] P IF MMZN+1 GENERATE FIRST BLOCK OF RECURSIVE RESULTS USING NUM LOOP 

[16] NUM-.T+N + 1 
[ 1 7 [ L2:U+MM+1 

[18] p PRIMARY FORMULA USED IN GENERATION OF FIRST BLOCK OF RECURSIVE RESULTS 
fi Q ] £i:*Tcn«- 7r £n -ftu-ti 

[20] p ASSIGNS' NEW DECREMENTED VALUE TO U AND TESTS IF T < THIS NEW U 

[21] +L1* \ (T<U+U-1 ) 

[22] +L2*\ (P2T+T+1 ) 

[23] p GENERATE FINAL RECURSIVE RESULTS USING DENOM LOOP 

[24] DENOMiS* 1 

[25] 34 : 7«-S+l 

[ 26 1 p PRIMARY FORMULA USED IN GENERATION OF FINAL RECURSIVE RESULTS 

[27] Z2:?LVZ +FZV1 -rFtV-Sl 

[28] ^33xi ( (MM+l)S7«-7+l) 

[29] -►34x 1 (QTzS+S+l ) 

[30] P CONVERT FREQUENCY TABLE TO CDF VALUES FOR FINAL OUTPUT 

[31] MAN+(+\F)*(N l (W+M)l 
7 



1 

;i°: 

{l2, 

,13. 

[14] 

[15] 



V Z+NORMCDP X:A\B-,C;D 

A EVALUATES NORMAL CDF AT VECTOR X. FOR |X<4, 26.2.11 JW ABRAMOWITZ AND 
a STBCOT7 - P. - 932“I5 - ffS£B.~FWrL;iBc;B XT THE CONTINUED PRACTTOWTN" W ALL", 
AP. "3577 92.11, iS"VSED"AT DEPTTTZB". APPEAR S"T0-GIVE" AT" LEAST" 1S"S TONI - 
rFICANT FIGURES. "PORTED TO MAINFRAME 17J^~,~ZINE ZB1~HAS BEFIT C BANGED" TO 
TTlTOTD UNDERFLOW PROBLEMS WITH x\ 

*Uo ,A*C\X<H)/Z*X*,X)=p ,X)/3+0 LC 
B*Uxzn)/x 
-*-((p,B) = p X)/BIC 

4-«-0, 5+C*i?x((*d*2)xo2)*0.S )x+/x\ CA*2 )® . + l+2x\ (C-^IOT F10xr/ 4)+l 

zc ' - 



:[(Ix<4)/idX]^' 
»((p,A)=pX)/g. 



BIG:C* 16589790 56295540 S20S0600 19934640 3680160 341952 15232 256 
D* 2027025 32432400 75675600 60540480 21621600 3843840 349440 15360 256 
B-«-l- (B + 2 x ( (o2)x*B*2)*0.S)x(+/((0.5x,B*2)®.*0 > x7)x((p,B),8)pC) + +/((0.5x,B* 



2)®.*0,t8)x((p 5),9)pD 

ZC(|X24)/ipX]-«-$ 



fli 

’7 ' 

.9] 
, 10 . 
. 11 . 
12 
‘l3 ' 
L 14 , 
.15, 
.16, 
,17„ 
-IS] 

.19] 



V Z+NORMPTH P:A:B iC;DiQ;T;SiRiF 

A IMPLEMENTS ALGORITHM AS 111 BY BEASLEY SPRINGER. APPLIED ST AT , 1977 
a FOR A VEGTOR~TNPUT~OF FRACTIONS ."EETURNS CORRESPONDING" NORMA Z~OuAN- 
a TILE ST CLAIMED AGGURAGY TS~ BETTER TRAN~T~.S*ZU"E~.-TOR~GREATER"AcOURACY , 
a ESPECIAZZT--OR-F-TRENE P" VAZGEN", ADD" ONE OR MORE NENTOTTP.APR SON" LOOPS . 

■'‘^/IP^U ) ,(P^Z)17ERR 

-<v/(( | Q->- ,P-Q . 5 )S0 . 42 ) )/3+0LC 

S+-Z +- , Q 

+EXT 

T+(Q .42210 )/Z«-.C 

+ (.P+((p ,T}=p ,p))/2+CLC 
S*-(.0 . 42< I O /Q 

A* 2-50662823884 18.6150006252 

B*- 8.4735109309 23 . 083^6743743 



41.39119773534 ' 
21.06224101826 



T+Tx ( ( (T*2 )® . *0,i3)+.x4)+l+( (T*2 )®.*i4)+.xB 
Z[ (0 . 422 1 <J )/ " 



'25.44106049637 

3.13082909833 



J/ip.4]«-r 

*lF= 1)^0 

EXT-.C+- 2.78718931138 2.29796479134 4.8 501412713 5 ’2.321212768 58 

D* 3.54388924762 1.63706781897 

S<-(xS)x({po.,o.l3) + .xC) + l + (((B^(|*0.5-|S)*0.5)®.* 1 2 )+.xD) 
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I 



>o] zcco.42<iQjApa]^s 

22] ERR\'OyZ OR MQZZ i ZMiQIS ARE OUT OF RANGE . 




illj 



7 P+PERM N:X;Y;Z 

p THIS FUNCTION IS CALLED BY SPEARP . 
p PERMUTATIONS OF N RANKS . N= SAMPLE 
-0 *\N=?^ 1 1 pi 
Z+PERM N ~ 1 
P+- iX+Q 

Ll:-*0x\J7<X«-X+l 

)eX)\Z 

M(x*iiv-i),jnp(,P),,i 

♦LI 

7 



IT GENERATES ALL POSSIBLE 
SIZE. 



Cl] 

3’ 

" 4 * 

Cs" 

* « a 

■o< 



1°1 

113 

121 



F*gJ 

Eisi 



7 JMV 30A//C QiI;J;K:L;M;U 

p TPJS FUNCTION GIVES A CHOICE OF NONPARAMETRIC CONFIDENCE INTERVALS 

p pop tpp gyp ommtlz op a continuous population . n=sample size, 

p Q‘ QUANTILE. IT CALLS THE SUBPROGRAM CBIN . 

p CENTERS CHOICES AT 'THE * ORDER STATISTIC NEAREST ESTIMATE . 

1+ LO . 5^iV xg-ri 20 

p c/ IS 1.5 STANDERD DEVIATIONS (APPROX . ) 

. 5t 1 . 3x (Jxi-^rlOC )*0 . 3 

p 1 RUN 1 OUT K BOTH WAYS. 

K+ l+tl+L/(/C23 JV-IClJ ) 
p CONFIDENCE COEFFICIENTS 
LW$(2,pP)o(tf .C.OlxC) CSJ^((Jl1]+P-1), (JC23-P)) 

Ji- 1 ORDER STATISTICS I COEFFICIENTS ' 



p p<apr op pop^Trjivo output n 

EO.i: 29o ( (loo f -» ) , f * ' ) 

o- 7 o 6 o 3<5(2 ,p/:)o(i-t((j:2:-/0,jci: 



7 



,((pP),4)0‘ ’ ’’T > x 

,C13 U, 12 6 *((pS),l)pL 



+X-1)) 



15 



5 . 

7 

8 , 

9] 

101 

11J 

12 J 



7 SPN+N SPAPROX X:Y 

P THIS FUNCTION IS A SUBPROGRAM OF SPEARMAN'S R (SPMANP) 
p IT APPROXIMATES THE CUMMULATIVE PROB FOR R WHEN PASSED THE SAMPLE 
p SIZE IN THE LEFT ARGUMENT AND THE ABSOLUTE VALUE OF R IN THE RIGHT 
p ARGUMENT . SUBPROGRAMS OF THIS FUNCTION INCLUDE : TDISTN 

p CALCULATE THE CONTINUITY CORRECTION 
y-*-6rWx~i+A7*2 

p TRANSFORM THE STATISTIC R INTO ONE THAT CAN BE USED WITH THE STUDENT 
p T DISTRIBUTION 
X>(X-y)x((JV- 2 )*l-(X-7)*2)*0. 5 
p CALL THE T DIST FUNCTION TO CALCULATE THE P-VALUE 
SPN+1- (N-2 ) TDISTN X 



7 SPP£P«-SPP4PP W;Cl;A;P:CiL:P:M:/V:Ll;L2:LJM:P;CLP 
[11 P THIS FUNCTION IS A SUBPROGRAM OF SPEARMAN'S R (SPMANP ) . IT 
12 J p CALCULATES THE EXACT CUMULATIVE DIST. FOR R FOR THE SAMPLE SIZE 
[3] p PASSED AS THE RIGHT ARGUMENT. BECAUSE OF THE LARGE COMPUTER MEMORY 
p REQUIREMENTS , N IS LIMITED TO SIX ON THE PC AND 1 ON THE MAINFRAME. 

£5] p SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: PERM 
[ 6 1 A+B + 0 

C7 1 p INITIALIZE VARIABLES, VECTORS, AND MATRICES. 

£91 Cl+D l-o 0 
,9. ;*♦(//, :V)oO 

- q i (2—p - (2/- * ) 

!ll[ a IRIS ZCd? GENERATES AN N*N ARRAE OF TEE POSSIBLE VALUES OF DIFFERENCES 
12 a BETWEEN ANY TWO PAIRED RANKS BETWEEN SAMPLES. 

,13, Ll:B+3+l 
14. MiB-.l+C-A 
[15, A*A+ 1 
[16, *(B<N)/L 1 

17 a NOW CALCULATE TEE SQUARES OF ALL POSSIBLE DIFFERENCES. 

[18, M+M* 2 

19 a CALL PERM TO LIST ALL POSSIBLE PERMUTATIONS OF N NUMBERS 
, 20 , D+PERM N 

[21 a CALCULATE SIZE LIMIT OF FINAL VECTOR OF R STATS., 
l22j LIM+1 +(((N*3)-N)*12) 

i,23 j a INITIALIZE VALUES BEFORE INDEXING OUT COMBINATIONS OF ALL POSSIBLE 
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[24 

h 25 
,26 
C 27 
28 
29 
3 0 

31 * 

32 

33 

2u1 

3 5 ^ 

36 * 

37 
38] 
3 9* 
40 
41] 
U2] 

C u 3 J 
L 44 J 



p SQUARED VALUES. 

A+0 

E+{N,IN) P0 

p THIS LOOP CALCULATES ALL POSSIBLE COMBINATIONS OF THE SQUARED VALUES. 

sd.A*MZA;DZ;A11 

*lA<N)/L 2 

fl ADD DOWN ALL ROWS FOR EACH COLUMN TO SUM UP SQUARES COMBINATIONS . 

D2*+VE 

a ADD UP NUMBER OF DUPLICATED SUMS OF D-SQUARED VALUES AND COMPRESS 
s VECTOR DOWN TO UNIQUE VALUES. 



L3:Cl*Cl,(.i/D2 ) 
D1+D1.U/((L/D 

D2+Ul/D2)*D2)/D2 



: (t /D2 )-D2 ) ) 



( (pCl )<LIM) / L3 

a TRANSFORM SUM OF SQUARES VALUES TO SPEARMAN ' S R STATISTIC. 
R+l- (6xCl )*(Nx ( (W*2 )-l ) ) 

fl CALCULATE CDF VALUES ASSOCIATED WITE TEE R STATISTIC. 
CDF*(+\D 1 )* liV 

a FORM TWO ROW MATRIX FOR OUTPUT OF R STATS AND CDF VALUES. 
SPEAR * ( 2 . ( pCl } ) p (R , CDF ) 

V 



I 



12 j 

EiS 

if: 

:ii 

:IS: 

22 

[ 23 , 

.24, 

25 

[ 26 ] 

[27 

. 28 , 

.29. 

:If: 

:1s; 

:If: 

II; 

.40. 

:«h 

43 

’45j 

’U7‘ 

43] 

UQ] 

’SO] 

5 1] 

52 ] 

53 ] 



7 SPM+X SPMANP 7iC‘,DiDD;Dl;D2;N;DEN0MR;XX;Il;NSiNUMR;P;PV;PVALiSU;SViRB0 s 
U‘V'ARR0’U1'V1'X‘I‘WW 

o ’ ’ THIS ’ FUNCTION Is A SUBPROGRAM OF NONPAR LINEAR REGRESSION {NPLR ) 
a AND SPEARMAN'S R ( SPMAN ) . IT COMPUTES THE SPEARMAN R STATISTIC 

fl p-valuss. the left argument that is passed is the X 

a SAMPLE; THE RIGHT ARGUMENT IS THE I SAMPLE. 

fl SUBPROGRAMS OF THIS FUNCTION INCLUDE; TIES, TIESK, SPEARP , SPAPROX , 
fl INTERP, AND THE VARIABLE PMATSP . 



V+IZ&X1 

*U+XZ&X1 



ORDER I IN INCREASING ORDER OF X 
ORDER X IN INCREASING ORDER 
COMPUTE CURRENT RANKING OF I 
NOW ORDER I RANKS IN INCREASING ORDER 



D1+VZ&V1 

P IF TIES EXIST IN EITHER X OR I RANKED VECTOR USE MID -RANK METHOD 

DD+1 TIES D 1 
XX-<-l TIES U 

P FIND ORIGINAL RANKING OF J WITH TIES RESOLVED 

Zl+DDZCl 

P RECORD SIZE OF INPUT VECTOR 
N+ pX 

P CALCULATE DIFFERENCES BETWEEN RANKS OF X AND I VECTORS 
D+XX-Y 1 

p DETERMINE THE SUN OF SQUARES OF THE DIFFERENCES 
D2++/(D*2 ) 

P OBTAIN THE NUMBER OF TIES IN EACH VECTOR USING THE TIESK FUNCTION 
Ul+TIESK U 
Vl+TIESK D1 



((+/< 




)- C+/U1 ) 


)*12 


((+/< 


[71*3 ) 


j-(+/Vl) 


5 + 12 



NS+Nx ( (N* 2 ) - 1 . 

P CALCULATE TQE R STATISTIC INCLUDING THE CORRECTION FOR TIES 
NUMR+(NS)+(( 6 )*D2 ) + ( ( 6)x(Stf+S7)) 
DENOMR+UNS-(l2xSU))*0.5)xUNS-(l2xSV))*0.5) 

RHO+NUMR* DENOMR 

ARHO+\RHO 

+{N26)/L1 

p CALL SPEARP TO CALCULATE THE RIGHT TAIL OF THE CDF OF R 

P+SPEARP N 
+L2 

Ll:+(N>10)/L3 
P+PMAT SP t (iV - 5 ) : : ] 

p CHANGE SIZE OF P TO AN MxN MATRIX 
P+PZ 1; :] 

p CALL INTERP TO CALCULATE P-VALUE BY INTERPOLATION 

L2 iPVAL+ARHO INTERP P 
-<PVAL=~i)/L7 
PVAL+Q . 5 
~L 7 

p CALCULATE P VALUE USING STUDENT T APPROX. 

L3:PVAL+N SPAPROX ARHO 

L7:SPM+(RHO),PVAL 

V 



1 04 



7 P+K TDISTM X',7 
fl THIS FUNCTION IS A _ _ 
b POP SPEARMAN ' S R (SPEASP) 
s TEE STUDENT'S T VIST WI r E K 
n THIS FUNCTION CALLS ON TEE 
V+(X+,X)20 

P-*-0.5x(l JO FDISTN X*2 
*3x iO=v // 

P[V/ipX]*0. 5+7/P 
* 0 x; 1 =a /7 

?C(~7)/ipX3+0.3-<~7)/P 

7 



SUBPROGRAM OF TEE CUMULATIVE PROBABILITY GENERATOR 

Jr CALCULATES TEE CDF AT X USING 

DEGREES OF FREEDOM. 

DISTRIBUTION FUNCTION (.FDISTN). 



iis 

"a 

. 10 ] 

tiij 

f 14 1 

r 

!is3 

rial 

C 20 J 

m 

*2U* 
25! 
26 
27! 
.23, 
29 
*3 0, 

31, 

32, 

33, 
3 4 

35, 

36, 

ta 3 7 

.33, 
.39, 
[40. 
[41 
.42 
*43 
,44 
45 
>6 
47 
*4 3 

49 

50 

51 



r 



7 TI+33 TIES B;C;D;I;N;T;Y;Z;K:M;L;PP;NR;MM 
fl THIS FUNCTION IS A SUBPROGkAM OF KENDALL'S B (KEN). SPEARMAN'S 
fl R ( SPi’lAN? ) , KRU SEAL -WALLIS ( KRWL ) , MANN-WHITNEY (MANW) 
fl ^77 WIZCOXON (WISIG). IT CHECKS THE RIGHT ARG . VECTOR FOR TIES AND 
fl CHANGES THE TIED POSITIONS OF THE LEFT ARG. 31 THE MIDRANK METHOD. 

N+c3 

fl J7 iV0 VECTOR OF RANKS IS PASSED; GENERATE ONE 

+ ( (poBB )-0 )/76 
33+ \N 
L 6 : I+pO 
7«-iVp0 

a* * CHECKING FOR TIES 31 INCREMENTING TBRU THE VECTOR 

^3 . ) = 3 " TT " ! 

' ” COUNT NUMBER OF TIES ; 17 M2 TIES GO TO Z2 

D+- / C 
+(D= 0 )/Z2 

fl RECORD WHERE TIES STARTED AND HOW MANX RANKS INVOLVED 
I+I ,T t (D+ 1 ) 

fl INCREMENT NEXT T 31 THE NUMBER OF TIES ENCOUNTERED PLUS 1 

L2iT+T+(D±l ) 

flJ7 T LESS THAN SIZE OF ORIGINAL VECTOR GO TO L 3 AND START AGAIN AT NEW T 
~(T<N)/Z 3 

fl* * ASSIGN THE RANKS OF THE LEFT ARG. TO TI 

TI+BB 

fl IF NO TIES FOUND QUIT 

+(i- 0)/0 

z^o 

fl LOCATE THE INDEXED POSITIONS OF TIED RANKS 
LS :PP+( ( J[l+Z] )-l)+i((ICl+Z3+J[2+Z] )-l) 
fl 72W7 777 VALUE OF THESE RANKS 

NR+(+/TILPP2 )+I[2+Z3 

fl SET UP VECTOR WITH ZEROS AND ONES; ONES WHERE TIE RANKS INVOLVED 
K+0 

MM+NpO 

L*:K+K+1 

M+(TI=TI[PP[KJJ ) 

MM+MM+M 

-►(7< J[2 + Z] )/Z4 

fl SET UP A VECTOR WITH ZEROS WHERE TIED RANKS OCCUR 
L+~MM 

fl TRANSFORM ONES OF MM VECTOR TO MIDRANK VALUE 
MM+NRxMM 

fl TRANSFORM ONES OF L VECTOR TO REMAINING UNCHANGED RANKING VECTOR 
L+TIxL 

fl FILL IN MIDRANK VALUES 
TI+L+MM 
Z+Z+2 

fl DO THE SAME FOR ANT OTHER TIES INVOLVED BUT WITH NEWLY COMPUTED TI 
+(T>Z)/L 5 
7 



7 TIE+TIESK A A xAA :B ;C;D; I: N;T 

THIS FUNCTION IS A SUBPROGRAM OF KENDALL'S 



4 ] 

I 

. 10 ] 



__ _ , 3 (KEN), SPEARMAN'S 

R ( SPMAN ) AND (SPMAN1 ) , AND KRUSKAL -WALLIS t KRWL ). IT CHECKS THE 
RIGHT ARGUMENT FOR TIES AND RECORDS THE NUMBER OF OCCURENCES OF EACH 
TIE AND THE TOTAL NUMBER OF TIES IN THE VECTOR. 



5 ASSIGN ORDERED 

N+oAA 
B+AAL&AA1 
TIE+p 0 
T+l 



VECTOR TO B AND INITIALIZE VALUES 
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fll] a CHECKING FOR TIES 31 INCREMENTING THRU THE VECTOR 

12] L3:C<-(T+3) = 3lT1 

[13] a COUNT NUMBER OF TIES ; IF NO TIES GO TO L 2 

14] D++/C 

15] +(D=Q)/L2 

16] TIE+TIE, (P+1 ) 

17] o INCREMENT NEXT T BI THE NUMBER OF TIES ENCOUNTERED PLUS 1 

[ 18 j L2 :T+T+(D+1 ) 

19] pJF T LESS THAN SIZE OF ORIGINAL VECTOR GO TO L 3 AND START AGAIN AT NEW T 
[20] +{T<(N-1))/L3 

7 



7 VAR+VARMW 3 \C ;D ;E ; D1;E1 

p THIS FUNCTION IS A SUBPROGRAM OF THE MANN -WHITNEY TEST (MANW ) . 

p IT GENERATES THE RANKING SCHEME USED IN CALCULATING THE 
p DIFFERENCES IN SCALE (1 ASSIGNED SMALLEST , 2 ASSIGNED LARGEST . 3 
p ASSIGNED NEXT LARGEST . 4 SECOND SMALLEST .ETC . BY TWOS TILL PROPER 
p SAMPLE SIZE IS REACHED). SAMPLE SIZE IS PASSED IN THE RIGHT ARG. 

E+El+p 0 
D+0 

p FIND FLOOR OF MIDPOINT OF VECTOR AND ASSIGN TO C 

C+L (5*2 ) 

p LOOPS GENERATE RANKING VALUES LEFT HALF FIRST 
L2-.D+D+1 
E+E ,D 

+((qE)=C)/L3 
D+D + 2 
E*~S , D 

-*■( (oE)<C)/L2 

NOW GENERATE RIGHT HALF 

L3:Dl*i 

Dl+Dl+1 

ei+ei ,m 

■*■( (sEl )=C)/L5 
L6-.D1+D1 + 1 

*t(oEi)lc)/L5 
Dl+Dl+3 
El+El yDl 
-►( (pEl )<C)/L6 

p IF SIZE OF VECTOR IS 'ODD VALUE CONCATENATE MIDDLE RANK IN BETWEEN HALFS 
L5:-K(2T2)*o)/L7 

p IF SIZE IS EVEN CONCATENATE LEFT HALF WITH THE REVERSE OF THE RIGHT 
VAR+E, (BE 1 ) 

L7 : VAR+E, B, (051) 

7 



7 WIL+WTLP NN * N * A * P * T * NN % W * PP * NM 

a THIS FUNCTION’ ik 1 SUBPROGRAM OF THE WILCOXON SIGNED RANK TEST 
a (WISIG). IT GENERATES THE CUMULATIVE DIST. FOR THE TEST STATISTIC 
A THE GENERATOR USES A RECURSIVE FORMULA. NN= SAMPLE SIZE. 

JVM-KL((+/tAW)*2 )) + l 
N* 2 

a SET P EQUAL TO PROB. DIST. WHEN N EQUALS 2. 

P-*-4 pi 
L3-.N+N+1 
A* 1 

a SET T VECTOR TO PROPER SIZE OF ZEROS. 

r+(+/(iW))pO 

a IF ASN USE TRUNCATED RELATION TO COMPUTE OCCURRENCES. 

LU: + (.ASN)/L1 

a IF A>N USE FULL FORMULA TO COMPUTE OCCURRENCES. 

*(A>N)/L2 

a WHILE ( A-N ) IS NEGATIVE TRUNCATE FORMULA TO AVOID A NEGATIVE INDEX. 
Ll’.TZAl+PZAl 
A+A+l 
*L 4 

a IF A IS LARGER THAN THE LENGTH OF P GO TO LB. 

L2:+(A>(oP) )/L6 

a ONCE (A-N) BECOMES POSITIVE; THE RECURSIVE FORMULA CAN BE USED. 
TZAl-Pi(A-N)j-PiA: 

+L7 

a ONCE A IS LARGER THAN THE LENGTH OF P TRUNCATE FUNCTION AGAIN. 

LS-.TZAl+PZZA-N)! 

L7:A-*-4 + l 

a IF A AS AN INDEX HAS NOT EXCEEDED N(N+1 )/2 GO AGAIN. 

■*(AS ( +/ ( \N) ) )/L4 

a CONVERT T INTO P AND CONCATENATE 1 FOR USE IN NEXT ITERATION OR OUTPUT. 
LS:P*((+/(iJV))pD,l 
PP+NM+P 

WIL+(+\PP)*(.2*NN) 

a CHEQK IF LENGTH OF INPUT VECTOR EXCEEDS NUMBER OF N' S GENERATED. 
*(NN>N)/L3 
7 
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APPENDIX G 



LISTINGS OF PROGRAMS USED TO GENERATE C.D.F. COMPARISON 
TABLES 






7 KENTEST-.N: ALPHA iB;A-,C;TAU;P’,NUM; DEN iZ;ZZiERRZZiMiZZC;DiAA\PP;NUMC;I;R;P 
;FS;S;J; ERRZZC \ZC\KR 

sTHIS PROCRAM GENERATES TABLES OF C.D.F. COMPARISONS FOR KENDALL'S 3 



E13 

[10 

[12 

Ub 

m 

; 185 
.19] 
J|20 J 
r 21 < 

c?5: 

[25 

C 26 

*27 

> 28 

C29 

.30 

.31 

.32 

.33 

.34 

*35 

36 

>37. 

38 

*39^ 

40 

41 

£ U2 3 

*43< 

44 

45 

f46; 

I 47 
*43 
>49 
50j 

51, 

52. 

53 * 

54 
55 1 

:l?3 

;Ic: 

*i^ 

>63 

t64 

65 

[66 

*67 

68 

69 

70 
*71 

72 

73 



SET SAMPLE SIZE AND ALPHA VALUES. 

,01 0.02 0.03 0.05 0.03 0.13 0.13 
THIS LOOP INCREMENTS SAMPLE SIZE. 



N+ll 
ALPHA + : 

Ll:PP*pQ 
M+p 0 
XK+pO 
FS+oO 
ZZ*-o 0 
ZZC^o 0 
N+N + 1 
D + C 
3*-G 

fl COMPUTE CUMULATIVE DIST AND ASSOCIATED STATS. 

P+KENDALP j V 

p TF/S LOOP CALCULATES ALPHA VALUES AND APPROXIMATIONS . 

L2:B+3+l 
C+ALP3AZ32 



A*(+/(?t2:]SC)) 

2VU/~?Cl:*l} 



TK* KK.TAU 
AA+6*Nx 1+N*2 
PP*PP ,PZ2;A1 

A COMPUTE NORMAL APPROXIMATION . 

NUM+3x(TAU)*<.<.Nx(.N-1))*0.5) 

DEN+(.2x(. (2xtf)+5))*0.S 
Z+NUM*DEN 

ZZ+ZZ ( 1 “NORMCDF Z) 

A ’ COMPUTE NORMAL APPROXIMATION WITH CONTINUIT1 CORRECTION. 

NUMC+3* (.TAU-AA )x ( (tfx (tf-1 ) )*0 . S ) 

ZC+-NUMC *DEN 
ZZC+-ZZC , (1 -NORMCDF ZC) 

*(B<7)/L2 

a COMPUTE ERROR DIFFERENCES. 

ERR ZZ+PP-ZZ 
ERRZZC+PP-ZZC 

fl PRINT OUT TABLE OF VALUES. 

1+0 . 1 
•N-19 5 

S+'TEST ST AT. VALUE « 

M+M,' PROBIB 2 g]; FOR SAMPLE SIZE EQUAL TO » , (2 0 *N),' 

+L3 

I5:t7^194 
L3:M+M. [J] 

♦ (Z? = 4)/£4 
1 18 pS 
-»(£>*0)/£6 

F«-'F4C<T 9.9999 >» QFMr(l 7 p0) 

F<- 1 62 pF 
FS+H,F 
+L7 

LS:FS+pO 
H+l 7 oS 
FS-FS.F 
C-Q 

Z3 : C^C-1 

I 1 , (7 5 sPPCC]) 

FS+FS.F 
+(C<7)/L9 
L7:M+M,lll FS 
D+D+l 
„ , ^1 

.671 ♦(iVl f JV2 1 JV3,£5)[03 

[68] flfi:s+ f rxicr 1 

>7-*-198 
*L 3 

N2:S*-' ERROR ; NORMAL ' 

PP+ERR ZZ 
*L3 



Q57C(l8pl97 ),«7,61p(8pl97 ) ,«7] 



1 07 



N3:S+' ERROR: NORM . W/CC « 
PP+ERRZZC 
+L 3 
L*:M 

+(N<1*)/L1 

7 



7 KWTEST ; A:3;C ; ?P ; PC ; NN : ALPHA :N;K;P\F;E; PVALUZ ; PVAL ; PF 1 ; PF ; P7F ; P7F1 ; D ; Pflfl 
B;ERRF;ERRF1;CDF;CC;KK;I;J;FF;S;M 

p TPJS PflCCflAM CfliVPMrflS TABLES of c.d.f. comparisons for tee eruskal- 

p WALLIS TEST . 

^444 

CC«-2 

4LPfl4-<- 0,01 0,02 0.03 0.05 0.08 0.13 0.18 
+L2 

LliA.+NN 

p CALL KRUWALP TO GENERATE EXACT DISTRIBUTION FOR SELECTED SAMPLE SIZE . 
L2:P+KRUWALP A 
3+0 
N++/A 
CC+CC+1 
K+pA 
M+pO 
PP+oO 
KK+o 0 
PC+o 0 
PF+oO 
PFl+oO 
CDF^P [2;] 

L^:B+B+1 

p LOCATE POSITION OF EXACT CDF VALUE £ ALPS A. 

D + ( + / ( CDF&ALPHA CP ] )) 

+(D=0)/L3 
D+ 1 

p DETERMINE CORRESPONDING TEST STATISTIC VALUE. 

L2:3+Pll;Dl 
KK+KK E 

R RECORD EXACT VALUE OF CDF. 

PVALUE+P [ 2 : Z>] 

PP+PP,P74Ltffl 

fl COMPUTE CORRESPONDING P -VALUE USING CBISQ APPROX . 

PVAL+U- 1) CflJSC fl 
PC+PC, (l-P^L) 

p COMPUTE CORRESPONDING P-VALUE USING F APPROXIMATION 
F+UN-K)xH)*C(K-l)x((N-l)-in) 

PVF+UK-DAN-K)) FDISTN F 
PF+PF. (1-PVF) 

fl COMPUTE CORRESPONDING P-VALUE USING F APPROX W/ 1 LESS D.F . IN DENOM 
PVF1+(UZ-1)A(N-K)-1)) FDISTN F 
PFl+PFl , ( 1 -PVFl ) 

-*(fl<7)/£4 
D+ 0 

ERRH+PP-PC 

ERRF+PP-PF 

ERRFl+PP-PFl 

p PRINT OUT TABLE OF VALUES. 

I+'0 . 1 
J+ x 

S+'TEST ST AT. VALUE 1 
M+M, ' PflOflCfl ;> x UAVm6lA2; 

, 4 , AND 1 , (sCC), 1 OBS. 1 

+L 8 

LS:J+'-' 

L8 : M+M . [ J] ( 17p 1 - 1 ) , J , 61p ( 8p 1 - * ) f J 
+ (D=5 )/L9 
F+ p 0 
fl^!6pS 
F+F ,B 
C+ 0 

^Cfl*0)/£10 
Cll :C-K>1 

?p- ' Y 1 , (7 4 ) 

F+F »FF 
+ {C< 1 )/L 11 
+L12 

L10 lC+C+1 

FF+' M, (7 5 vKKLCl) 

F+F t FF 
+ (C<7 )/L 10 

L12:M+M,IH F 
D+D+l 
1 + 1 

+ <in.W2 J73,tf4,L5 )[Z>] 
tfl.’S+ULPfl^ 74LCP 1 

KK+PP 



POP A CPC C/P CP 3 S4MP£PS CONSISTING ■ OF 4 



1 08 



C7S] +L 8 

80l N2:S+' ERROR; CEI SQUARE' 
811 KK+ERRB 
8 2 J ~*L 8 

'S3 j N3 : S* 1 ERROR ; F DIST ' 
]8u] KK*ERRF 

85] *L8 

86] N^iS*' ERROR; F V/ 1 OF' 

87] KK*ERRFl 
, 88 ] *L 8 

[ 39 ] L3-.M 
L90j *(CC<3)/L1 
7 



7 KWTESTSM;A;R2;PC;NN;ALPBA;N;K;1;W;B;TT;PVAL;PF1;BB;MM;M;X;ERRC;ERRF1;F; 
PV;B;KK;b; PvFl ;Q;FF ;U;C;I;J;S 

fl TBIS PROGRAM GENERATES C.D.F. COMPARTSONS FOR TBE KRUSKAL -WALLIS TEST 
a 3 AS ED ON 300 OC RANDOM VECTORS OF U 0 RANKINGS WHICH SIMULATE COMPARING 
fl 5 SAMPLES OF 8 OBSERVATIONS EACH. 

A* 2 

W * 10000 
K *5 
HB*oO 
N* 4 0 

NN*0 

] 1 2 ] Ll:NN*NN*l 

Cii] MM * 0 

.15] a GENERATE A RANDOM VECTOR OF RANKINGS. 

]l 6 . X*40?40 

LI 7] I?*d0 

[18] fl SUM TBE RANXS OF EACB SAMPLE. 

]l9] L 2 :Y**/(MM+(.Q + X')) 

. 'V+rv v 

riil «*c*s’ 

’22] MM*RM J - 3 

[23] *(CS40)/L2 

C24] J?2*+/((XT*2 ) + 8 ) 

25 a CALCULATE THE B STATISTIC FOR EACB VECTOR OF RANKINGS. 

C 2 6 J B* ( ( 3*410 )xi?2 ) ”12 3 

27 BB*BH.B 

28 *( (flW*10000 ) = 1 , 2 ,3 ,4 )/L4 

[29 fl CBANGE RANDOM SEED EVERT 10000 ITERATIONS. 

[30, URL* 0 60 60 1000 l"4+DrS 

31] L'i:*(NN<W)/Ll 

32] b ORDER ALL H'S FROM LARGEST TO SMALLEST. 

[33] BB*BHIVHB1 

[34] B * 0 
35] D * 0 

.36] M*o 0 

37] ALPBA* 0.01 0.02 0.03 0.05 0.08 0.13 0.18 
[38] KK * pO 
39] PF l*oO 
.40] PC * DO 
41 L3:B*B+1 

>2] fi DETERMINE DESIRED B VALUE CORRESPONDING TO ALPBA VALUE. 

43] E*HB'ANN*ALPHAiBl )] 

* 44 ] KK+KK R 

|45 j b COMPETE CORRESPONDING P -VALUE USING CBISQ APPROX. 

46 PVAL*(K- 1 ) CBISQ B 

47 PC*PC. (, 1 -PVAL) 

48] b COMPUTE CORRESPONDING P-VALUE USING F APPROX W/ 1 LESS D.F. IN DENOM 
> 9 ] J?*((W-X)xffU(4x(39-Hj) 

50] PVF 1 *((K- 1 ), \{N-K)- 2 )) FDISTN F 

[ 51 ] PF 1 *PF 1 .( 1 -PVF 1 ) 

[52] *(B<7 ) /L3 

.53. ERRC*ALPBA-PC 
54] EP.P.F1*ALPHA-PF1 

.55, s PRINT OUT TABLE OF VALUES. 

.56] 1*0.1 

*58* 3+' TEST ST AT . VALUE 1 

159! M-M.' FRO 3 ZB Z : , ZAVZ^Sl , ’ ] : BASED ON ’,(5 3 GENERATED H"S FOR 5 

SAkFLES OF 8 OBS. EACH. * 

+L 8 

L5:J +' -» 

L8 [J] (17p 1 - 1 ) , J , 6 lp ( 8p 1 - 1 ) , J 

-*(I? = 0 /£ 9 
F«-p0 
H+I 60 S 
F+F,B 
C+ 0 N 
+ (D*0 ) /L10 
Ll 1 : C+C+l 
FF + 1 M,(7 
F+F , Ft 



3 9KKIC1) 



1 09 



+{C<7)/L11 
+L 12 

L10 :C+C+1 

FF+ ' T* , (7 5 xKKICI ) 
F+F , FF 
«*(C<7)/L10 
Ll2:M+M,lll F 
D+D+l 
1+1 

** (Nl r N2 ,N3 ,L5 ) CD] 
Nl:S+ J ALPHA VALUE 
KK+ALPHA 
«/*«- 1 + 1 
+L8 

N2:S+' ERROR; CHI SQUARE' 
KK+ERRC 
+L 8 

tf3:S«- f £7?P0P; P W/~ 1 DF ' 
KK+ERRF 1 
->£8 
L3:tf 

D+ 20000 30000 1 

(A< 4 ) / Ll 
V 



V MANTEST ; U : NNiMMiALPHA :3xA;C;T; ?:NUM ; DEW : Z : ZZ ; ERRZZ ; A? ; ZZCiDiPP : AWMC : 
P : PS : 5 : «7 : 2*2* : T ; Z’i’C ; LC ; OTC ; NUMT ; bENOM ; DENOMC ; ERR A VE ; ERRTT ; ERRTC ; PPPZ 
C ; ERRTCZC ; KK 



[ 3 ] 

L101 

. 11 . 

:i§: 

:is; 

:i7. 

.201 
,213 
, 22 ] 
[2 3 ] 
,24* 
25 
[ 26 , 
, 27 , 
, 28 , 
[29 
.30 
. 31 , 
, 32 , 
, 33 , 
. 34 , 
. 35 , 
,36 
. 37 , 
[38 
[ 39 , 
No. 
Ni- 

. 42 , 

.43 
r ul ' 
*4 5" 

N 6 4 

. 47 ] 

. 48 . 

49 

‘ 50 ‘ 

:si: 

. 52 . 

.53 

. 54 . 

55 

- 

►ila 



COMPARISONS FOR THE MANN - 
SET SAMPLE SIZE AND ALPHA VALUES . 



a rtfis PPDCT4M GENERATES TABLES OF C.D.F. 
A WHITNEY TEST. 



01 0.02 0.03 0.05 0.08 0.13 0.18 
THIS LOOP INCREMENTS SAMPLE SIZE. 



AW- 7 
MM+7 

ALPHA + 0. 

Ll :PP«-pO 
M«-p0 
IC/C^pO 
FS+oO 
Z Z^p 0 
ZZOpO 
rr**po 
TTC+oO 
NN+NN+1 
MM+MM+1 
D+0 
B + 0 

« „„ CUMULATIVE DIST AND ASSOCIATED STATS . 

P*NN MANWP MM 

« c _ „ THIS LOOP CALCULATES ALPS A VALUES AND APPORXIMATIONS . 

L 2 • D’+'D + 1 

C*ALPSA CS3 
A*(+/(PZC)') 

U*A- 1 
KK*KK , U 
PP*PP,P[A] 

« ,, COMPUTE NORMAL APPROXIMATION. 

NUM*U-UNNxMM)*2) 

DEN*( (NNx-MM* (NN+MM+1 ) ) + 12 )*0 . 5 

Z*NUM*DEN 

ZZ*ZZ ,NORMCDF Z 

COMPUTE NORMAL APPROXIMATION WITS CONTI NUITI CORRECTION. 
NUMC*(U+0 . 5 ) - ( (.NN*MM)*2 ) 

ZC*NUMC*DEN 
ZZC*ZZC .NORMCDF ZC 

« n „„ , COMPUTE STUDENT T APPROXIMATION 

DENOM * ( ( (NN+MM-1)- (Z*2 ) )+ (.NN+MM-2 ) )*0 . 5 
T*Z± DENOM 

TT*TT . (NN+MM-Z) TDISTN T 

a „„„F 0M ?, U \ S STUDENT I APPROXIMATION WITS CONTINUITI CORRECTION . 

NUMT+’i NUM ) + 0 . 5 

TC^NUMTtDENOM^' 1 ( ' CEN * 2 } ^*^N+MM-2 ))•(.(. (.NUM+ 0.5 )*2 )*(NN+MM-2 ) ) )*0 . 5 

TTC*TTC, ( (AW+AW-2) TDISTN TC ) 

*(.B<7 )/L2 

* „„ COMPUTE ERROR DIFFERENCES. 

ERRZZ+PP-ZZ 
ERRZZC*-PP-ZZC 
ERRTT*PP-TT 
ERRTC*PP-TTC 
ERRAVE+PP- ((ZZ+2T)+2 ) 

ERRTCZC+PP- ( ( ZZC+TTC )+2 ) 
fi PflJWr OC/Z 1 TABLE OF VALUES. 

1*0 . 1 
«/«-ig5 



1 1 0 



t^N 



I 2’J7£'Tt n VALU * 

X * PROBLU <£ Ul; FOR SAMPLE SIZES N EQUAL TO 1 , (2 0 9NN),' AND 

M EQUAL TO ' , (2 0 *MM) , T . ' 

*^L 3 

Z3 :5m& rj] CU7[(l8pl97),J,61p(8pl97),«;j 
-►a?=8 5 /lu. 

P«- 1 18 dS 
-(£=0)/L6 

F*-' BG< | ZZ9 >■ DF«r(l 7 pM) 

p-** 1 6 2 pF 

FS+3.F 

+L7 

£6 : FS+ p 0 
B+17oS 
FS+FS , P 
OO 

£9:00+1 , 

F+' L , x .(7 

FS+FS.F 
**“ ( C< 7 ) /L 9 
L7:JMi,[l] FS 
00 + 1 

♦(ifl .JV2 .JV3 ,iV4 .NS ,i76 ,JV7 .£5 ) [I?] 
tfl:S«-‘FXlcr C.fi.F. r 

019 8 
+L 3 

7V2 : S+ 1 ERROR ; NORMAL « 

PP+ERRZZ 

Nils*-' ERROR: NORM . W/CC' 

^D4.r33"7/- 



5 ®PP[0] ) 



• 95 1 

961 



iV4 : S<- ' ERROR ; T PTSP 1 
PP+ERRTT 
+L 3 

iV5 :S+' ERROR; T W/CC 1 
PP+ERETC 
+L 3 

N6:S+ l ERROR; AVE T/Z ' 

| PP+ERRA V E 
+L 3 

1 N7:S+' ERROR; AVE TC/ZC 1 
j PP+ERRTCZC 
+L 3 
£4:M 

I •►(ilW<9)/£l 
V 



a; 

!■ 

'4' 

“ 5 * 

oi 

its, 

■:ii: 

ii: 

,1*. 

£15. 

15 

^17 

LI 8 g 

Cis 
, 20 1 
2^ 

23 

24 
[25 
[26 

v 27 . 

23, 
[29 
[3 0, 
[31, 
L 3 2, 
[33, 
[ 3 4 j 

J 35 

*3 6 j 



V SIGNTEST;N;ALPHA;B;A;C;K;P;Z;ZZ;ERRZZ;M;ZZC;D;PP;I;H;F;FS;S;J;ERRZZC;ZC 
; KK 

p TSIS PROGRAM GENERATES TABLES OF C.D.F. COMPARISONS FOR TEE SIGN TEST . 



N * 23 
ALPHA* 



SET SAMPLE SIZE AND ALPHA VALUES . 



COMPUTE CUMULATIVE DIST AND ASSOCIATED STATS . 



01 0.03 0.06 0.12 0.22 0.35 0.5 
p THIS LOOP INCREMENTS SAMPLE SIZE . 

LUPP+pQ 
M+pO 
KK+ P 0 
ZZ+pO 
ZZO-p 0 
N+N + 1 
D+Q 
B + 0 

P+BINOM N 

p THIS £00P CALCULATES ALPHA VALUES AND APPORXIMATIONS . 

£2:05+1 
04£PP£C5! 

0( -/ J ) 

O* - 1 

77+VV £ 

PP+PP* f PlA3 

p COMPUTE NORMAL APPROXIMATION . 

Z«-(F-(0.5xtf) )*(0. 5x(tf*0.5)) 

ZZ+ZZ , (NORMCDF Z) 

p COMPUTE NORMAL APPROXIMATION WITH CONTINUITI CORRECTION . 

ZO( (F+Q. 5 )-(0 . 5x/7) U(0. 5x (tf*0.5 ) ) 

ZZC+ZZC ANORMCDF ZC) 

+ (B<7 ) /tl 

p COMPUTE ERROR DIFFERENCES. 

ERRZZ+PP-ZZ 

ERRZZC+PP-ZZC 

p PPJiVr oar TABLE OF VALUES. 

I* 0.1 



. 7^195 

S+'TEST ST AT. 
M+M, ■ 

*-L3 

L5:<7*- 194 



VALUE 



PROBIK Z Kl; FOR SAMPLE SIZE EQUAL TO ' , (2 0 



L3:M+M.ZI1 
*-(L = 4)/L4 



QA7C(18pl97 ) ,J ,61p(8pl97 ) ,*7] 



H*- i 18 pS 
-(L*0)/L6 
F*- 1 9G < I ZZ9 
F*- 1 5 2 pF 
FS*-H , F 
+£7 

L6:FS*-pO 
H+17pS 
PS+FS ,H 
C*-0 

L9:C*-C+1 

F*-' T * (7 5 sPPCC]) 
FS+FS.F 
-(C<7)/L9 
L7:M+M,Zll FS 
D+D + 1 
-T+-1 

-►(JVl.tf2.tf3,L5)CL] 
tfl:S*- r EXACT C.L.F. 
•7*-198 
*-L3 

N2 :S+' ERROR; NORMAL 
PP+ERRZZ 
-L 3 

7/3 : S+ ' FHHCH ; NORM. W/C< 
PP+ERRZZC 
-L 3 
L4:M 

*-(tf<26)/Ll 

V 



>' UFMT ( 1 7 pEH ) 



£13 

* 

"5“ 

"6" 

[7" 

s 

;9] 

10 

:ii: 

12 

13 

.14* 

15 

16 
17 

is 

19 

L20 

,21 

.22 

,23 

,24 

r 2 5 
!26 
27 
,28 
29 
,30 

:«3 

fhj 

£P 

38 

i:“i 

.39, 

, 4 0. 

[41. 

.42, 

■52- 

:ui 

: 52 ; 

:so: 



7 SPMTEST:iVMLPHA;5:A:C?AA:P:Z;ZZ;FHHZZ:tf:ZZC;L;PP;2’;H;P;PS;S;«/;rT;; 
TC ; ERR A VE ; EHHZT ; ERRTC ; EHH ZZC ;ZC\PC ;RHO ;KK 

a THIS PHCCHAM GENERATES TABLES OF C.D.F. COMPARISONS FOR SPEARMAN'S 

a SET SAMPLE SIZE AND ALPHA VALUES . 

tf*-8 

ALPHA*- 0,01 0.02 0.03 0.05 0.08 0.13 0.18 
a THIS LOOP INCREMENTS SAMPLE SIZE. 

Ll :PP*-p0 
M+pO 
EH*-p0 
FS*-p 0 
ZZ*-p 0 
ZZC*-p0 
rr*-po 
rrc*-po 

N+N+l 
D+ 0 
B*-0 

AA*-6+tfx 1+N*2 

a COMPUTE CUMULATIVE VIST AND ASSOCIATED STATS. 

PI : P+PMATSPlN - 5 : ; ] 

PC*- (PC 2 ; ]*0 )/Pt2 : 3 

a THIS LOOP CALCULATES ALPHA VALUES AND APPORXIMATIONS . 

L2 : H*-H+l 
C^-ALPHAZBl 
A*- ( +/ (PCSC ) ) 

HHC*-PL1 : A] 

KK+KK ,Rho 
PP*-PP,PC2;A3 

a COMPUTE NORMAL APPROXIMATION. 

Z+RHOx ( (N-l )*0 . 5 ) 

ZZ*-ZZ, ll-NORMCDF Z) 

a COMPUTE NORMAL APPHCXZMAriC/7 WITH CONTINUITY CORRECTION . 

ZC*-(HHC-AA )x ( (JV-i )*0 . 5 ) 

ZZC+ZZC Al-NORMCDF ZC) 
a COMPUTE STUDENT T APPROXIMATION 

TT+TT.N SPAPROX 1 RRO 

a COMPUTE STUDENT T APPROXIMATION WITH CONTINUITY CORRECTION. 
TTC+TTC.N SPAPROX RHO 
*-(S<7)/L2 

a COMPUTE ERROR DIFFERENCES. 

ERRZZ+PP-ZZ 

ERRZZC+PP-ZZC 

ERRTT+PP-TT 

ERRTC+PP-TTC 

a PPIJVT OUT TABLE OF VALUES. 

I*-0 . 1 
c7*-19 5 

S*- ' TEST ST AT. VALUE ' 

M+M,' PROBZR 2 fi]; FCH SAJfPLE SIZE EQUAL TO 1 , (2 0 



*N ) , 



T . ffHVQ 

R. 



9N) , * 



1 1 2 



9999 >' UFMTd 7 pKK) 



5 ®FPCC] ) 



.39' 



+L 3 

L 5:«/«-194 

L 3 :M+M.ZI1 CU7C (18p 19 7 ) ,J, 61p (8pl97),«7] 
*(p=6 )/L4 
H+ 1 18 dS 
*(P*0)/L6 
F+'K*G<\ 9 
F* 1 62 pF 
FS+3.F 
+L7 

L6:FS+-o 0 
P«-17pS 
FS+FS , 3 
O 0 

L9:C«-C+1 
F+ 1 ]» (7 
FS*FS.F 
+(C<7 ) /L9 
L7 :M+M, [1] FS 
D+D + l 
1*1 

*(tfl f JV2 wV 3 ,A/4 JV5 ,L5 )CF] 
tfl:S«- r EXACr C .D .F . ‘ 

<7*1 98 
*£3 

N2'.S+' ERROR ; NORMAL ' 

PP+ERRZZ 

iV3:5^'FFF0F; NORM. W/CC 1 
PP+ERRZZC 
-*Z3 

N*:S+ l ERROR; T DIET » 

PP+ERRTZ 
*£3 

N5:S+ l ERROR; T W/CC ' 

PP+ERRTC 

->u3 

*OV<10)/il 
7 



Fij 

*3* 

¥ 

> 7 : 

* 8* 
10 ] 

14 , 

15 

16. 

17, 

183 

19 

20 , 

[ 21 . 

122, 

[23 

24 ] 
25] 
.26 
1 2 7 ; 
28] 
29] 
301 

*3 - 

■*o?- 

[34 

[35 

36 

.37 

t 3 8 

C 3 9 
40 
,41 
, 4 2 
, 4 3 
[44 

r 45 

[46 



7 WILTEST;N; ALPHA; B:A;C:TiP:NUM: DEN ;Z;ZZ:ERRZZ;M;ZZC;DiPP;NUMC; I;B;F;FS:S 
;J;TT:T; TTC ; TC ; NUMC ; NUMT ; DENOM ; DENOMC ; ERRA VE ; ERRTT ; ERRTC ; ERRZZC ; ZC ; ERRT 
CZCiKK 

ft THIS PROGRAM GENERATES TABLES OF C.D.F. COMPARISONS FOR TEE WILCOXON 
R SIGNED-RANK TEST. 



AN- 8 
ALP3A < 



SET SAMPLE SIZE AND ALPS A VALUES. 



COMPUTE CUMULATIVE DIST AND ASSOCIATED STATS. 

TEIS LOOP CALCULATES ALP3A VALUES AND APPORXIMATIONS . 



,01 0.02 0.03 0.05 0.08 0.13 0.18 
fi TEIS LOOP INCREMENTS SAMPLE SIZE. 

LllPP+pQ 
0 

FF*pO 
FS* pO 
ZZ*pO 
ZZOpO 

rr*po 
rrc*po 

N+N+l 
D+Q 
B+0 

*P+WILP N 

L2 : B*S+1 
C+ALPBACB 1 
A*- ( +/ (PSC) ) 

T+A-l 
KK+KK , T 
PP+PP,PZA 3 

R COMPUTE NORMAL APPROXIMATION . 

NUM+T-((Nx(N+l))**) 

DEN*-( ( N * (N+l )x(( 2 */V)+l)) + 24)*o.5 

Z+NUM*DEN 

ZZ+-ZZ t NORMS DF Z 

fl COMPUTE NORMAL APPROXIMATION WITS CONTINUITY CORRECTION . 

NUMC+(T+ 0 . 5 )- ( (Nx (N+l ) )+4 ) 

ZC+NUMC+DEN 
ZZC+ZZC , NORMCDF ZC 

R COMPUTE STUDENT T APPROXIMATION 

DENOM'*- ( ( (Nx (DEN* 2 ) )* (N-l ) ) - ( (NUM* 2 )+ (iV-l ) ) )*0 . 5 
T+NUM* DENOM 
TT+TT AN-1) TDISTN T 

R COMPUTE STUDENT T APPROXIMATION WITS CONTINUITY CORRECTION . 
NUMT+(NUM)+0.5 

DENOMC+(((Nx(DEN*2))*(N-l))-((((NUM)+Q.S)*2)*(N-l)))*Q.$ 

TC+NUMT* DENOMC 
TTC+TTCAiN- 1) TDISTN TC) 
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*(B < 7 )/L 2 

a COMPUTE ERROR DIFFERENCES. 

ERRZZ+PP-ZZ 
ERRZZC+PP-ZZC 
ERRTT*PP-TT 
ERRTC+PP-TTC x s 

ERRAVE*PP-( (ZZ+TT )+2 ) , 

ERP.TCZC+PP- 1 (ZZC+TTC 5 + 2 ) 
fi PRINT OUT TABLE OF VALUES. 

I* 0.1 

U+19S 

S +- ' TEST ST AT. VALUE » 

M+M,' PROBZW S W1 FOR SAMPLE SIZE EQUAL TO 

*L 3 

L 5 :« 7-<-194 

L3 :M*-M. CJ] Qfl7C(18pl97 ),«7,61p(8pl97),«7] 

-*•(£>=8 )/ L 4 
1 18 dS 
■(.D-0 )/L6 



A 



F* ' BC < 1 
F+ 1 62 
FS*H,F 

+L7 

L 6 :FS«-pO 
S-*-17 pS 
FS+FS.B 
C* 0 

L9:OC+l 



ZZ 9 
PF 



(7 5 



>' OFOTCl 7 pFF) 



) 

•tc< 7 }/L9 
JC7:^W,CI3 FS 
£>«-£+ 1 
J«-l 

-►(M .JY3 J75 f tf6 ,tf7 . JG5 ) CD] 
/n:s^ r Fxicr c.b.F. f 

*7-^198 

2 ^S-*- 1 FFF0F ; * 

PP+ERRZZ 

+L3 

N2:S+' ERROR; NORM. W/CC' 
PP+ERRZZC 
-+L2 

m iS+' ERROR; T DIST f 

PP+ERRTT 
+L 3 

N5:S+' ERROR; T W/CC * 

PP+ERRTC 
+L 3 

J76:S«- ! FFF0i?; A7F T/2 * 

PP+FFFA7F 
■+L3 

N7 :S+' ERROR; AVE TC/ZC 1 
PP+FFFTCZC 
+L 3 
£4:M 

+(tf<10)/£l 

V 



2 0 
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